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Express 



FedEx Express 
Customer Support Trace 
3875 Airways Boulevard 
Module H. 4th Floor 
Memphis, TN 381 16 



U.S. Ma* PO Box 727 
Memphis, TN 38194-4643 

Telephone: 901-369-3600 



11/19/2002 



Dear Customer 

Here is the proof of delivery for the shipment with tracking number 819623119682. Our records 
reflect the following information. 




Deliver y Information: 



Signed For By: A.BARRY 



re? 




Delivery Location: 16 PRESCOTT DR 
Delivery Date: June 19, 2002 
Delivery Time: 1103 



Shipping Information: 



Tracking No: 819623119682 



Ship Date: June 18, 2002 



Recipient: 

PAUL F CANTRELL 

15 PRESCOTT DR 
CHELMSFORD, MA 01824 
US 



Shipper: 

DAVID J POWSNER 
NUTTER MCCLENNEN & FISH LLP 
1 INTERNATIONAL PL FL 15 
BOSTON, MA 021102699 



Shipment Reference Information: 



0102323-00087 



Thank you for choosing FedEx Express. We look forward to working with you in the future. 

FedEx Worldwide Customer Service 
1-800-Go-FedEx® 

Reference No.: R20021 11900064455270 



Exhibit A 



http://www. fedex . com/cgi-bin/spod 



11/19/2002 



HI Nutter 




David J. Powsner 
Direct Line: 617-439-2717 
Fax: 617-310-9717 
E-mail: djp@nutter.com 



MERCURY COMPUTER SYSTEMS CONFIDENTIAL 



June 17, 2002 



Paul E. Cantrell 
15 Prescott Drive 
Chelmsford, MA 01863 

Re U.S. Patent Applications for 

WIRELESS COMMUNICATIONS SYSTEMS WIRELESS COMMUNICATIONS SYSTEMS 



AND METHODS FOR MULTIPLE 

PROCESSOR BASED MULTIPLE USER 

DETECTION 

Serial No.: 10/099,883 

Filing Date: 3/14/02 

OurRef.: 0102323-00099 



AND METHODS FOR CONTIGUOUSLY 
ADDRESSABLE MEMORY ENABLED 
MULTIPLE PROCESSOR BASED 
MULTIPLE USER DETECTION 
Serial No.: 10/099,912 
Filing Date: 3/14/02 
OurRef.: 0102323-00100 



WIRELESS COMMUNICATIONS SYSTEMS 
AND METHODS FOR CACHE ENABLED 
MULTIPLE PROCESSOR BASED 
MULTIPLE USER DETECTION 
Serial No.: 10/099,879 
Filing Date: 3/14/02 
OurRef.: 0102323-00101 



WIRELESS COMMUNICATIONS SYSTEMS 

AND METHODS FOR MULTIPLE 

OPERATING SYSTEM MULTIPLE USER 

DETECTION 

Serial No.: 10/099,918 

Filing Date: 3/14/02 

OurRef.: 0102323-00103 



WIRELESS COMMUNICATIONS SYSTEMS 
AND METHODS FOR NONVOLATILE 
STORAGE OF OPERATING PARAMETERS 
FOR MULTIPLE PROCESSOR BASED 
MULTIPLE USER DETECTION 
Serial No.: 10/099,545 
Filing Date: 3/14/02 
OurRef.: 0102323-00102 

WIRELESS COMMUNICATIONS SYSTEMS 

AND METHODS FOR DIRECT MEMORY 

ACCESS AND BUFFERING OF DIGITAL 

SIGNALS FOR MULTIPLE USER 

DETECTION 

Serial No.: 10/099,889 

Filing Date: 3/14/02 

OurRef.: 0102323-00104 



Dear Paul: 



Per my telephone message, I would appreciate your signing the enclosed patent applications on 
which you are named as an inventor by Mercury Computer Systems, Inc. These cover aspects of the 
MUD wireless telecommunications system you worked on while at Mercury. As with other Mercury 
proprietary information, please treat the applications as confidential. 



Nutter McClennen & Fish llp » Attorneys at Law 
One International Place ■ Boston, MA 021 1 0-2699 ■ 617-439-2000 ■ Fax:617-310-9000 ■ www.nutter.com 




Paul E. Cantrell 
June 17, 2002 
Page 2 



Once you have reviewed the applications for completeness and accuracy, please sign and date 
the "Combined Declaration and Power of Attorney" attached to the front of each of the applications, as 
well as the accompanying Assignments. Please sign the Assignment before a Notary Public . 

It is best to refrain from marking the patent applications, apart from signing each of the 
attached "Combined Declaration and Power of Attorney". Accordingly, if there are any points that you 
feel should be raised with respect to the content of the applications, please contact me. 

Also, please do not detach the "Combined Declaration and Power of Attorney" from each 
application. 

Once the papers have been executed, please use the enclosed Fedex envelopes to return them to 
me for filing with the U.S. Patent & Trademark Office. If you have any questions, please do not 
hesitate to give me a call. 



Sincerely yours, 




David J. Powsner 



cc: Rebecca Dowse 

DJP:crm 

Enclosures 



1115834.1 




^DECLARATION AND POWER OF ATTORNEY FOR 
' UNITED STATES LETTERS PATENT APPLICATION 



Attorney Docket No.: 0102323-00100 



As a below-named inventor, I hereby declare that: 

My residence, post-office address and citizenship are as stated below next to my name. 

I believe I am the original, first and sole inventor (if only one name is listed below) or an 
original, first and joint inventor (if plural names are listed below) of the subject matter which 
is claimed and for which a patent is sought on the invention entitled: 

WIRELESS COMMUNICATIONS SYSTEMS AND METHODS FOR 

CONTIGUOUSLY ADDRESSABLE MEMORY ENABLED 
MULTIPLE PROCESSOR BASED MULTIPLE USER DETECTION 



the specification of which 
(check one) 

EI is attached hereto. 

□ was filed on: 3/14/02 

as Application No.: 10/099,912 

and was amended on: 
(if applicable). 

In the event that the filing date and/or Application No. are not entered above at the time I 
execute this document, and if such information is deemed necessary, I hereby authorize and 
request my attorneys/agent(s) at Nutter, McClennen & Fish, LLP, One International Place, 
Boston, MA 02110-2699, to insert above the filing date and/or Application No. of said 
application. 

I hereby state that I have reviewed and understand the contents of the above-identified 
application specification, including the claims, as amended by any amendment specifically 
referred to herein. 

I acknowledge the duty to disclose all information known to me that is material to patentability 
in accordance with Title 37, Code of Federal Regulations, §1.56. 



FOREIGN PRIORITY CLAIM 

I hereby claim foreign priority benefits under Title 35, United States Code §119(a)-(d) of any 
foreign application(s) for patent or inventor's certificate listed below and have also identified 
below any foreign application for patent or inventor's certificate having a filing date before that 
of the application on which priority is claimed: 

(check one) 

0 no such foreign applications have been filed. 

such foreign applications have been filed as follows: 

EARLIEST FOREIGN APPLICATION(S), IF ANY FILED WITHIN 12 MONTHS 
(6 MONTHS FOR DESIGN) PRIOR TO THIS U,S. APPLICATION 



Country 


Application Number 


Date of Filing 
(month, day, year) 


Priority Claimed 
Under 35 USC 119 








Yes No , 








Yes No 








Yes No 








Yes No 








Yes No 



ALL FOREIGN APPLICATION(S), IF ANY FILED MORE THAN 12 MONTHS 
(6 MONTHS FOR DESIGN) PRIOR TO THIS U.S. APPLICATION 
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CLAIM FOR BENEFIT OF EARLIER U.S. PROVISIONAL APPLICATION(s) 

I hereby claim priority benefits under Title 35, United States Code §1 19(e), of any United 
States provisional patent application(s) listed below: 

(check one) 

□: no such U.S. provisional applications have been filed. 



0 such U.S. provisional applications have been filed as follows: 



Application Number 


Date of Filing 
(month, day, year) 


Priority Claimed 
Under 35 USC 119(e) 


60/275,846 


3/14/01 


X Yes No 


60/289,600 


5/7/01 


X Yes No 


60/295,060 


6/1/01 


X Yes No 



CLAIM FOR BENEFIT OF EARLIER U.S./PCT APPLICATION(s) 

I hereby claim the benefit under Title 35, United States Code §120, of the United States 
application(s) listed below and, insofar as the subject matter of each of the claims of this 
application is not disclosed in the prior United States application in the manner provided by the 
first paragraph of Title 35, United State Code, §112, I acknowledge the duty to disclose all 
information that is material to patentability in accordance with Title 37, Code of Federal 
Regulations, §1.56, and which became available to me between the filing date of the prior 
application and the national or PCT international filing date of this application: 

(check one) 

0 no such U.S./PCT applications have been filed. 

□: such U.S./PCT applications have been filed as follows: 



Application Number 



Date of Filing 
(month, day, year) 



Status 

(Patented/Pending/Abandoned) 



I hereby declare that all statements made herein of my own knowledge are true and that all 
statements made on information and belief are believed to be true; and further that these 
statements were made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment or both under Section 1001 of Title 18 of the United States 
Code and that such willful false statements may jeopardize the validity of the application or any 
patent issued thereon. 



I hereby appoint: 

Ronald E. Cahill 
Thomas J. Engellenner 
William C. Geary III 
Lisa J. Michaud 
Reza Mollaaghababa 



Reg. No. 38,403 

Reg. No. 28,711 

Reg. No. 31,359 

Reg. No. 44,238 

Reg. No. 43,810 



Herbert A. Newborn 
Tram Anh T. Nguyen 
David J. Powsner 
Jasbir Sagoo 



Reg. No. 42,031 

Reg. No. 47,257 

Reg. No. 31,868 

Reg. No. 51,177 



all of Nutter McClennen & Fish, LLP, One International Place, Boston, Massachusetts 02110- 
2699, jointly, and each of them severally, my attorneys at law/patent agent(s), with full power 
of substitution, delegation and revocation, to prosecute this application, to make alterations and 
amendments therein, to receive the patent, and to transact all business in the Patent and 
Trademark Office connected therewith. 

Please mail correspondence to: David J. Powsner 
at Customer Number 021J25, whose address is: 

Nutter McClennen & Fish, LLP 
One International Place 
Boston, Massachusetts 02110-2699 

Please direct telephone calls to: David J. Powsner 
at (617)439-2717. 



Please direct facsimiles to: (617)310-9717 



Full name of sole or first joint inventor 




John H. Oates 




Inventor's Signature 


Date 


Residence 




59B Seaverns Bridge Road, Amherst, New Hampshire 03031 




Country of Citizenship 




US 




Post Office Address (required) 




same as above 





Full name of second joint inventor 




Steven R. Imperiali 




Inventor's Signature 


Date 


Residence 




43 Haynes Road, Townsend, Massachusetts 01469 




Country of Citizenship 




US 




Post Office Address (required) 




same as above 





Full name of third joint inventor 




Alden J. Fuchs 




Inventor's Signature 


Date 


Residence 




160 Pine Hill Road, Nashua, New Hampshire 03063 




Country of Citizenship 




US 




Post Office Address (required) 




same as above 





Full name of fourth joint inventor 




Kathleen J. Jacques 




Inventor's Signature 


Date 


Residence 




10 Juniper Street, Jay, Maine 04239 




Country of Citizenship 




US 




Post Office Address (required) 




same as above 





Full name of fifth joint inventor 




Jonathan E. Greene 




Inventor's Signature 


Date 


Residence 




83 Hollenbeck Avenue, Great Barrington, Massachusetts 01230 




Country of Citizenship 




US 




Post Office Address (required) 




same as above 





Full name of sixth joint inventor 




William J. Jenkins 




Inventor's Signature 


Date 


Residence 




479 Whitepond Road, Leominster, Massachusetts 01453 




Country of Citizenship 




US 




Post Office Address (required) 




same as above 





Full name of seventh joint inventor 




Frank P. Lauginiger 




Inventor's Signature 


Date 


Residence 




772 Russell Station Road, Francestown, New Hampshire 03043 




Country of Citizenship 




US 




Post Office Address (required) 




same as above 





Full name of eighth joint inventor 




David E. Majchrzak 




Inventor's Signature 


Date 


Residence 




33 Jessica Drive, Merrimack, New Hampshire 03054 




Country of Citizenship 




US 




Post Office Address (required) 




same as above. 





Full name of ninth joint inventor 




Paul E. Cantrell 




Inventor's Signature 


Date 


Residence . 




15 Prescott Drive, Chelmsford, Massachusetts 01863 




Country of Citizenship 




US 




Post Office Address (required) 




same as above 






Full name of tenth joint inventor 




Mirza Cifric 




Inventor's Signature 


Date 


Residence 




180 Arborway, Jamaica Plain, Massachusetts 02130 




Country of Citizenship 




US 




Post Office Address (required) 




same as above 





Full name* of eleventh inint inventor 

A Uil JlUlIlv v/ 1 viv v vi ill 1 jv/lili ill » viiivi 




Ian N Dunn 




Inventor's Signature 


Date 


Residence 




48 North Court Street, #1, Providence, Rhode Island 02903 




Country of Citizenship 




US 




Post Office Address (required) 




same as above 





Full name of twelfth joint inventor 




Michael J. Vinskus 




Inventor's Signature 


Date 


Residence 




5 Cranberry Lane, Litchfield, New Hampshire 03052 




Country of Citizenship 




US 




Post Office Address (required) 




same as above 
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Attorney Docket: 0102323-00100 

ASSIGNMENT 

WHEREAS, we, John H. Oates, of 59B Seaverns Bridge Road, Amherst, New 
Hampshire 03031, Steven R. Imperiali, of 43 Haynes Road, Townsend, Massachusetts 
01469, Alden J. Fuchs, of 160 Pine Hill Road, Nashua, New Hampshire 03063, Kathleen 
J. Jacques, of 10 Juniper Street, Jay, Maine 04239, Jonathan E. Greene, of 83 
Hollenbeck Avenue, Great Barrington, Massachusetts 01230, William J. Jenkins, of 479 
Whitepond Road, Leominster, Massachusetts 01453, Frank P. Lauginiger, of 772 Russell 
Station Road, Francestown, New Hampshire 03043, David E. Majchrzak, of 33 Jessica 
Drive, Merrimack, New Hampshire 03054, Paul E. Cantrell, of 15 Prescott Drive, 
Chelmsford, Massachusetts 01863, Mirza Cifric, of 180 Arborway, Jamaica Plain, 
Massachusetts 02130, Ian N. Dunn, of 48 North Court Street, #1, Providence, Rhode 
Island, 02903, and Michael J. Vinskus, of Litchfield, New Hampshire 03052, having 
invented certain new and useful improvements in WIRELESS COMMUNICATIONS 
SYSTEMS AND METHODS FOR CONTIGUOUSLY ADDRESSABLE MEMORY 
ENABLED MULTIPLE PROCESSOR BASED MULTIPLE USER DETECTION, for which 
we have made application for United States Letters Patent, said application filed March 14, 
2002 and accorded U.S. Serial No. 10/099,912 and 

WHEREAS, Mercury Computer Systems, Inc., a Massachusetts corporation, having 
a place of business at 199 Riverneck Road, Chelmsford, Massachusetts 01824, and its 
successors and assigns (collectively hereinafter called "ASSIGNEE"), is desirous of acquiring 
our entire right, title and interest in and to the invention, and in and to the said applications and 
any Letters Patent that may issue thereon; 

NOW THEREFORE, for good and valuable consideration, the receipt and sufficiency 
of which is hereby acknowledge, we do hereby sell assign and transfer unto said ASSIGNEE, 
its successors, assigns and legal representatives the full and exclusive right, title and interest in 
and to said invention and in and to said applications and all patents which may granted 
therefore, and all divisions, reissues, substitutions, continuations, continuations-in-part and 
extensions thereof; and we hereby authorize and request the Commissioner of Patents and 



Trademarks to issue all patents for said invention, or patents resulting therefrom, insofar as our 
interest is concerned, to the said ASSIGNEE of our entire right, title and interest. 

We also hereby sell and assign to said ASSIGNEE, its successors, assigns and legal 
representatives the full and exclusive rights, title and interest to the invention disclosed in said 
applications throughout the world, including the right to file applications and obtain patents, 
utility models, industrial models and designs for said invention in its own name throughout the 
world, including all rights of priority, all rights to publish cautionary notices reserving 
ownership of said invention, and all rights to register said invention in appropriate registries; 
and we further agree to execute any and all powers of attorney, applications, assignments, 
declarations, affidavits, and any other papers in connection therewith necessary to perfect such 
rights, title and interest in ASSIGNEE, its successors, assigns and legal representatives. 

We hereby further agree that we will communicate to said ASSIGNEE, or to its 
successors, assigns and legal representative, any facts known to us respecting any 
improvements; and, at the expense of said ASSIGNEE, to testify in any legal proceedings, sign 
all lawful papers, execute all divisional, continuation, continuation-in-part, reissue and 
substitute applications, make all lawful oaths, and generally do everything possible to vest title 
in said ASSIGNEE and to aid said ASSIGNEE, its successors, assigns and legal 
representatives to obtain and enforce proper protection for said invention in all countries. 



in witness whereof, I hereto set my hand and seal at 

this day of , 20 . 

(City, State) (Day) (Month) (year) 



John H. Oates 

STATE OF 

SS. 

COUNTY OF 



Before me this day of , 20 , personally appeared John H. 

Oates known to me to be the person whose name is subscribed to the foregoing Assignment 
and acknowledged that he executed the same as his free act and deed for the purposes therein 
contained. 



NOTARY PUBLIC 

[Notary's Seal Here] My Commission Expires: 



IN WITNESS WHEREOF, I hereto set my hand and seal at 




this day of , 20 


(City, State) 


(Day) (Month) (year) 




Steven R. Imperiali 


STATE OF 






ss. 


COUNTY OF 




Before me this 


day of , 20 personally appeared Steven R. 


Imperiali known to me to be the person whose name is subscribed to the foregoing 


Assignment and acknowledged that he executed the same as his free act and deed for the 


purposes therein contained. 




NOTARY PUBLIC 


[Notary's Seal Here] 


My Commission Expires: 



in witness whereof, I hereto set my hand and seal at 




this day of , 20 


(City, State) 


(Day) (Month) (year) 


STATE OF 


Alden J. Fuchs 


COUNTY OF 


ss. 


Before me this day of , 20 personally appeared Alden J. 
Fuchs known to me to be the person whose name is subscribed to the foregoing Assignment 
and acknowledged that he executed the same as his free act and deed for the purposes therein 
contained. 




NOTARY PUBLIC 


[Notary's Seal Here] 


My Commission Expires: 



IN witness WHEREOF, I hereto set my hand and seal at 




this day of , 20 


(City, State) 


(Day) (Month) (year) 




Kathleen J. Jacques 


STATE OF 






SS. 


COUNTY OF 




Before me this 


day of , 20 personally appeared Kathleen J. 


Jacques known to me to be the person whose name is subscribed to the foregoing 


Assignment and acknowledged that he executed the same as his free act and deed for the 


purposes therein contained. 




NOTARY PUBLIC 


[Notary's Seal Here] 


My Commission Expires: 



in witness whereof, I hereto set my hand and seal at 




this day of ,20 


(City, State) 


(Day) (Month) (year) 


STATE OF 


Jonathan E. Greene 


COUNTY OF 


ss. 


Before me this _ day of .20 . personally appeared Jonathan F 
Greene known to me to be the person whose name is subscribed to the foregoing Assignment 
and acknowledged that he executed the same as his free act and deed for the purposes therein 
contained. 




NOTARY PUBLIC 


[Notary's Seal Here] 


My Commission Expires: 




IN WITNESS WHEREOF, 


I hereto set my hand and seal at 
this day of ,20 


(City, State) 


(Day) (Month) (year) 


STATE OF 


William J. Jenkins 


COUNTY OF 


SS. 


Before me this day of . 20 . personally app^reH William t 
Jenkins known to me to be the person whose name is subscribed to the foregoing 
Assignment and acknowledged that he executed the same as his free act and deed for the 
purposes therein contained. 




NOTARY PUBLIC 


[Notary's Seal Here] 


My Commission Expires: 



in witness whereof, I hereto set my hand and seal at 

this day of , 20 . 

(City, State) (Day) ~ (Month) (year) 



Frank P. Lauginiger 

STATE OF 

SS. 

COUNTY OF 



Before me this day of , 20 personally appeared Frank P. 

Lauginiger known to me to be the person whose name is subscribed to the foregoing 
Assignment and acknowledged that he executed the same as his free act and deed for the 
purposes therein contained. 



NOTARY PUBLIC 

[Notary's Seal Here] My Commission Expires: 



IN witness WHEREOF, I hereto set my hand and seal at 




this day of , 20 


(City, State) 


(Day) (Month) (year) 


STATE OF 


David E. Majchrzak 


COUNTY OF 


SS. 


Before me this 


da y of , 20 personally appeared David E. 


Majchrzak known to me to be the person whose name is subscribed to the foregoing 
Assignment and acknowledged that he executed the same as his free act and deed for the 
purposes therein contained. 




NOTARY PUBLIC 


[Notary's Seal Here] 


My Commission Expires: 



IN witness whereof, I hereto set my hand and seal at 




this day of ,20 


(City, State) 


(Day) (Month) (year) 


STATE OF 


Paul E. Cantrell 


COUNTY OF 


ss. 


Before me this 


day of , 20 , personally appeared Paul E. 


Cantrell known to me to be the person whose name is subscribed to the foregoing 
Assignment and acknowledged that he executed the same as his free act and deed for the 
purposes therein contained. 




NOTARY PUBLIC 


[Notary's Seal Here] 


My Commission Expires: 



IN 



witness WHEREOF, I hereto set my hand and seal at 



this 



(City, State) 



day of 



(Day) 



,20 



(Month) 



(year) 



STATE OF 



Mirza Cifric 

ss. 



COUNTY of 



Before me this 



day of 



20 , personally appeared Mirza 



Cifric known to me to be the person whose name is subscribed to the foregoing Assignment 
and acknowledged that he executed the same as his free act and deed for the purposes therein 
contained. 



I [Notary's Seal Here] 



NOTARY PUBLIC 
My Commission Expires: 
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IN witness whereof, I hereto set my hand and seal at 




this day of , 20 


(City, State) 


(Day) (Month) (year) 
Ian N. Dunn 


STATE OF 


ss. 


COUNTY OF 




Before me this day of , 20 personally appeared Ian N. 
Dunn known to me to be the person whose name is subscribed to the foregoing Assignment 
and acknowledged that he executed the same as his free act and deed for the purposes therein 
contained. 




NOTARY PUBLIC 


[Notary's Seal Here] 


My Commission Expires: 



IN witness WHEREOF, I hereto set my hand and seal at 

this day of _, 20 . 

(City, State) (Day) ~ (Month) (year) 



Michael J. Vinskus 

STATE OF 

SS. 

COUNTY OF 



Before me this day of , 20 , personally appeared Michael J. 

Vinskus known to me to be the person whose name is subscribed to the foregoing 
Assignment and acknowledged that he executed the same as his free act and deed for the 
purposes therein contained. 



NOTARY PUBLIC 

[Notary's Seal Here] My Commission Expires: 



U05366.1 



-8- 
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IN THE UNITED 



STATES PATENT AND TRADEMARK OFFICE 



United States Patent Application 



for 



JohnH. Oates 

59B Seavems Bridge Road 

Amherst, New Hampshire 03031 

Steven R. Imperial! 
43 Haynes Road 

Townsend, Massachusetts 0146V 

Alden J. Fuchs 

160 Pine Hill Road 

Nashua, New Hampshire 03063 

Kathleen J. Jacques 
10 Juniper Street 
Jay, Maine 04239 



Inventors ( 

FrankP.Lauginiger 
772 Russell Station Road 
Francestown, New Hampshire 03043 

David E.Majchrzak 

11 Regine Street m|M1 

Hudson, New Hampshire 03051 

Paul E. Cantrell 
15PrescottDrive 

Chelmsford, Massachusetts 01863 
MirzaCifric 

180Arborway n0 
Jamaica Plain, Massachusetts 02130 



Jonathan E. Greene 

83 Hollenbeck Avenue 

Seat Barrington, Massachusetts 01230 

William J. Jenkins 
479 Whitepond Road 
Leominster, Massachusetts 01453 



Ian N.Dunn 

48 North Court Street #1 
Providence, Rhode Island 02903 

Michael J. Vinskus 

5 Cranberry Lane 

Litchfield, New Hampshire 03052 
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Background of the Invention 

This application claims the benefit of priority of (i) US Provisional Application Serial 
No. 60/275,846 filed March 14, 2001, entitled "Improved Wireless Communications Systems 
5 and Methods"; (ii) US Provisional Application Serial No. 60/289,600 filed May 7, 2001, enti- 
tled "Improved Wireless Communications Systems and Methods Using Long-Code Multi-User 
Detection"' and (iii) US Provisional Application Serial Number. 60/295,060 filed June 1, 2001 
entitled "Improved Wireless Communications Systems and Methods for a Communications 
Computer," the teachings all of which are incorporated herein by reference. 

10 

The invention pertains to wireless communications and, more particularly, by way of 
example, to methods and apparatus providing multiple user detection for use in code division 
multiple access (CDMA) communications. The invention has application, by way of non-lim- 
iting example, in improving the capacity of cellular phone base stations. 

15 

Code-division multiple access (CDMA) is used increasingly in wireless communica- 
tions. It is a form of multiplexing communications, e.g., between cellular phones and base 
stations, based on distinct digital codes in the communication signals. This can be contrasted 
with other wireless protocols, such as frequency-division multiple access and time-division 
20 multiple access, in which multiplexing is based on the use of orthogonal frequency bands and 
orthogonal time-slots, respectively. 

A limiting factor in CDMA communication and, particularly, in so-called direct 
sequence CDMA (DS-CDMA), is interference — both that wrought on individual transmissions 
30 by buildings and other "environmental" factors, as well that between multiple simultaneous 
communications, e.g., multiple cellular phone users in the same geographic area using their 
phones at the same time. The latter is referred to as multiple access interference (MAI). Along 
with environmental interference, it has effect of limiting the capacity of cellular phone base 
stations, driving service quality below acceptable levels when there are too many users. 

35 

A technique known as multi-user detection (MUD) is intended to reduce multiple 
access interference and, as a consequence, increases base station capacity. It can reduce inter- 
ference not only between multiple transmissions of like strength, but also that caused by users 
so close to the base station as to otherwise overpower signals from other users (the so-called 
40 near/far problem). MUD generally functions on the principle that signals from multiple simul- 
taneous users can be jointly used to improve detection of the signal from any single user. Many 
forms of MUD are discussed in the literature; surveys are provided in Moshavi, "Multi-User 
Detection for DS-CDMA Systems," IEEE Communications Magazine (October, 1996) and 

1 



Express Mail Label: US EV 093 931 8 06 US 



Duel-Hallen et al, "Multiuser Detection for CDMA Systems," IEEE Personal Communications 
(April 1995). Though a promising solution to increasing the capacity of cellular phone base 
stations, MUD techniques are typically so computationally intensive as to limit practical appli- 
cation. 

5 

An object of this invention is to provide improved methods and apparatus for wireless 
communications. A related object is to provide such methods and apparatus for multi-user 
detection or interference cancellation in code-division multiple access communications. 

10 A further related object is to provide such methods and apparatus as provide improved 

short-code and/or long-code CDMA communications. 

A further object of the invention is to provide such methods and apparatus as can be 
cost-effectively implemented and as require minimal changes in existing wireless communica- 
1 5 tions infrastructure. 

A still further object of the invention is to provide methods and apparatus for executing 
multi-user detection and related algorithms in real-time. 

20 A still further object of the invention is to provide such methods and apparatus as 

manage faults for high-availability. 



30 



35 



40 



2 
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Summary of the Invention 

The foregoing and other objects are among those attained by the invention which pro- 
vides methods and apparatus for multiple user detection (MUD) processing. These have appli- 
5 cation, for example, in improving the capacity CDMA and other wireless base stations 

Wireless Communications Systems And Methods For 
Multiple Processor Based Multiple User Detection 

10 One aspect of the invention provides a multiuser communications device for detecting 

user transmitted symbols in CDMA short-code spread spectrum waveforms. A first processing 
element generates a matrix (hereinafter, "gamma matrix") that represents a correlation between 
a short-code associated with one user and those associated with one or more other users. A set 
of second processing elements generates, e.g., from the gamma matrix, a matrix (hereinafter, 

1 5 "R-matrix") that represents cross-correlations among user waveforms based on their ampli- 
tudes and time lags. A third processing element produces estimates of the user transmitted 
symbols as a function of the R-matrix. 

In related aspects, the invention provides a multiuser communications device in which 
20 a host controller performs a "partitioning function," assigning to each second processing ele- 
ment within the aforementioned set a portion of the R-matrix to generate. This partitioning can 
be a function of the number of users and the number of processing elements available in the set. 
According to related aspects of the invention, as users are added or removed from the spread 
spectrum system, the host controller performs further partitioning, assigning each second pro- 
30 cessing element within the set a new portion of the R-matrix to generate. 

Further related aspects of the invention provide a multiuser communications device as 
described above in which the host controller is coupled to the processing elements by way of a 
multi-port data switch. Still further related aspects of the invention provide such a device in 
35 which the first processing element transfers the gamma-matrix to the set of second processing 
elements via a memory element. 

Similarly, the set of second processing elements place the respective portions of the R- 
matrix in memory accessible to the third processing element via the data switch. Further 
40 related aspects of the invention provide devices as described above in which the host controller 
effects data flow synchronization between the first processing element and the set of second 
processing elements, as well as between the set of second processing elements and the third 
processing element. 

3 
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Wireless Communications Systems And Methods For 
Contiguously Addressable Memory Enabled Multiple 
Processor Based Multiple User Detection 

5 Another aspect of the invention provides a multiuser communications device for detect- 

ing user transmitted symbols in CDMA short-code spread spectrum waveforms in which a set 
of first processing elements generates a matrix (hereinafter the "R-matrix") that represents 
cross-correlations among user waveforms based on their amplitudes and time lags. The first 
processing elements store that matrix to contiguous locations of an associated memory. 

10 

Further aspects of the invention provide a device as described above in which a second 
processing element, which accesses the contiguously stored R-matrix, generates estimates of 
the user transmitted symbols. 

15 Still further aspects of the invention provide such a device in which a third processing 

element generates a further matrix (hereinafter, "gamma-matrix") that represents a correlation 
between a CDMA short-code associated with one user and those associated with one or more 
other users; this gamma-matrix used by the set of first processing elements in generating the 
R-matrix. In related aspects, the invention provides such a device in which the third process- 

20 ing element stores the gamma-matrix to contiguous locations of a further memory. 

In other aspects, the invention provides a multiuser device as described above in which 
a host controller performs a "partitioning function" of the type described above that assigning 
to each processing element within the set a portion of the R-matrix to generate. Still further 
30 aspects provide such a device in which the host controller is coupled to the processing elements 
by way of a multi-port data switch. 

Other aspects of the invention provide such a device in which the third processing ele- 
ment transfers the gamma-matrix to the set of first processing elements via a memory ele- 
35 ment. 

Further aspects of the invention provide a multiuser communications device as 
described above with a direct memory access (DMA) engine that places elements of the R- 
matrix into the aforementioned contiguous memory locations. 

40 

Further aspects of the invention provide methods for operating a multiuser communica- 
tions device paralleling the operations described above. 
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Wireless Communications Systems And Methods For 
Cache Enabled Multiple Processor Based 
Multiple User Detection 

5 Other aspects of the invention provide a multiuser communications device that makes 

novel use of cache and random access memory for detecting user transmitted symbols in 
CDMA short-code spectrum waveforms. According to one such aspect, there is provided a 
processing element having a cache memory and a random access memory. A host controller 
places in the cache memory data representative of characteristics of the user waveforms. The 

10 processing element generates a matrix as a function of the data stored in the cache, and stores 
the matrix in either the cache or the random access memory. 

Further aspects of the invention provide a device as described above in which the host 
controller stores in cache data representative of the user waveforms short-code sequences. The 
1 5 processing element generates the matrix as a function of that data, and stores the matrix in 
random access memory. 

Still further aspects of the invention provide such a device in which the host controller 
stores in cache data representative of a correlation of time-lags between the user waveforms 
20 and data representative of a correlation of complex amplitudes of the user waveforms. The 
host controller further stores in random access memory data representing a correlation of short- 
code sequences for the users waveforms. The processing element generates the matrix as a 
function of the data and stores that matrix in RAM. 

30 Further aspects of the invention provide a device as described above in which a host 

controller stores in cache an attribute representative of a user waveform, and stores in random 
access memory an attributes representing a cross-correlation among user waveforms based on 
time-lags and complex amplitudes. The processing element generates estimates of user trans- 
mitted symbols and stores those symbols in random access memory. 

35 

Other aspects of the invention provide such a device in which the host controller trans- 
mits the matrix stored in the cache or random access memory of a processing element to the 
cache or random access memory of a further processing element. 

40 Further aspects of the invention provide a multiuser communications device as 

described above with a multi-port data switch coupled to a short-code waveform receiver 
system and also coupled to a host controller. The host controller routes data generated by the 
receiver system to the processing element via the data switch. 

5 
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Further aspects of the invention provide methods for operating a multiuser communica- 
tions device paralleling the operations described above. 

Wireless Communications Systems And Methods For 
5 Nonvolatile Storage Of Operating Parameters For 

Multiple Processor Based Multiple User Detection 

Another aspect of the invention provides a multiuser communications device for detect- 
ing user transmitted symbols in CDMA short-code spectrum waveforms in which fault and 

10 configuration information is stored to a nonvolatile memory. A processing element, e.g. that 
performs symbol detection, is coupled with random access and nonvolatile memories. A fault 
monitor periodically polls the processing element to determine its operational status. If the 
processing element is non-operational, the fault monitor stores information including configu- 
ration and fault records, as well at least a portion of data from the processing element's RAM, 

1 5 into the nonvolatile memory. 

According to further aspects according to the invention, following detection of the non- 
operational status, the fault monitor sends to a host controller a reset-request interrupt together 
with the information stored in the nonvolatile RAM. In turn, the host controller selectively 
20 issues a reset command to the processing element. In related aspects, the processing element 
resets in response to the reset command and transfers (or copies) the data from the nonvolatile 
memory into the RAM, and therefrom continues processing the data in the normal course. 

Further aspects of the invention provide a device as described above in which the pro- 
30 cessing element periodically signals the fault monitor and, in response, the fault monitor polls 
the processing element. If the fault monitor does not receive such signaling within a specified 
time period, it sets the operational status of the processing element to non-operational. 

According to a related aspect of the invention, the fault monitor places the processing 
35 elements in a non-operational status while performing a reset. The fault monitor waits a time 
period to allow for normal resetting and subsequently polls the processor to determine its 
operational status. 

Still further aspects of the invention provide a device as described above in which there 
40 are a plurality of processing elements, each with a respective fault monitor. 

Yet still further related aspects of the invention provide for the fault monitoring a data 
bus coupled with the processing element. 
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Further aspects of the invention provide methods for operating a multiuser communica- 
tions device paralleling the operations described above. 

Wireless Communications Systems And Methods For 
5 Multiple Operating System Multiple User Detection 

Another aspect of the invention provides a multiuser communications device for detect- 
ing user transmitted symbols in CDMA short-code spectrum waveforms in which a first pro- 
cess operating under a first operating system executes a first set of communication tasks for 
10 detecting the user transmitted symbols and a second process operating under a second operat- 
ing system — that differs from the first operating system — executes a second set of tasks for like 
purpose. A protocol translator translates communications between the processes. According to 
one aspect of the invention, the first process generates instructions that determine how the 
translator performs such translation. 

15 

According to another aspect of the invention, the first process sends a set of instructions 
to the second process via the protocol translator. Those instructions define the set of tasks 
executed by the second process. 

20 In a related aspect of the invention, the first process sends to the second process instruc- 

tions for generating a matrix. That can be, for example, a matrix representing any of a correla- 
tion of short-code sequences for the user waveforms, a cross-correlation of the user waveforms 
based on time-lags and complex amplitudes, and estimates of user transmitted symbols embed- 
ded in the user waveforms. 

30 

Further aspects of the invention provide a device as described above in which the first 
process configures the second process, e.g., via data sent through the protocol translator. This 
can include, for example, sending a configuration map that defines where a matrix (or portion 
thereof) generated by the second process is stored or otherwise directed. 

35 

Still further aspects of the invention provide a device as described above in which the 
first process is coupled to a plurality of second processes via the protocol translator. Each of 
the latter processes can be configured and programmed by the first process to generate a respec- 
tive portion of a common matrix, e.g., of the type described above. Further aspects of the inven- 
40 tion provide methods for operating a multiuser communications device paralleling the 
operations described above. 
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Wireless Communications Systems And Methods For 
Direct Memory Access And Buffering Of 
Digital Signals For Multiple User Detection 

5 Another aspect of the invention provides a multiuser communications device for detect- 

ing user transmitted symbols in CDMA short-code spectrum waveforms in which a program- 
mable logic device (hereinafter "PLD") enables direct memory access of data stored in a digital 
signal processor (hereinafter "DSP"). The DSP has a memory coupled with a DMA controller 
that is programmed via a host port. The PLD programs the DMA controller via the host port to 

1 0 allow a buffer direct access to the memory. 

In a related aspect according to the invention, the PLD programs the DMA controller to 
provide non-fragmented block mode data transfers to the buffer. From the buffer, the PLD 
moves the blocks to a data switch that is coupled to processing devices. In a further related 
15 aspects according to the invention, the PLD programs the DMA controller to provide frag- 
mented block mode data transfers utilizing a protocol. The PLD provides the protocol which 
fragments and unfragments the blocks prior to moving them to the data switch. 

In further aspects provided by a device as described above, the PLD is implemented as 
20 a field programmable gate array that is programmed by a host controller coupled with the data 
switch. In a related aspect, the PLD is implemented as a application specific integrated circuit 
which is programmed during manufacture. In still aspects, a device as described above pro- 
vides for a buffer implemented as a set of registers, or as dual-ported random access memory. 

30 Further aspects of the invention provide methods for operating a multiuser communica- 

tions device paralleling the operations described above. 
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Improved Wireless Communications Systems And Methods 
For Short-code Multiple User Detection 

Still further aspects of the invention provide methods for processing short code spread 
5 spectrum waveforms transmitted by one or more users including the step of generating a matrix 
indicative of cross correlations among the waveforms as a composition of (i) a first component 
that represents correlations among time lags and short codes associated with the waveforms 
transmitted by the users, and (ii) a second component that represents correlations among 
multipath signal amplitudes associated with the waveforms transmitted by the users. The 
10 method further includes generating detection statistics corresponding to the symbols as a func- 
tion of the correlation matrix, and generating estimates of the symbols based on those detection 
statistics. 

Related aspects of the invention provided methods as described above in which the first 
1 5 component is updated on a time scale that is commensurate with a rate of change of the time 
lags associated with the transmitted waveforms, and the second component is updated on a dif- 
ferent time scale, i.e., one that is commensurate with a rate of change of the multipath ampli- 
tudes associated with these waveforms. In many embodiments, the updating of the second 
component, necessitated as a result of change in the multipath amplitudes, is executed on a 
20 shorter time scale than that of updating the first component. 

Other aspects of the invention provide methods as described above in which the first 
component of the cross-correlation matrix is generated as a composition of a first matrix com- 
ponent that is indicative of correlations among the short codes associated with the respective 
30 users, and a second matrix component that is indicative of the waveforms transmitted by the 
users and the time lags associated with those waveforms. 

In a related aspect, the invention provides methods as above in which the first matrix 
component is updated upon addition or removal of a user to the spread spectrum system. This 
35 first matrix component (referred to below as T-matrix) can be computed as a convolution of the 
short code sequence associated with each user with the short codes of other users. 

According to further aspects of the invention, elements of the T-matrix are computed in 
accord with the relation: 

40 , 

r„[m]s— J^c,[n]c t [n-m] 

n=0 

wherein 
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c i M represents the complex conjugate of a short code sequence associated with the P h 
user, 

c k [n-m] represents a short code sequence associated with the k h user, 

5 

N represents a length of the short code sequence, and 

N f represent a number of non-zero length of the short code sequence. 

^ In further aspects, the invention provides a method as described above in which the first 

component of the cross-correlation matrix (referred to below as the C matrix) is obtained as a 
function of the aforementioned T-matrix in accord with the relation: 



15 



20 



30 



m 

wherein 

g is a pulse shape vector, 

Nc is the number of samples per chip, 

x is a time lag, and 

T represents the T matrix, e.g., defined above. 

In a related aspect, the cross-correlation matrix (referred to below as the R-matrix) can 
be generated as a function of the C matrix in accord with the relation: 



35 



9=1 9 =1 I J 



40 



wherein 

a is an estimate of % , the complex conjugate of one multipath amplitude component 
of the 1 th user, 

a kq . is one multipath amplitude component associated with the k lh user, and 
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C denotes the C matrix, e.g., as defined above. 

In further aspects, the invention provides methods as described above in which the 
detection statistics are obtained as a function of the cross-correlation matrix (e.g., the R-matrix) 
5 in accord with the relation: 

K K K 

fl[m] = r f [Oft[m]+Xr.[-l^ 

*=1 *=1 *=1 



10 wherein 

y[m] represents a detection statistic corresponding to m* symbol transmitted by the P h 
user, 

1 5 ^/[0]^/[ m ] represents a signal of interest, and 

remaining terms of the relation represent Multiple Access Interference (MAI) and 
noise. 

20 In a related aspect, the invention provides methods as described above in which esti- 

mates of the symbols transmitted by the users and encoded in the short code spread spectrum 
waveforms are obtained based on the computed detection statistics by utilizing, for example, a 
multi-stage decision-feedback interference cancellation (MDFIC) method. Such a method can 
provide estimates of the symbols, for example, in accord with the relation: 

30 f K & % ~ 

Mm] = */^Mm]-£r.[-l]M^^ 

[ *=] A=1 *=l 

wherein 

35 b { [m] represents an estimate of the m th symbol transmitted by the 1 th user. 

Further aspects of the invention provide logic carrying out operations paralleling the 
methods described above. 

40 
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Load Balancing Computational Methods In A Short-code Spread-spectrum Communi- 
cations System 

In further aspects, the invention provides methods for computing the cross-correlation 
5 matrix described above by distributing among a plurality of logic units parallel tasks — each for 
computing a portion of the matrix. The distribution of tasks is preferably accomplished by 
partitioning the computation of the matrix such that the computational load is distributed sub- 
stantially equally among the logic units. 

10 In a related aspect, a metric is defined for each partition in accord with the relation 

below. The metric is utilized as a measure of the computational load associated with each logic 
unit to ensure that the computational load is distributed substantially equally among the logic 
units: 

15 B=A r A,-l 
wherein 

A. represents an area of a portion of the cross-correlation matrix corresponding 
20 to the /'* partition, and 

i represents an index corresponding to the number of logic units over which the 
computation is distributed. 

30 In another aspect, the invention provides methods as described above in which the 

cross-correlation matrix is represented as a composition of a rectangular component and a tri- 
angular component. Each area, represented by A. in the relation above, includes a first portion 
corresponding to the rectangular component and a second portion corresponding to the triangu- 
lar component. 

35 

Further aspects of the invention provide logic carrying out operations paralleling the 
methods described above. 

Hardware And Software For Performing Computations In A 
40 Short-code Spread-spectrum Communications System 

In other aspects, the invention provides an apparatus for efficiently computing a Y- 
matrix as described above, e.g., in hardware. The system includes two registers, one associated 
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with each of P h and k* users. The registers hold elements of the short code sequences associated 
with the respective user such that alignment of the short code sequence loaded in one register 
can be shifted relative to that of the other register by m elements. Associated with each of the 
foregoing registers is one additional register storing mask sequences. Each element in those 
5 sequences is zero if a corresponding element of the short code sequence of the associated reg- 
ister is zero and, otherwise, is non-zero. The mask sequences loaded in these further registers 
are shifted relative to the other by m elements. A logic performs an arithmetic operation on the 
short code and mask sequences to generate, for rrf h transmitted symbol, the (/, k) element of the 
T-matrix, i.e., TJm] 

10 

In a related aspect, the invention provides an apparatus as described above in which the 
arithmetic operation performed by the logic unit includes, for any two aligned elements of the 
short code sequences of the 1 th and k* user and the corresponding elements of the mask 
sequences, (i) an XOR operation between the short code elements, (ii) an AND operation 
1 5 between the mask elements, (iii) an AND operation between results of the step (i) and step (ii). 
The result of step (iii) is a multiplier for the aligned elements, which the logic sums in order to 
generate the (/, k) element of the T-matrix. 

Further aspects of the invention provide methods paralleling the operations described 

20 above. 



Improved Computational Methods For Use In A 
Short-code Spread-spectrum Communications System 



30 In still further aspects, the invention provides improved computational methods for 

calculating the aforesaid cross-correlation matrix by utilizing a symmetry property. Methods 
according to this aspect include computing a first one of two matrices that are related by a sym- 
metry property, and calculating a second one of the two matrices as a function of the first com- 
ponent through application of the symmetry property. 

35 

According to related aspects of the invention, the symmetry property is defined in 
accord with the relation: 

R lk (m) = ^R kJ (-m). 

40 

wherein 
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\(m) and R^m) refer to (/, k) and (k, I) elements of the cross-correlation matrix, 
respectively. 

Further aspects of the invention provide methods as described above in which calcula- 
5 tion of the cross-correlation matrix further includes determining a C matrix that represents 
correlations among time lags and short codes associated with the waveforms transmitted by the 
users, and an R-matrix that represents correlations among multipath signal amplitudes associ- 
ated with the waveforms transmitted by the users. In related aspects the step of determining the 
C matrix includes generating a first of two C-matrix components related by a symmetry prop- 
1 0 erty. A second of the components is then generated by applying the symmetry property. 

Related aspects of the invention provide a method as described above including the step 
of generating the r-matrix in accord with the relation: 

J am 

15 r tt [/w] = — rJTcjM c k [n-m] 

wherein 

c][n] represents complex conjugate of the short code sequence associated with the 1th 
20 user, 

c k [n-m\ represents the short code sequence associated with kth user, 

N represents the length of the code, and 

30 

N { represent the number of non-zero length of the code. 

Further aspects of the invention provide logic carrying out operations paralleling the 
methods described above. 

35 
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Wireless Communications Systems And Methods For Virtual User Based Multiple 
User Detection Utilizing Vector Processor Generated Mapped Cross-correlation 
Matrices 

5 Still further aspects of the invention provide methods for detecting symbols encoded in 

physical user waveforms, e.g., those attributable to cellular phones, modems and other CDMA 
signal sources, by decomposing each of those waveforms into one or more respective virtual 
user waveforms. Each waveform of this latter type represents at least a portion of a symbol 
encoded in the respective physical user waveforms and, for example, can be deemed to "trans- 

10 mit" a single bit per symbol period. Methods according to this aspect of the invention deter- 
mine cross-correlations among the virtual user waveforms as a function of one of more 
characteristics of the respective physical user waveforms. From those cross-correlations, the 
methods generate estimates of the symbols encoded in the physical user waveforms. 

1 5 Related aspects of the invention provide methods as described above in which a physi- 

cal user waveforms is decomposed into a virtual user waveform that represents one or more 
respective control or data bits of a symbol encoded in the respective physical user waveform. 

Other related aspects provide for generating the cross-correlations in the form of a first 
20 matrix, e.g., an R-matrix for the virtual user waveforms. That matrix can, according to still 
further related aspects of the invention, be used to generate a second matrix representing cross- 
correlations of the physical user waveforms. This second matrix is generated, in part, as a 
function of a vector indicating the mapping of virtual user waveforms to physical user wave- 
forms. 

30 

Further aspects of the invention provide a system for detecting symbols encoded in 
physical user waveforms that has multiple processors, e.g., each with an associated vector pro- 
cessor, that operates in accord with the foregoing methods to generate estimates of the symbols 
encoded in the physical user waveforms. 

35 

Still other aspects of the invention provide a system for detecting user transmitted sym- 
bols encoded in short-code spread spectrum waveforms that generates cross-correlations 
among the waveforms as a function of block-floating integer representations of one or more 
characteristics of those waveforms. Such a system, according to related aspects of the inven- 
40 tion, utilizes a central processing unit to form floating-point representations of virtual user 
waveform characteristics into block-floating integer representations. A vector processor, 
according to further related aspects, generates the cross-correlations from the latter representa- 
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Hons. The central processing unit can "reformat" the resulting block-floating point matrix into 
floating-point format, e.g., for use in generating symbol estimates. 

Still further aspects of the invention provide methods and apparatus employing any and 
5 all combinations of the foregoing. These and other aspects of the invention, which includes 
combinations of the foregoing, are evident in the illustrations and in the text that follows. 
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Brief Description of the Illustrated Embodiment 

A more complete understanding of the invention may be attained by reference to the 
drawings, in which: 

5 

Figure 1 is a block diagram of components of a wireless base-station utilizing a multi- 
user detection apparatus according to the invention; 

Figure 2 is a block diagram of components of a multiple user detection processing card 
1 0 according to the invention; 

Figure 3 is a more detailed view of the processing board of Figure 2; 

Figure 4 depicts a majority-voter sub-system in a system according to the invention; 

15 

Figure 5 is a block diagram of an integrated direct memory access (DMA) engine of the 
type used in a system according to the invention; 

Figures 6 and 7 depict power on/off curves for the processor board in a system accord- 
20 ing to the invention; 

Figure 8 are an operational overview of functionality within the host processor and 
multiple compute nodes in a system according to the invention; 

30 Figure 9 is a block diagram of an external digital signal processor apparatus used to 

supply digital signals to the processor board in a system according to the invention; 

Figure 10 illustrates an example of loading the R matrices on multiple compute nodes 
in a system according to the invention; 

35 

Figure 11 depicts a short-code loading implementation with parallel processing of the 
matrices in a system according to the invention; 

Figure 12 depicts a long-code loading implementation utilizing pipelined processing 
40 and a triple-iteration of refinement in a system according to the invention; 

Figure 13 illustrates skewing of multiple user waveforms; 

17 
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Figure 14 is a graph illustrating MUD efficiency as a function of user velocity in units 
ofKm/hr. 

Figure 15 schematically illustrates a method for defining a common interval for three 
5 short-code streams utilized in a FFT calculation of the T-matrix; 

Figure 16 schematically illustrates the T-matrix elements calculated upon addition of a 
new physical user to a system according to the invention; 

1 0 Figures 17,18 and 1 9 depict hardware calculation of the T-matrix in a system according 

to the invention; 

Figure 20 illustrates parallel computation of the R and C matrices in a system according 
to the invention; 

15 

Figure 21 depicts a use of a vector processor using integer operands for generating a 
cross-correlation matrix of virtual user waveforms in a system according to the invention. 

20 
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Detailed Description of the Illustrated Embodiment 

Code-division multiple access (CDMA) waveforms or signals transmitted, e.g., from a 
user cellular phone, modem or other CDMA signal source, can become distorted by, and 

5 undergo amplitude fades and phase shifts due to phenomena such as scattering, diffraction 
and/or reflection off buildings and other natural and man-made structures. This includes 
CDMA, DS/CDMA, IS-95 CDMA, CDMAOne, CDMA2000 IX, CDMA2000 lxEV-DO, 
WCDMA (or UTMS), and other forms of CDMA, which are collectively referred to hereinafter 
as CDMA or WCDMA. Often the user or other source (collectively, "user") is also moving, 

10 e.g., in a car or train, adding to the resulting signal distortion by alternately increasing and 
decreasing the distances to and numbers of building, structures and other distorting factors 
between the user and the base station. 

In general, because each user signal can be distorted several different ways en route to 
15 the base station or other receiver (hereinafter, collectively, "base station"), the signal may be 
received in several components, each with a different time lag or phase shift. To maximize 
detection of a given user signal across multiple tag lags, a rake receiver is utilized. Such a 
receiver is coupled to one or more RF antennas (which serve as a collection point(s) for the 
time-lagged components) and includes multiple fingers, each designed to detect a different 
20 multipath component of the user signal. By combining the components, e.g., in power or 
amplitude, the receiver permits the original waveform to be discerned more readily, e.g., by 
downstream elements in the base station and/or communications path. 

A base station must typically handle multiple user signals, and detect and differentiate 
30 among signals received from multiple simultaneous users, e.g., multiple cell phone users in the 
vicinity of the base station. Detection is typically accomplished through use of multiple rake 
receivers, one dedicated to each user. This strategy is referred to as single user detection 
(SUD). Alternately, one larger receiver can be assigned to demodulate the totality of users 
jointly. This strategy is referred to as multiple user detection (MUD). Multiple user detection 
35 can be accomplished through various techniques which aim to discern the individual user sig- 
nals and to reduce signal outage probability or bit-error rates (BER) to acceptable levels. 

However, the process has heretofore been limited due to computational complexities 
which can increase exponentially with respect to the number of simultaneous users. Described 
40 below are embodiments that overcome this, providing, for example, methods for multiple user 
detection wherein the computational complexity is linear with respect to the number of users 
and providing, by way of further example, apparatus for implementing those and other methods 
that improve the throughput of CDMA and other spread-spectrum receivers. The illustrated 
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embodiments are implemented in connection with short-code CDMA transmitting and receiver 
apparatus; however those skilled in the art will appreciate that the methods and apparatus 
therein may be used in connection with long-code and other CDMA signalling protocols and 
receiving apparatus, as well as with other spread spectrum signalling protocols and receiving 
5 apparatus. In these regards and as used herein, the terms long-code and short-code are used in 
their conventional sense: the former referring to codes that exceed one symbol period; the 
latter, to codes that are a single symbol period or less. 

Figure 1 depicts components of a wireless base station 100 of the type in which the 
10 invention is practiced. The base station 100 includes an antenna array 114, radio frequency/ 
intermediate frequency (RF/IF) analog-to-digital converter (ADC), multi-antenna receivers 
110, rake modems 112, MUD processing logic 118 and symbol rate processing logic 120, 
coupled as shown. 

1 5 Antenna array 1 14 and receivers 1 10 are conventional such devices of the type used in 

wireless base stations to receive wideband CDMA (hereinafter "WCDMA") transmissions 
from multiple simultaneous users (here, identified by numbers 1 through K). Each RF/IF 
receiver (e.g., 1 10) is coupled to antenna or antennas 1 14 in the conventional manner known in 
the art, with one RF/IF receiver 110 allocated for each antenna 114. Moreover, the antennas 

20 are arranged per convention to receive components of the respective user waveforms along dif- 
ferent lagged signal paths discussed above. Though only three antennas 114 and three receiv- 
ers 110 are shown, the methods and systems taught herein may be used with any number of 
such devices, regardless of whether configured as a base station, a mobile unit or otherwise. 
Moreover, as noted above, they may be applied in processing other CDMA and wireless com- 

30 munications signals. 

Each RF/IF receiver 110 routes digital data to each modem 112. Because there are 
multiple antennas, here, Q of them, there are typically Q separate channel signals communi- 
cated to each modem card 112. 

35 

Generally, each user generating a WCDMA signal (or other subject wireless communi- 
cation signal) received and processed by the base station is assigned a unique short-code code 
sequence for purposes of differentiating between the multiple user waveforms received at the 
basestation, and each user is assigned a unique rake modem 1 12 for purposes of demodulating 
40 the user's received signal. Each modem 1 1 2 may be independent, or may share resources from, 
a pool. The rake modems 1 1 2 process the received signal components along fingers, with each 
receiver discerning the signals associated with that receiver's respective user codes. The 
received signal components are denoted here as denoting the channel signal (or wave- 

20 
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form) from the fc 01 user from the 9 th antenna, or r k [t] denoting all channel signals (or wave- 
forms) originating from the 4 th user, in which case r k [t] is understood to be a column vector 
with one element for each of the Q antennas. The modems 112 process the received signals 
r k D] to generate detection statistics y^im] for the A* user for the mth symbol period. To this 

5 end, the modems 122 can, for example, combine the components ^[t] by power, amplitude or 
otherwise, in the conventional manner to generate the respective detection statistics y\ 
In the course of such processing, each modem 1 12 determines the amplitude (denoted herein as 
a ) of and time lag (denoted herein as t) between the multiple components of the respective 
user channel. The modems 112 can be constructed and operated in the conventional manner 

'0 known in the art, optionally, as modified in accord with the teachings of some of the embodi- 
ments below. 

The modems 112 route their respective user detection statistics y^lm] , as well as the 
amplitudes and time lags, to common user detection (MUD) 1 1 8 logic constructed and oper- 

' 5 ated as described in the sections that follow. The MUD logic 1 1 8 processes the received signals 
from each modem 1 1 2 to generate a refined output, y { k ]) [m] , or more generally, y { k n) [m] , where 
n is an index reflecting the number of times the detection statistics are iteratively or regenera- 
tively processed by the logic 118. Thus, whereas the detection statistic produced by the 
modems is denoted as J^ 0) [w] indicating that there has been no refinement, those generated by 

20 processing the y* 0) [«] detection statistics with logic 1 18 are denoted , those generated 

by processing the y^l***] detection statistics with logic 1 18 are denoted ^[m] , and so forth. 
Further waveforms used and generated by logic 118 are similarly denoted, e.g., r (n) [/]. 

Though discussed below are embodiments in which the logic 1 1 8 is utilized only once, 
^® i.e., to generate ^l 0 !^] from j4 0) [w] , other embodiments may employ that logic 118 multiple 
times to generate still more refined detection statistics, e.g., for wireless communications appli- 
cations requiring lower bit error rates (BER). For example, in some implementations, a single 
logic stage 1 18 is used for voice applications, whereas two or more logic stages are used for 
data applications. Where multiple stages are employed, each may be carried out using the same 
35 hardware device (e.g., processor, co-processor or field programmable gate array) or with a suc- 
cessive series of such devices. 

The refined user detection statistics, e.g., .y^t"*] or more generally y\ n) [m], are com- 
municated by the MUD process 118 to a symbol process 120. This determines the digital 
^ information contained within the detection statistics, and processes (or otherwise directs) that 
information according to the type of user class for which the user belongs, e.g., voice or data 
user, all in the conventional manner. 
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Though the discussion herein focuses on use of MUD logic 1 18 in a wireless base sta- 
tion, those skilled in the art will appreciate that the teachings hereof are equally applicable to 
MUD detection in any other CDMA signal processing environment such as, by way of non- 
limiting example, cellular phones and modems. For convenience, such cellular base stations 
5 other environments are referred to herein as "base stations." Multiple User Detection Process- 
ing Board 

Figure 2 depicts a multiple user detection (MUD) processing card according to the 
invention. The illustrated processing card 118 includes a host processor 202, an interface block 
10 204, parallel processors 208, a front panel device 210, and a multi-channel cross-over device 
206 (hereinafter "Crossbar"). Although these components are shown as separate entities, one 
skilled in the art can appreciate that different configurations are possible within the spirit of the 
invention. For example, the host processor 202 and the interface block 204 can be integrated 
into a single assemble, or multiple assemblies. 

15 

The processing card 1 1 8 processes waveform and waveform components received by 
a base station, e.g., from a modem card 1 12 or receiver 110 contained within the base station, 
or otherwise coupled with the base station. The waveform typically includes CDMA wave- 
forms, however the processing card 118 can also be configured for other protocols, such as 
20 TDMA and other multiple user communication techniques. The processing card 1 1 8 performs 
multiple user detection (MUD) on the waveform data, and generates a user signal correspond- 
ing to each user, with includes less interference than within the received signals.. 

The illustrated processing card 1 1 8 is a single board assembly and is manufactured to 
30 couple (e.g., electrically and physically mate) with a conventional base station (e.g., a modem 
card 1 12, receiver 1 10 or other component). The board assembly illustrated conforms to a 3 A 
form factor modem payload card of the type available in the marketplace. The processor card 
1 1 8 is designed for retrofitting into existing base stations or for design into new station equip- 
ment. In other embodiments, the processing card can be either single or multiple assemblies. 

35 

The host processor 202 routes data from the interface block 204 to and among the paral- 
lel processors 208, as well as performs fault monitoring and automated resets, data transfer, and 
processor loading of the parallel processors 208. The host processor 202 also processes output 
received from the parallel processors 208, and communicates the processed output to the inter- 
40 face block 204 for subsequent return to the base station. 

The parallel processors 202 process waveforms and waveform components routed from 
the host processor 206. Typically, the parallel processors 202 process the waveform compo- 
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nents, and communicate the processed data back to the host processor 202 for further process- 
ing and subsequent transmission to the base station, however, the intermediate processed 
waveforms can be communicated to other parallel processors or directly to the base station. 

5 The crossbar 206 is a communication switch which routes messages between multiple 

devices. It allows multiple connection data ports to be connection with other data ports. In the 
illustrated embodiment, the crossbar 206 provides eight ports, where a port can be "connected" 
to any other port (or to multiple ports) to provide communication between those two (or indeed, 
multiple) ports. Here, the crossbar 206 is a RACEway™ switch of the type commercially 

10 available from the assignee hereof. In other embodiments, other switching elements, whether 
utilizing the RACEway™ protocol or otherwise, may be used, e.g., PCI, I2C and so on. Indeed, 
in some embodiments, the components communicate along a common bus and/or are distrib- 
uted via over a network. 

15 A front panel 210 is used to monitor the processor card and can be used to apply soft- 

ware patches, as well as perform other maintenance operations. Additionally, the front panel 
2 1 0 can be used to monitor fault status and interface connections through a series of LED indi- 
cators, or other indicators. Illustrated front panel interfaces with the board via the RACEway™ 
switch and protocol, though other interface techniques may be used as well. 

20 

Figure 3 depicts further details of the processor card of Figure 2. The illustrated pro- 
cessor card includes a host processor 202 in communication with an interface block 205 and a 
set of parallel processors 208 (hereinafter "compute elements") as described above, as well as 
a crossbar 206 and a front panel 210. Further, a power status/control device 240 is assembled 
30 on the processor card 118. However, in other embodiments, the power status/control device 
240 can be within the base station or elsewhere. 

The host processor 202 includes a host controller 203 with an integrated processor con- 
taining a peripheral logic block and a 32-bit processor core. The host controller 203 is coupled 
35 with various memory devices 205, a real time clock 206, and a protocol translator 208. In the 
illustrated embodiment, the host controller 203 can be a Motorola PowerPC 8240 commer- 
cially available, but it will be appreciated by one skilled in the art that other integrated proces- 
sors (or even non-integrated processors) can be used which satisfy the requirements herein. 

40 The host controller 203 controls data movement within the processor card 118 and 

between the processor card and the base station. It controls the crossbar device 206 by assign- 
ing the connection between connection ports. Further, the host controller 203 applies function- 
ality to the output generated by the parallel processors 208. The host controller 203 includes a 
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The multiple-port communication device 206 (hereinafter "crossbar") provides com- 
munication between all processing and input/output elements on the processing card 118. In 
the illustrated embodiment, the crossbar 206 is an EEPROM device which can be read and 
programmed by a RACEway™ compatible component (e.g., the front panel 210 or parallel 
5 processors 208), but it is typically programmed initially during manufacture. An embedded 
ASIC device controls the EEPROM programming, and hence, the function of the crossbar 
206. 

The crossbar 206 in the illustrated provides up to three simultaneous 266-Mbytes/Sec 
10 throughput data paths between elements for a total throughput of 798 Mbytes/Sec, however, in 
other embodiments the actual throughput varies according to processing speed. Here, two 
crossbar ports (e.g., ports 0 and 1) connect to a bridge FPGA which further connect to the front 
panel 210. Each of the multiple processors use an crossbar port (e.g., ports 2, 3, 5, and 6), and 
the interface block 224 and host controller 203 share one crossbar port (e.g., port 4) via the 
15 PXB 206. The number of ports on the crossbar 206 depends on the number of parallel proces- 
sors and other components that are in communication. 

The multiple processors 208 in the illustrated embodiment include four compute ele- 
ments 220a-220d (hereinafter, reference to element 220 refers to a general compute element, 

20 also referred to herein as a "processing element" or "CE"). Each processing element 220 
applies functionality on data, and generates processed date in the form of a matrix, vector, or 
waveform. The processing elements 220 can also generate scalar intermediate values. Gener- 
ated data is passed to the host controller 208, or to other processing elements 220 for further 
processing. Further, individual processing elements can be partitioned to operate in series 

30 (e.g., as a pipeline) or in parallel with the other processing elements. 

A processing element 220 includes a processor 228 coupled with a cache 230, a Joint 
Test Action Group (hereinafter "JTAG") interface 232 with an integrated programming port, 
and an application specific integrated circuit 234 (hereinafter "ASIC"). Further, the ASIC 234 
35 is coupled with a 1 28 Mbyte SDRAM device 236 and HA Registers 238. The HA Registers are 
coupled with 8 Kbytes of NVRAM 244. In the illustrated embodiment the compute elements 
220 are on the same assembly as the host controller 203. In other embodiments, the compute 
nodes 220 can be separate from the host controller 203 depending on the physical and electrical 
characteristics of the target base station. 

40 

The compute node processors 228 illustrated are Motorola PowerPC 7400, however in 
other embodiments the processor can be other processor devices. Each processor 228 uses the 
ASIC 234 to interface with a RACEway™ bus 246. The ASIC 234 provides certain features 
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of a compute node 220, e.g., a DMA engine, mail box interrupts, timers, page mapping regis- 
ters, SDRAM interface and the like. In the illustrated embodiment the ASIC is programmed 
during manufacture, however, it can also be programmed in the field, or even at system reset in 
other embodiments. 

5 

The cache 230 for each compute node 220 stores matrices that are slow-changing or 
otherwise static in relation to other matrices. The cache 230 is pipelined, single-cycle deselect, 
synchronous burst static random access memory, although in other embodiments high-speed 
RAM or similar devices can be used. The cache 230 can be implemented using various devices, 
10 e.g., multiple 64 Kbyte devices, multiple 256 Kbyte devices, and so on. 

Architecture Pairing of Processing Nodes with NVRAM and Watchdog; Majority Voter 

The HA registers 238 store fault status for the software and/or hardware of the compute 
15 element 220. As such, it responds to the watchdog fault monitor which also monitors the host 
controller 203 and other components. The NVRAM 244 is, much like the NVRAM coupled 
with the host controller 203, stores data from which the current state of the compute element 
220 can be recreated should a fault or reset occur. The SDRAM 236 is used for intermediate 
and temporary data storage, and is directly addressable from both the ASIC 234 and the proces- 
20 sor 228. These memory devices can be other devices in other embodiments, depending on 
speed requirements, throughput and computational complexity of the multiple user detection 
algorithms. 

NVRAM is also used to store computational variables and data such that upon reset of 
30 the processing element or host controller, execution can be re-started without the need to re- 
fresh the data. Further, the contents of NVRAM can be used to diagnose fault states and/or 
conditions, thus aiding to a determination of the cause of fault state. 

As noted above, a "watchdog" monitors performance of the processing card 118. In the 
35 illustrated embodiment, there are five independent "watchdog" monitors on the processing card 
118 (e.g., one for the host controller 203 and one each for each compute node 220a-220d, and 
so on). The watchdog also monitors performance of the PCI bus as well as the RaceWay bus 
connected with each processing element and the data switch. The RACEWay bus includes out- 
of-band fault management coupled with the watchdogs. 

40 

Each component periodically strobes its watchdog at least every 20 msec but not faster 
that 500 microseconds (these timing parameters vary among embodiments depending on over- 
all throughput of the components and clock speed). The watchdog is initially strobed approxi- 
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mately two seconds after the initialization of a board level reset, which allows for start-up 
sequencing of the components without cycling erroneous resets. Strobing the watchdog for the 
processing nodes is accomplished by writing a zero or a one sequence to a discrete word (e.g., 
within the HA Register 212) originating within each compute element 220a-220d, the host 
5 controller 203, and other components). The watchdog for the host controller 203 is serviced by 
writing to the memory mapped discrete location FFF D027 which is contained within the HA 
Registers 212. 

The watchdog uses five 8-bit status registers within the HA registers 2 1 2, and additional 
10 registers (e.g., HA registers 238) within each compute node 220. One register represents the 
host controller 203 status, and the other four represent each compute node 220a-220d status. 
Each register has a format as follows: 



Bit 


Name 


Description | 


0 


CHECKSTOP_OUT 


Checkstop state of CPU (0 = CPU in checkstop) 


1 


WDMFAULT 


WDM failed (0 = WDM failed, set high after reset and valid service) | 


2 


SOFTWARE_FAULT 


Software fault detected (Set to 0 when a software exception was 
detected) (R/W local) 


3 


RESETREQIN 


Wrap status of the local CPU's reset request 


4 


WDM INIT 


WDM failed in initial 2 second window ( 0 = WDM failed) 


5 


Software definable 0 


Software definable 0 


6 


Software definable I 


Software definable 1 


7 


Unused 


Unused 



The five registers reflect status information for all processors within the processing 
board 1 1 8, and allow the host controller 203 to obtain status of each without the need for poll- 
ing the processor individually (which would degrade performance and throughput). Addition- 
ally, the host controller 203 and each compute node processor 228 has a fault control register 
which contains fault data according to the following format: 



35 



40 



Bit 


Name 


Description 


0 


RESETREQOUT0 


Request a reset event (0 => forces reset) 


1 


CHKSTOPOUT0 


Request that node 0 enter checkstop state (0 => request checkstop) 


2 


CHKSTOPOUT I 


Request that node 1 enter checkstop state (0 => request checkstop) 


3 


CHKSTOPOUT_2 


Request that node 2 enter checkstop state (0 => request checkstop) 


4 


CHKSTOPOUT_3 


Request that node 3 enter checkstop state (0 => request checkstop) 


5 


CHKSTOPOUT_8240 


Request that the host controller enter checkstop state (0 => request 
checkstop) 


6 


Software definable 0 


Software definable 0 


7 


Software definable 1 


Software definable 1 
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A single write of any value will strobe the watchdog. Upon events such as power-up, 
the watchdogs are initialized to a fault state. Once a valid strobe is issued, the watchdog exe- 
cutes and, if all elements are properly operating, writes a no- fault state to the HA register 212. 
This occurs within the initial two-second period after board level reset. If a processor node fails 
5 to service the watchdog within the valid time frame, the watchdog records a fault state. A 
watchdog of a compute node 220 in fault triggers an interrupt to the host controller 203. If a 
fault is within the host controller 203, then the watchdog triggers a reset to the board. The 
watchdog then remains in a latched failed state until a CPU reset occurs followed by a valid 
service sequence. 

10 

Each processor node ASIC 234 accesses a DIAG3 signal that is wired to an HA register, 
and is used to strobe the compute element's hardware watchdog monitor. A D1AG2 signal is 
wired to the host processor's embedded programmable interrupt controller (EPIC) and is used 
by a compute element to generate a general purpose interrupt to the host controller 203. 

15 

A majority voter (hereinafter "voter") is a dual software sub-system state machine that 
identifies faults within each of the processors (e.g., the host controller 230 and each compute 
node 220a-220d) and also of the processor board 1 1 8 itself. The local voter can reset individual 
processors (e.g., a compute node 220) by asserting a CHECKSTOPJN to that processor. The 
20 board level voter can force a reset of the board by asserting a master reset, wherein all proces- 
sors are reset. Both voters follow a rule set that the output will follow the majority of non- 
checkstopped processors. If there are more processors in a fault condition than a non-fault 
condition, the voter will force a board reset. Of course, other embodiments may use other 
rules, or can use a single sub-system to accomplish the same purpose. 

30 

A majority voter is illustrated in Figure 4. Board level resets are initiated from a vari- 
ety of sources. One such source is a voltage supervisor (e.g., the power status/control device 
240) which can generate a 200 ms reset if the voltage (e.g., VCC) rises above a predetermined 
threshold, such as 4.38 volts (this is also used in the illustrated embodiment in a pushbutton 

35 reset switch 406, however, the push button can also be a separate signal). The board level voter 
will continue to drive a RESETJ) 408 until both the voltage supervisor 404 and the PCI_ 
RESET 0 410 are de-asserted. Either reset will generate the signal RESETJ) 412 which resets 
the card into a power-on state. RESETJ) 412 also generates HRESETJ) 414 and TRST 416 
signals to each processor. Further, a HRESETJ) and TRSTcan be generated by the JTAG ports 

40 using a JTAG HRESET 0 4 1 8 and JTAGTRST 420 respectively. The host controller 203 can 
generate a reset request, a soft reset (C_SRESET_0 422) to each processor, a check-stop 
request, and an ASIC reset (CE^RESET O 424) to each of the four compute element's ASIC. 
A discrete word from the 5v-powered reset PLD will generate the signal NPORESET 1 (not a 
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power on reset). This signal is fed into the host processor discrete input word. The host proces- 
sor will read this signal as logic low only if it is coming out of reset due to either a power con- 
dition or an external reset from off board. Each compute element, as well as the host processor 
can request a board level reset. These requests are majority voted, and the result RESET- 
VOTE_0 will generate a board level reset. 

Each compute node processor 228 has a hard reset signal driven by three sources gated 
together: a HRESETO pin 426 on each ASIC, a HRESETJ) 418 from the JTAG connector 
232, and a HRESETJ) 412 from the majority voter. The HRESETJ) pin 426 from the ASIC 
is set by the "node run" bit field (bit 0) of the ASIC MisconA register. Setting HRESETJ) 
426 low causes the node processor to be held in reset. HRESETJ) 426 is low immediately after 
system reset or power-up, the node processor is held in reset until the HRESETJ) line is pulled 
high by setting the node run bit to 1 . The JTAG HRESETJ) 41 8 is controlled by software when 
a JTAG debugger module is connected to the card. The HRESETJ) 412 from the majority voter 
is generated by a majority vote from all healthy nodes to reset. 

When a processor reset is asserted, the compute processor 228 is put into reset state. 
The compute processor 228 remains in a reset state until the RUN bit 0 of the Miscon A regis- 
ter is set to 1 and the host processor has released the reset signals in the discrete output word. 
The RUN bit is set to 1 after the boot code has been loaded into the SDRAM starting at location 
0x0000_0100. The ASIC maps the reset vector 0xFFF0_0100 generated by the MPC7400 to 
address OxOOOOJHOO. 

Turning now to discuss memory devices 205 coupled with the host controller 203, the 
memory devices are addressable by the host controller 203 as follows. The host controller 203 
addresses the memory devices (e.g., the HA registers 212, SDRAM 214, Flash 216 and 
NVRAM 218) using two address mapping configurations designated as address map A and 
address map B, although other configurations are possible. Address map A conforms to the 
PowerPC reference platform (PreP) specification (however, if other host controllers are used, 
map A conforms with a native reference platform to that host controller). Address map B con- 
forms to the host controller 203 common hardware reference platform (CHRP). 

Support of map A is provided for backward compatibility, and further supports any 
retrofitting of existing base station configurations. The address space of map B is divided into 
four areas: system memory, PCI memory, PCI Input/Output (I/O), and system ROM space. 
When configured for map B, the host controller translates addresses across the internal periph- 
eral logic bus and the external PCI bus as follows: 
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5 



Processor Core Address Range 






PCI Address Range 


Definition 


Hex 












0000_0000 


0009_FFFF 


0 


64 UK. - I 


NO PCI CYCLE 


System memory 


000A_0000 


000F_FFFF 


640JC 


1M-I 


0O0A 0000 - 000F FFFF 


Compatibility hole 


0010_0000 


3FFFFFFF 


IM 


1G-1 


NO PCI CYCLE 


System memory 


4000,0000 


7FFFFFFF 


1G 


2G-1 


via r»/~"i (~*\.t i—* i r* 

NO PCI CYCLE 


Reserved 


8000_0000 


FCFFFFFF 


2G 


4G-48M-1 


8000_0000 - FCFF FFFF 


PCI memory 


FD00 0000 


FDFFFFFF 


4G-48M 


4G-32M-1 


0000_0000 - 00FFFFFF 


PCI/ISA memory 


FEOO 0000 


FE7FFFFF 


4G-32M 


4G-24M-1 


0000_0000 - 007FFFFF 


PCI/ISA I/O 


FE80 0000 


FEBFFFFF 


4G-24M 


4G-20M-1 


0080 0000 - O0BF_FFFF 


PCI I/O 


FECO 0000 


FEDFFFFF 


4G-20M 


4G-18M-I 


CONFIG_ADDR 


PCI configuration address 


FEE0_0000 


FEEFFFFF 


4G-18M 


4G-I7M-1 


CONFIG_DATA 


PCI configuration data \ 


FEFO 0000 


FEFF FFFF 


4G-17M 


4G-I6M-I 


FEFO 0000 - FEFF FFFF 


PCI interrupt acknowledge 


FFOO 0000 


FF7F_FFFF 


4G-16M 


4G-8M-1 


FF0O 0OO0 - FF7F FFFF 


32/64-bit Flash/ROM 


FF80 0000 


FFFFFFFF 


4G-8M 


4G-I 


FF80_0000 - FFFF FFFF 


8/32/64-bit Flash/ROM 



15 In the illustrated embodiment, hex address FF00 0000 through FF7F FFFF is not used, 

and hence, that bank of Flash ROM is not used. The address of FF80 0000 through FFFF_ 
FFFF is used, as the Flash ROM is configured in 8-bit mode and is addressed as follows: 

20 
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35 



40 
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Bank Select 


Processor Core Address Range 


Definition 


11111 








MHO- 


FFE0 0000 


FFEFFFFF 


Accesses Bank 0 


onooi 

vUUvl 


FFE0_0000 


FFEFFFFF 


Application code (30 pages) 


00000 


rrcu__uuuu 


rrcr_rrrr 


Annlication/boot roHp 




FFFO_0000 


FFFF_CFFF 


Application/boot code 




FFFF_D000 


FFFFD000 


Discrete input word 0 




FFFFD00 1 


FFFFJDOOI 


Discrete input word 1 




FFFFD002 


FFFFD002 


Discrete output word 0 




FFFFD003 


FFFFD003 


Discrete output word 1 




FFFF D004 


FFFF D004 


Discrete output word 2 




FFFFD010 


FFFFD010 


IC (Pending interrupt) 




FFFF D011 


FFFF D01 1 


IC (Interrupt mask low) 




FFFFD012 


FFFFD012 


IC (Interrupt clear low) 




FFFFD013 


FFFF D013 


IC (Unmasked, pending low) 




FFFFD014 


FFFFJD014 


IC (Interrupt input low) 


xxxx 


FFFFD015 


FFFF_D015 


Unused (read FF) 


FFFF D016 


FFFF DO 16 


Unused (read FF) 




FFFFD017 


FFFFD017 


Unused (read FF) 




FFFFD018 


FFFFD018 


Unused (read FF) 




FFFFD019 


FFFFD019 


Unused (read FF) 




FFFF D020 


FFFFD020 . 


HA (Local HA register) 




FFFFD021 


FFFF D021 


HA (Node 0 HA register) 




FFFF D022 


FFFFD022 


HA (Node 1 HA register) 




FFFF D023 


FFFFD023 


HA (Node 2 HA register) 




FFFFD024 


FFFFD024 


HA (Node 3 HA register) 




FFFFD025 


FFFFD025 


HA (8240 HA register) 




FFFFD026 


FFFFD026 


HA (Software Fail) 




FFFFD027 


FFFFD027 


HA (Watchdog Strobe) 




FFFFD028 


FFFFDFFF 


4068 Bytes Flash 




FFFFE000 


FFFFFFFF 


8K NVRAM 



Address FFEF 0000 through FFEF FFFF contains 30 pages, and is used for 
application and boot code, as selected by the Flash bank bits. Further, there a 2 Mbyte block 
available after reset. Data movement occurs on the PCI 211a and/or a memory bus. 



DMA Engine Supported by Host Controller andFPGA 

40 

Direct memory access (DMA) is performed by the host controller 203, and operates 
independently from the host processor 203 core, as illustrated in Figure 5. The host controller 
203 has an integrated DMA engine including a DMA command stack 502, a DMA state engine 
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504, an address decode block 506, and three FIFO interfaces 508, 510, 512. The DMA engine 
receives and sends information via the PXB 208 coupled with the crossbar 206. 

The command stack 502 and state machine 504 processes DMA requests and transfers. 
5 The stack 502 and state machine 504 can initiate both cycle stealing and burst mode, along with 
host controller interupts. The address decode 506 sets the bus address, and triggers transmis- 
sions of the data. 

The host controller 203 has two DMA I/O interfaces, each with a 64-byte queue to 
1 0 facilitate the gathering and sending of data. Both the local processor and PCI masters can initi- 
ate a DMA transfer. The DMA controller supports memory transfers between PCI to memory, 
between local and PCI memory, and between local memory devices. Further, the host control- 
ler 203 can transfer in either block mode or scatter mode within discontinuous memory. A 
receiving channel 510 buffers data that is to be received by the memory. A transmit channel 
15 512 buffers data that is sent from memory. Of course, the buffers can also send/receive infor- 
mation from other devices, e.g., the compute nodes 220, or other devices capable of DMA 
transfers. 

The host controller 203 contains an embedded programmable interrupt controller 
20 (EPIC) device. The interrupt controller implements the necessary functions to provide a flex- 
ible and general-purpose interrupt controller. Further, the interrupt controller can pool inter- 
rupts generated from the several external components (e.g., the compute elements), and deliver 
them to the processor core in a prioritized manner. In the illustrated embodiment, an OpenPIC 
architecture is used, although it can be appreciated by one skilled in the art that other such 
30 methods and techniques can be used. Here, the host controller 203 supports up to five external 
interrupts, four internal logic-driven interrupts, and four timers with interrupts. 

Data transfers can also take effect via the FPGA program interface 508. This interface 
can program and/or accept data from various FPGAs, e.g., the compute note ASIC 234, cross- 
35 bar 242, and other devices. Data transfers within the compute node processor 228 to its ASIC 
234 and RACEway™ bus 246 are addressed as follows: 



40 
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5 



15 



From Address 


To Address 


Function 


0x0000 0000 


OxOFFF FFFF 


Local SDRAM 256 Mb 


0x1000 0000 


OxlFFFFFFF 


crossbar 256 Mb map window 1 


0x2000 oooo ; 


0x2FFF FFFF 


crossbar 256 MB map window 2 


0x3000 0000 


0x3FFF FFFF 


crossbar 256 MB map window 3 


0x4000 0000 


0x4FFF FFFF 


crossbar 256 MB map window 4 


0x5000 0000 


0x5 FFF FFFF 


crossbar 256 MB map window 5 


0x6000 0000 


0x6FFF FFFF 


crossbar 256 MB map window 6 


0x7000 0000 


0x7FFFFFFF 


crossbar 256 MB map window 7 


0x8000 0000 


0x8FFF FFFF 


crossbar 256 MB map window 8 


0x9000 0000 


0x9FFF FFFF 


crossbar 256 MB map window 9 


OxAOOO 0000 


OxAFFFFFFF 


crossbar 256 MB map window A 


OxBOOO 0000 


OxBFFF FFFF 


crossbar 256 MB map window B 


OxCOOO 0000 


OxCFFF FFFF 


crossbar 256 MB map window C 


[JXUVvv VUVU 


UxJJrrr rrrr 


crossoar zjo ivjd map winaow u 


OxEOOO OOOO 


OxEFFF FFFF 


crossbar 256 MB map window E \ 


OxFOOO 0000 


OxFBFF FBFF 


Not used (CE reg replicated mapping) 


OxFBFF FCOO 


OxFBFF FDFF 


Internal CN ASIC registers 


OxFBFF FEOO 


OxFEFF FFFF 


Pre-fetch control 


OxFFOO 0000 


OxFFFF FFFF 


16 MB boot FLASH memory area 



The SDRAM 236 can be addressable in 8, 16, 32 or 64 bit addresses. The RACEway™ 
bus 246 supports locked read/write and locked read transactions for all data sizes. A 16 Mbyte 
boot flash area is further divided as follows: 
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From Address 


To Address 


Function 


OxFFOO 2006 


OxFFOO 2006 


Software Fail Register 


OxFFOO 2005 


OxFFOO 2005 


MPC8240 HA Register 


OxFFOO 2004 


OxFFOO 2004 


Node 3 HA Register 


OxFFOO 2003 


OxFFOO 2003 


Node 2 HA Register 


OxFFOO 2002 


OxFFOO 2002 


Node 1 HA Register 


OxFFOO 2001 


OxFFOO 2001 


Node 0 HA Register 


OxFFOO 2000 


OxFFOO 2000 


Local HA Register (status/control) 


OxFFOO 0000 


OxFFOO 1FFF 


NVRAM 



40 



Slave accesses are accesses initiated by an external RACEway™ device directed 
toward the compute element processor 238. The ASIC 234 supports a 256 Mbyte address 
space which can be partitioned as follows: 



From Address 


To Address 


Function 


0x0000 0000 


OxOFFF FBFF 


256 MB less 1 Kb hole SDRAM 


0XfTT_FC00 


0xFFF_FFFF 


PCE133 internal registers 
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There are 16 discrete output signals directly controllable and readable by the host con- 
troller 203. The 1 6 discrete output signals are divided into two addressable 8-bit words. Writing 
to a discrete output register will cause the upper 8-bits of the data bus to be written to the dis- 
crete output latch. Reading a discrete output register will drive the 8-bit discrete output onto the 
5 upper 8-bits of the host processor data bus. The bits in the discrete output word are defined as 
follows: 

There are 16 discrete input signals accessible by the host controller 203. Reads from 
the discrete input address space will latch the state of the signals, and return the latched state of 
10 the discrete input signals to the host processor. The bits in the discrete input word are as fol- 
lows: 







Output Word 2 


DH(0:7) 


Signal 


Description 


0 


ND0_FLASH_EN_1 


Enable the CE ASIC's FLASH port when 1 


1 


ND1_FLASH_EN_1 


Enable the CE ASIC's FLASH port when 1 


2 


ND2 FLASH EN 1 


Enable the CE ASIC's FLASH port when 1 


3 


ND3_FLASH_EN1 


Enable the CE ASIC's FLASH port when 1 


4 


Wrap 1 


Wrap to discrete input 


5 






6 






7 











Output Word 1 


DH(0:7) 


Signal 


Description 


0 


WRAP0 


Wrap to Discrete Input 


1 


I2C_RESETJ) 


Reset the I2C serial bus when 0 


2 


SWLED 


Software controlled LED 


3 


FLASHSEL4 


Flash bank select address bit 4 


4 


FLASHSEL3 


Flash bank select address bit 3 


5 


FLASHSEL2 


Flash bank select address bit 2 


6 


FLASHSEL1 


Flash bank select address bit 1 


7 


FLASHSEL0 


Flash bank select address bit 0 
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uutput wora u 


DH(0:7) 


Signal 


Description 


0 


C_SRESET3_0 


Issue a Soft Reset to CPU on Node 3 when 0 


1 


C_PRESET3_0 


Reset PCE133 ASIC Node 3 when 0 


2 


C_SRESET2_0 


Issue a Soft Reset to cpu on Node 2 when 0 


3 


C_PRESET2_0 


Reset PCE133 ASIC Node 2 when 0 


4 


C_SRESET1_0 


Issue a Soft Reset to cpu on Node 1 when 0 


5 


CJ>RESET1_0 


Reset PCE133 ASIC Node 1 when 0 


6 


C_SRESET0_0 


Issue a Soft Reset to cpu on Node 0 when 0 


7 


C_PRESETOJ) 


Reset PCE133 ASIC Node 0 when 0 



10 







Input Word 1 


DH(0:7) 


Signal 


Description 


0 


WRAP1 


Wrap from discrete output word 


1 






2 


V3.3_FAIL_0 


Latched status of power supply since last reset 


3 


V2.5_FAIL_0 


Latched status of power supply since last reset 


4 


VCORE1_FAIL_0 


Latched status of power supply since last reset 


5 


VCORE0_FAIL_0 


Latched status of power supply since last reset 


6 


RIOR_CNF_DONE_1 


RIO/RACE++ FPGA configuration complete 


7 


PXB0_CNF_DONE_1 


PXB++ FPGA configuration complete 



20 
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Input Word 0 


DH(0:7) 


Signal 


Description 


0 


WRAP0 


Wrap from discrete output word 


1 


WDMSTATUS 


MPC8240's watchdog monitor status (0 = failed) 


2 


NPORESETJ 


Not a power on reset when high 


3 






4 






5 






6 






7 







The host controller 203 interfaces with an 8-input interrupt controller external from 
processor itself (although in other embodiments it can be contained within the processor). The 
interrupt inputs are wired, through the controller to interrupt zero of the host processor external 
^ interrupt inputs. The remaining four host processor interrupt inputs are unused. 

The Interrupt Controller comprises the following five 8-bit registers: 
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Resister 


Description 


Pending Register 


A low bit indicates a falling edge was detected on that interrupt 
(read only); 


Clear Register 


Setting a bit low will clear the corresponding latched interrupt 
(write only); 


Mask Register 


Setting a bit low will mask the pending interrupt from generating 
a processor interrupt; 


Unmasked Pending 
Register 


A low bit indicates a pending interrupt that is not masked out 


Interrupt State Register 


indicates the actual logic level of each interrupt input pin. 



The interrupt input sources and their bit positions within each of the six registers are as 
follows: 



15 



Bit 


Signal 


Description 


0 


SWFAILJ) 


8240 Software Controlled Fail Discrete 


1 


RTCJNTJ) 


Real time clock event 


2 


NODE0_FAIL_0 


WDFAILJ) or IWDFAILJ) or SWFAILJ) active 


3 


NODE1_FAIL_0 


WDFAILJ) or IWDFAILJ) or SWFAILJ) active 


4 


NODE2 FAIL 0 


WDFAILJ) or IWDFAILJ) or SWFAILJ) active 


5 


NODE3_FAIL_0 


WDFAILJ) or IWDFAILJ) or SWFAILJ) active 


6 


PCIJNTJ) 


PCI interrupt 


7 


XB_SYS_ERR_0 


crossbar internal error 



A falling edge on an interrupt input will set the appropriate bit in the pending register 
low. The pending register is gated with the mask register and any unmasked pending interrupts 
will activate the interrupt output signal to the host processor external interrupt input pin. Soft- 
ware will then read the unmasked pending register to determine which interrupt(s) caused the 
exception. Software can then clear the interrupt(s) by writing a zero to the corresponding bit in 
the clear register. If multiple interrupts are pending, the software has the option of either servic- 
ing all pending interrupts at once and then clearing the pending register or servicing the highest 
priority interrupt (software priority scheme) and the clearing that single interrupt. If more inter- 
rupts are still latched, the interrupt controller will generate a second interrupt to the host pro- 
cessor for software to service. This will continue until all interrupts have been serviced. 

An interrupt that is masked will show up in the pending register but not in the unmasked 
pending register and will not generate a processor interrupt. If the mask is then cleared, that 
pending interrupt will flow through the unmasked pending register and generate a processor 
interrupt. 
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ramp down for the 3.3V supply 704a occurs before the ramp down for the 2.5V supply 704c, 
which occurs before the ramp down for the 1.8V supplies 704c. 

Also, The host processor requires the core supply to not exceed the I/O supply by more 
5 than 0.4 volts at all times. Also, the I/O supply must not exceed the core supply by more than 
2 volts. Therefore, to achieve an acceptable power-up and power-down sequencing, e.g., to 
avoid damage to the components, a circuit containing diodes is used in conjunction with the 
power supplied within the base station. 

10 The power status/control device 240 is designed from a programmable logic device 

(PLD). The PLD is used to monitor the voltage status signals from the on board supplies. It is 
powered up from +5V and monitors +3.3V, +2.5V, 1 .8V_1 and +1 .8V 2. This device monitors 
the power jgood signals from each supply. In the case of a power failure in one or more sup- 
plies, the PLD will issue a restart to all supplies and a board level reset to the processor board. 

15 A latched power status signal will be available from each supply as part of the discrete input 
word. The latched discrete can indicate any power fault condition since the last off-board reset 
condition. 

In operation, the processor board inputs raw antenna data from the base station modem 
20 card 1 12 (or other available location of that data), detects sources of interference within that 
data, and produces a new stream of data which has reduced interference subsequently transmit- 
ting that refined data back to the modem card (or other location) for further processing within 
the base station. 

30 As can be appreciated by one skilled in the art, such interference reduction is computa- 

tionally complex; hence, the hardware must support throughputs sufficient for multiple user 
processing. In a preferred embodiment, characteristics of processing are a latency of less than 
300 microseconds handing data in the 110 Mbytes/Sec range, however, in other embodiments 
the latency and data load can vary. 

35 

In the illustrated embodiment, data from the modem board is supplied via the PCI bus 
211b through the PCI bridge 222. From there, the data traverses the crossbar 206 and is loaded 
into the host controller memory 205. Output data flows in the opposite direction. Additionally, 
certain data flows between the host controller 203 and the compute elements 220. 

40 
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Hybrid Operating System 

The compute elements 220 operate, in some embodiments, under the MC/OS operating 
system available commercially from the assignee herein, although different configurations can 
5 run under different operating systems suited for such. Here, one aspect is to reduce the use of 
non-POSIX system calls which can increase portability of the multiple user detection software 
among different hardware environments and operating system environments. The host proces- 
sor is operated by the Vx Works operating system, as is required by MC/OS and suitable for a 
Motorola 8240 PowerPC. 

10 

Figure 8 shows a block diagram of various components within the hardware/software 
environment. An MC/OS subsystem 802 is used as an operating system for the compute ele- 
ments 220. Further, a MC/OS DX 804 provides APIs acceptable overhead and latency access 
to the DMA engines which in turn provide suitable bandwidth transfers of data. DX 804 can be 
15 used to move data between the compute elements 220 during parallel processing, and also to 
move data between the compute elements 220, the host controller 203, and the modem card 
112. As described above, each compute element 220 continues an application 806, and a 
watchdog 808. Further, the HA registers provide the bootstrap 810 necessary for start-up. 

20 The host controller 203 runs under the Vx Works operating system 812. The host pro- 

cessor 202 contains a watchdog 814, application data 816, and a bootstrap 818. Further, the 
host processor 202 can perform TCP/IP stack processing 820 for communication through the 
Ethernet interface 224. 

30 Input/output between the processor card 118 and the modem card 112 takes place by 

moving data between the Race++ Fabric and the PCI bus 211b via the PCI bridge 222. The 
application 806 will use DX to initialize the PXB++ bridge, and to cause input/output data to 
move as if it were regular DX IPC traffic. For example, there are several components which can 
initiate data transfers and choose PCI addresses to be involved with the transfers. 

35 

One approach to increasing available on the processor card 11 8 is to balance host-pro- 
cessing time against application execution. For example, when the system comes up, the appli- 
cation determines which processing resources are available, and the application determines a 
load mapping on the available resources and record certain parameters in NVRAM. Although 
40 briefs interruptions in service can occur, the application does not need to know how to continue 
execution across faults. For instance, the application can make an assumption that the hardware 
configuration will not change without the system first rebooting. If the application is in a state 
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which needs to be preserved across reboots, the application checkpoints the data on a regular 
basis. The system software provides an API to a portion of the NVRAM for this purpose 

The host controller 203 is attached to an amount of linear flash memory 216 as dis- 
5 cussed above. This flash memory 216 serves several purposes. The first purpose the flash 
memory serves is as a source of instructions to execute when the host controller comes out of 
reset. Linear flash can be addressed much like normal RAM. Flash memories can be organized 
to look like disk controllers; however in that configuration they generally require a disk driver 
to provide access to the flash memory. Although such an organization has several benefits such 
10 as automatic reallocation of bad flash cells, and write wear leveling, it is not appropriate for 
initial bootstrap. The flash memory 216 also serves as a file system for the host and as a place 
to store permanent board information (e.g., such as a serial number). 

When the host controller 203 first comes out of reset, memory is not turned on. Since 
15 high-level languages such as C assume some memory is present (e.g., for a stack) the initial 
bootstrap code must be coded in assembler. This assembler bootstrap contains a few hundred 
lines of code, sufficient to configure the memory controller, initialize memory, and initialize the 
configuration of the host processor internal registers. 

20 After the assembler bootstrap has finished execution, control is passed to the processor 

HA code (which is also contained in boot flash memory). The purpose of the HA code is to 
attempt to configure the fabric, and load the compute element CPUs with HA code. Once this 
is complete, all the processors participate in the HA algorithm. The output of the algorithm is 
a configuration table which details which hardware is operational and which hardware is not. 

30 This is an input to the next stage of bootstrap, the multi-computer configuration. 

MC/OS expects the host controller system to configure the multi-computer (e.g., com- 
pute elements 220). A configmc program reads a textual description of the computer system 
configuration, and produces a series of binary data structures that describe the system configu- 
35 ration. These data structures are used in MC/OS to describe the routing and configuration of the 
multi-computer. 

The processor board 118 will use almost exactly the same sequence to configure the 
multi-computer. The major difference is that MC/OS expects configurations to be static, 
40 whereas the processor board configuration changes dynamically as faulty hardware cause vari- 
ous resources to be unavailable for use. 
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One embodiment of the invention uses binary data structures produced by configmc to 
modify flags that indicate whether a piece of hardware is usable. A modification to MC/OS 
prevents it from using hardware marked as broken. Another embodiment utilizes the output of 
the HA algorithm to produce a new configuration file input to configmc, the configmc execu- 
5 tion is repeated with the new file, and MC/OS is configured and loaded with no knowledge of 
the broken hardware whatsoever. This embodiment can calculate an optimal routing table in the 
face of failed hardware, increasing the performance of the remaining operational components. 

After the host controller has configured the compute elements 220, the runrnc program 
1 0 loads the functional compute elements with a copy of MC/OS. Because access to the proces- 
sor board 118 from a TCP/IP network is required, the host computer system acts as a connec- 
tion to the TCP/IP network. The Vx Works operating system contains a fully functional TCP/IP 
stack. When compute elements access network resources, the host computer acts as proxy, 
exchanging information with the compute element utilizing DX transfers, and then making the 
1 5 appropriate TCP/IP calls on behalf of the compute element. 

The host controller 203 needs a file system to store configuration files, executable pro- 
grams, and MC/OS images. For this purpose, flash memory is utilized. Rather than have a 
separate flash memory from the host controller boot flash, the same flash is utilized for both 
20 bootstrap purposes and for holding file system data. The flash file system provides DOS file 
system semantics as well as write wear leveling. 

There are in particular, two portions of code which can be remotely updated; the boot- 
strap code which is executed by the host controller 203 when it comes out of reset, and the rest 
30 of the code which resides on the flash file system as files. 

When code is initially downloaded to the processor board 1 18, it is written as a group 
of files within a directory in the flash file system. A single top-level index tracks which direc- 
tory tree is used for booting the system. This index continues to point at the existing directory 
35 tree until a download of new software is successfully completed. When a download has been 
completed and verified, the top-level index is updated to point to the new directory tree, the 
boot flash is rewritten, and the system can be rebooted. 

Fault detection and reporting 820, 822 is performed by having each CPU in the system 
40 gather as much information about what it observed during a fault, and then comparing the 
information in order to detect which components could be the common cause of the symptoms. 
In some cases, it may take multiple faults before the algorithm can detect which component is 
at fault. 
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communicate via DMA with processor ^ throughput can be increased, 
transfers can be faster than bus transfers, and hence, thro g 
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or dual-ported RAM is smaller that the data transfer srze. 
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A FPGA 904 is coupled with both the buffer 902 and an crossbar 906 (which can be the 
same crossbar coupled with the compute elements 220 and host controller 203). The FPGA 
904 moves data from the buffer 902 to the crossbar 906, which subsequently communicates the 
data to further devices, e.g., a RACEway™ or the host controller 203 or compute elements 220. 
5 The FPGA 904 also perform data transfers directly from the DSP 900 to the crossbar 906. This 
method is utilized in some embodiments where data transfer sizes can be accommodated with- 
out buffering, for instance, although either the buffer or direct transfers can be used. 

The DSP 900 contains at least one external memory interface (EMIF) 908 device, 
10 which is connected to the buffer 902 or dual-ported RAM. RACEway™ transfers actually 
access the RAM, and then additional processing takes place within the DSP to move the data 
to the correct location in SDRAM within the DSP. In embodiments where the RAM is smaller 
that the data transfer size, then there is a massaging protocol between two endpoint DSPs 
exchanging messages, since the message will be fragmented to be contained within the buffer 
15 or RAM. 

As more RACEway™ endpoints are added (for instance, to increase speed or through- 
put), the size of the dual-port RAM can be increased to a size of 2*F*N*P buffers of size F, 
where F is the fragment size, N is the number of RACEway™ endpoints in communication 

20 with the DSP, and P is the number of parallel transfers which can be active on an endpoint. The 
constant 2 represents double buffering so one buffer can be transferred to the RACEway™ 
simultaneously with a buffer being transferred to the DSR One skilled in the art can appreciate 
that the constant can be four times rather than two times to emulate a full-duplex connection. 
With a 4 mode system, this could be, for example, 4*8K*4*4 or 512 Kbytes, plus a overhead 

30 factor for configuration and data tracking. 

The FPGA 904 can program the DMA controller 910 within the DSP 900 to move data 
between the buffer 902 and the DSP/SDRAM 912 directly from a DSP host port 914. The host 
port 914 is a peripheral like the EMIF 908, but can master transfers into the DSP data-paths, 
35 e.g., it can read and write any location within the DSP. Hence, the host port 914 can access the 
DMA controller, 910 and can be used to initiate transfers via the DMA engine. One skilled in 
the art can appreciate that using this architecture, RACEway™ transfers can be initiated with- 
out the cooperation of the DSP, the thus, the DSP is free to continue processing while transfers 
take place and further, there is no need for protocol messaging within the buffer. 

40 

The FPGA 904 can also perform fragmentation of data. In embodiments where the 
buffer device is a dual-port RAM, the FPGA 904 an program the DMA controller within the 
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DSP to move fragments into or out-of the DSP. This method can be used to match throughput 
of the external transfer bus, e.g., the RACEway™. 

An example of the methods and systems described for a DSP, is as follows. In an 
5 embodiment where the RACEway™ reads date out of the DSP memory 912, this example 
assumes that another DSP is reading the SDRAM of the local DSP. The FPGA 904 detects a 
RACEway™ data packet arriving, and decodes the packet to determine that is contains instruc- 
tions for a data-read at, for example, memory location 0x10000. The FPGA 904 writes over 
the host port interface 914 to program the DMA controller 910 to transfer data starting at 
10 memory location 0x10000, which refers to a location in the primary EMIF 908 corresponding 
to a location in the SDRAM 912, and to move that data to a location in the secondary EMIF 
(e.g., the buffer device) 902. As data arrives in the buffer 902, the FPGA 904 reads the data out 
of the buffer, and moves it onto the RACEway™ bus. When a predetermined block of data is 
moved, the DMA controller 910 finishes the transfer, and the FPGA 904 finishes moving the 
1 5 data from the buffer 902 to the RACEway™. 

Another example assumes that another DSP is requesting a write instruction to the local 
DSP. Here, the FPGA 904 detects a data packet arriving, and determines that is it a write to 
location 0x20000, for instance. The FPGA 904 fills some amount of the buffer 902 with the 
20 data from the RACEway™ bus, and then writes over the host port 914 interface to program the 
DMA controller 910. The DMA controller 910 then transfers data from the buffer device 902 
and writes that data to the primary EMIF 908 at address 0x20000. At the conclusion of the 
transfer, an interrupt can be sent to the DSP 900 to indicate that a data packet has arrived, or a 
polling of a location in the SDRAM 912 can accomplish the same requirement. 

30 

These two examples are non-limiting example, and other embodiments can utilize dif- 
ferent methods and devices for the transfer of data between devices. For example, if the DSP 
900 utilizes RapidlO interfaces, the buffer 902 and FPGA 904 can be modified to accommodate 
this protocol. Also, the crossbar 906 illustrated may be in common with a separate bus struc- 
35 ture, or be in common with the processor board 118 described above. Even further, in some 
embodiments, the FPGA 904 can be directly coupled with the board processor, or be config- 
ured as a compute node 220. 

Therefore, as can be understood by one skilled in the art, the methods and systems 
40 herein are suited for multiple user detection within base stations, and can be used to accom- 
modate both short-code and long-code receivers. 
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Short-Code Processing 
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this period. Performance is improved under high input BER if the sign detector in is replaced 
by the hyperbolic tangent detector. This detector has a single slope parameter which is variable 
from iteration to iteration. Similarly, performance is improved if only a fraction of the total 
estimated interference is cancelled (e.g., partial interference cancellation), owing to channel 
5 and symbol estimation errors. 

Multiple Processors Generating Complementary R-Matrices 

The three R-matrices (R[-l], R[0] and R[l]) are each KvxKv in size. The total number 
1 0 of operation then is 6K* per iteration. The computational complexity of the multistage MDFIC 
algorithm depends on the total number of virtual users, which depends on the mix of users at 
the various spreading factors. For Kv = 200 users (e.g. 100 low-rate users) this amounts to 
240,000 operations. In the current implementation two iterations are used, requiring a total of 
480,000 operations. For real-time operation these operations must be performed in 1/15 ms. 
15 The total processing requirement is then 7.2 GOPS. Computational complexity is markedly 
reduced if a threshold parameter is set such that IC is performed only for values l^/Ml below 
the threshold. The idea is that if \y\ m \ is large there is little doubt as to the sign of b t [m] , and 
IC need not be performed. The value of the threshold parameter is variable from stage to 
stage. 

20 

Although three R matrices are output from the R matrix calculation function, only half 
of the elements are explicitly calculated. This is because of symmetry that exists between R 
matrices: 

30 R i,k=^ R kA~ m \ 

Therefore, only two matrices need to be calculated. The first one is a combination of 
R(l) and R(-l). The second is the R(0) matrix. In this case, the essential R(0) matrix elements 
have a triangular structure to them. The number of computations performed to generate the raw 
35 data for the R(l)/R(-1 ) and R(0) matrices are combined and optimized as a single number. This 
is due to the reuse of the X-matrix outer product values across the two R-matrices. Since the 
bulk of the computations involve combining the X-matrix and correlation values, they domi- 
nate the processor utilization. These computations are used as a cost metric in determining the 
optimum loading of each processor. 

40 
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Processor Loading Optimization 

The optimization problem is formulated as an equal area problem, where the solution 
results in each partition area to be equal. Since the major dimensions of the R-matrices are in 
5 terms of the number of active virtual users, the solution space for this problem is in terms of the 
number of virtual users per processor. By normalizing the solution space by the number of 
virtual users, the solution is applicable for an arbitrary number of virtual users. 

Figure 10 shows a model of the normalized optimization scenario. The computations 
1 0 for the R(l)/R(-1) matrix are represented by the square HJKM, while the computations for the 
R(0)matrix are represented by the triangle ABC. From geometry, the area of a rectangle of 
length b and height h is: 

A r =bh 

15 

For a triangle with a base width b and height h , the area is calculated by: 

A t =-bh 
' 2 

20 When combined with a common height a, the formula for the area becomes: 



Af - A ri + A ti 



1 2 

= aa+-a. 
30 ' e 2 ' 



The formula for A gives the area for the total region below the partition line. For exam- 
ple, the formula for A2 gives the area within the rectangle HQRM plus the region within tri- 
angle AFG. For the cost function, the difference in successive areas is used. That is: 

1 2 1 2 



35 



40 



For an optimum solution, the B must be equal for / =1 ,2,...JI , where TV is the number of 
processors performing the calculations. Because the total normalized load is equal to AN, the 
loading per processor load is equal to AN IN , 
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fori = l,2...,N. 



5 



By combining the two equations for B, the solution for a. is found by finding the roots 
of the equation: 



10 



— a* +a. -— af , = 0 

2 ' 1 2 " l IN < 



The solution for a is: 



15 



a. =-l± 11+ al Vl +2a._ x + — , fori-1,2, ...,N 



Since the solution space must fall in the range [0,1 ], negative roots are not valid solu- 
tions to the problem. On the surface, it appears that the a must be solved by first solving for case 
20 where =1 . However, by expanding the recursions of the a and using the fact that aO equals zero, 
a solution that does not require previous a ,=0,1,. -1 exists. The solution is: 



As shown in the following table, the normalized partition values for two, three, and four 
processors. To calculate the actual partitioning values, the number of active virtual users is 
multiplied by the corresponding table entries. Since a fraction of a user cannot be allocated, a 
35 ceiling operation is performed that biases the number of virtual users per processor towards the 
processors whose loading function is less sensitive to perturbations in the number of users. 
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Location 


Two Processors 


Three Processors 


Four Processors 


a l 

1 


-1 + ^(0.5811) 


r~ 

-1 + V2(0.4142) 


-1 + ^(0.3229) 


a 2 




-1 + V3 (0.7321) 


-1 + ^(0.5811) 


a 3 






-1 + ^(0.8028) 



10 



15 



One skilled in the art can appreciate that the load balancing for the R-matrix results in 
a non-uniform partitioning of the rows of the final matrices over a number of processors. The 
partition sizes increase as the partition starting user index increases. When the system is run- 
ning at full capacity (e.g., all co-processors are functional, and the maximum number of users 
is processed while still within the bounds of real-time operation), and a co-processor fails, the 
impact can be significant. 



20 



This impact can be minimized by allocating the first user partition to the disabled node. 
Also the values that would have been calculated by that node are set to zero. This reduces the 
effects of the failed node. By changing which user data is set to zero (e.g., which users are 
assigned to the failed node) the overall errors due to the lack of non-zero output data for that 
node are averaged over all of the users, providing a "soft" degradation. 
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R,C Values Contiguous in MP1C Processor Memory 



35 



Further, via connection with the crossbar multi-port connector, the multi-processor ele- 
ments calculating the R-matrix (which depends on the C-matrix, which in turn depends on the 
gamma-matrix) can place the results in a processor element performing the MPIC functions. 
For one optimal solution, the values can be placed in contiguous locations accessable (or local 
with) the MPIC processor. This method allows adjacent memory addresses for the R and C 
values, and increases throughput via simply incrementing memory pointers rather that using a 
random access approach. 



^® As discussed above, the values of the r-matrix elements which are non-zero need to be 

determined for efficient storage of the r-matrix. For high data rate users, certain elements 
c/nj are zero, even within the interval n = 0:N-1, N = 256. These zero values reduce the inter- 
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val over which TJm] is non-zero. In order to determine the interval for non-zero values con- 
sider the following relations: 
1 N ~ l 

5 

The index j t for the /th virtual user is defined such that c/nj is non-zero only over the 
interval n = j l N l :j J N t + N l -1. Correspondingly, the vector cjn] is non-zero only over the 
interval n = j k N k : j k N k + N k - 1 . Given these definitions, r /A /m/ can be rewritten as 

1 " H 

The minimum value of m for which TJm] is non-zero is 

15 

and the maximum value of m for which Tjjm] is non-zero is 
20 The total number of non-zero elements is then 

m total = m max2 ~ m min2 +1 

= N,+N k -\ 

30 

The table below provides a sample of the the number of bytes per l,k virtual-user pair 
based on 2 bytes per element - one byte for the real part and one byte for the imaginary part. 
In other embodiments, these values vary. 



35 





N =256 


128 


64 


32 


16 


8 


4 


N. = 256 


1022 


766 


638 


574 


542 


526 


518 


•128 


766 


510 


382 


318 


786 


770 


767 


64 


638 


382 


254 


190 


158 


147 


134 


32 


574 


318 


190 


176 


94 


78 


70 


16 


547 


786 


158 


94 


67 


46 


38 


8 


526 


270 


142 


78 


46 


30 


22 


4 


518 


262 


134 


70 


38 


22 


14 



40 

The memory requirements for storing the T-matrix for a given number of users at each 
spreading factor can be determined as described below. For example, for K virtual users at 
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spreading factor N q = 2*~ q , q = 0:6, where AT is the q\h element of the vector K (some ele- 
ments of K may be zero), the storage requirement can be computed as follows. Let the table 
above be stored in matrix M with elements M . For example, = 1022, and M Q/ = 766. The 
total memory required by the T matrix in bytes is then given by the following relation 

5 

q=0 [ 9-9+1 J 

JO «-° J (27) 

Then, continuing the example, for 200 virtual users at spreading factor N g = 256, A" ? = 
2008 ?0 , which in turn results in = VJCflL, + \)M m = 1 00(20 1)( 1022) = 20.5 MB. For 10 
384 Kbps users , K q = Kjb^ + K 6 h q6 with K Q = 10 and K 6 = 640, which results in a storage 
j 5 requirement that is given by the following relations: 



= + W„ + W<* + '>*W + " 5(1 1)(1022) 

+ 10(640)(518) + 320(641)(14) = 6.2 MB. 

20 

The T-matrix data can be addressed, stored, and accessed as described below. In par- 
ticular, for each pair (Ik), k>= I , there are 1 complex rjmj values for each value of w, where 
m ranges from m mjn2 to m maxV and the total number of non-zero elements is m totgl = m max2 - m min2 
+ 7. Hence, for each pair (l,k), k>= I , there exists 2m (otaJ time-contiguous bytes. 

30 

In one embodiment, an array structure is created to access the data, as shown below: 
struct { 

int m_min2; 
int m_max2; 
int mtotal; 
35 char * Glk; 

} G_info[N_VU_MAX][ N_VU_MAX]; 

The C-matrix data can then be retrieved by utilizing the following exemplary algo- 
rithm: 

40 

m mm2 = G_mfo[l][k].m_min2 
m m 3X 2 ~ G_info[l][k].m max2 
N g = L g /N c 

Nl =m'*N-L/(2N) 
form' = 0:l * 
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forq = 0:L-l 

forq ,: =0:L-l 



10 



m fljifl = max[ m,*, ,"U;J 
suml = 0.0; 

ptrl^&GjnfomM.GIktmJ 

m span 



end 



end 

CK][l][k][q][q1 = suml 



end 



end 



end 



15 
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20 



, SOC U 2 elements to calculate. Assuming all users are at 
Duetosymmet^therearel^K L) g ^ probability that a 

30 SF 256, each calculation requues 256 cmacs c , Q5 if all users are at 1 20 

-pathch^ 
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One skilled in the art can appreciate that a fast fourier transform (FFT) can be used to 
calculate the correlations for a range of offsets, tau, using: 



10 



= C.[W*1] 



nN c +m'T + Tii-Xi q - 



c',[n] 



1 



A A 



The length of the waveform sk[t] is Lg + 255N C - 1068 for L g = 48 and N c = 4. This is 
represented as N c waveforms of length Lg/N c + 255 = 267. One advantage of this approach is 
15 that elements can be stored for a range of offsets tau so that calculations do not need to be per- 
formed when lags change. For delay spreads of about 4 micro-seconds 32 samples need to be 
stored for each m\ 

The C-matrix elements need be updated when the spreading factor changes. The 
20 spreading factor can change du to AMR codec rate changes, multiplexing of the dedicated 
channels, or multiplexing of data services, to name a few reasons. It is reasonable to assume 
that 5% of the users, hence 10% of the elements, change every 10 ms. 
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Gamma-Matrix Generated in FPGA 

The C-matrix elements can be represented in terms of the underlying code correlations 



using: 



35 
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A A 



nNc + m'T+Tiq-Tiq- 

m 

TsmT + T/,-Ti,' 



c,[n) 
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If the length of g[t] is Lg = 48 and Nc = 4, then the summation over m requires 48/4 =12 
macs for the real part and 12 macs for the imaginary part. The total ops is then 48 ops per ele- 
ment. (Compare with 2048 operations for the direct method.) Hence for the case where there 
are 200 virtual users and 20% of the C-matrix needs updating every 10 ms the required com- 

5 plexity is (960000 el)(48 ops/el)(0.20)/(0.010 sec) = 921.6 MOPS. This is the required com- 
plexity to compute the C-matrix from the Tau-matrix. The cost of computing the Tau-matrix 
must also be considered. The Tau-matrix can be efficiently computed since the fundamental 
operation is a convolution of codes with elements constrained to be +/-1 +/-j. Further, the Tau- 
matrix can be calculated using modulo-2 addition (e.g., XOR) using several method, e.g. reg- 

10 ister shifting, XOR logic gates, and so on. 

The Gamma matrix (r) represents the correlation between the complex user codes. 
The complex code for user 1 is assumed to be infinite in lenght, but with only Nj non-zero 
values. The non-zero values are constrained to be ±l±j. The T -matrix can be represented in 
1 5 terms of the real and imaginary parts of the complex user codes, and is based on the relation- 
ship: 




n 

^[m]^(yf[n]®y Y k [n^m]) mf[n\ m Y k [n-m] 



which can be performed using a dual-set of shift registers and a logical circuit contain- 
ing modulo-2 (e.g., Exclusive-OR "XOR") logic elements. Further, one skilled in the art can 
appreciate that such a logic device can be implemented in a field programmable gate array, 
which can be programmed via the host controller, a compute element, or other device including 
an application specific integrated circuit. Further, the FPGA can be progammed via the RACE- 
way™ bus, for example. 

The above shift registers together with a summation device calculates the functions 
Mi r [m] and N? [m].The remaining calculations to form T^m] and subsequently T lk [m] 
can be performed in software. Note that the four functions [m] corrsponding to X, Y = R, I 
which are components of can be calculated in parallel. For K v = 200 virtual users, and assum- 
ing that 10% of all (I, k) pairs must be calculated in 2 ms, then for real-time operation we must 
calculate 0.10(200) 2 = 4000 elements (all shifts) in 2 ms, or about 2M elements (all shifts) per 
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second. For K v = 128 virtual users the requirement drops to 0.8192M elements (all shifts) per 
second. 

In what has been presented the elements are calculated for all 512 shifts. Not all of these 
5 shifts are needed, so it is possible to reduce the number of calculations per elements. The cost 
is increased design complexity. 

Therefore, a possible loading scenaio for performing short-code multiple user detection 
on the hardware described herein is illustrated in Figure 11. A processor board 118 with four 

10 compute elements 220 can be used as shown. Three of the compute nodes (e.g., 220a - 220c) 
can be used to calculate the C-matrix and R-matrix. One of the compute nodes (e.g., 220d) can 
be used for multiple-stage decision-feedback interference cancellation (MDF1C) techniques. 
The Tau-Matrix and R-Matrix is calculated using FPGA's that can be programmed by the host 
controller 203, or ASICs. Further, multiuser amplitude estimation is performed within the 

15 modem card 112. 

Long-Code Processing 

Therefore it can be appreciated by one skilled in the art that short-code MUD can be 
20 performed using the system architecture described herein. Figure 12 shows a preferred 
embodiment for long-code MUD processing. In this embodiment, each frame of data is pro- 
cessed three times by the MUD processor, although it can be recognized that multiple proces- 
sors can perform the iterative nature of the embodiment. During the first pass, only the control 
channels are respread which the maximum ratio combination (MRC) and MUD processing is 
30 performed on the data channels. During subsequent passes, data channels are processed exclu- 
sively. New y (i.e., soft decisions) and b (i.e., hard decisions) data are derived as shown in the 
diagram. 

Amplitude ratios and amplitudes are determined via the DSP (e.g., element 900, or a 
35 DSP otherwise coupled with the processor board 118 and receiver 110), as well as certain 
waveform statistics. These values (e.g., matrices and vectors) are used by the MUD processor 
in various ways. The MUD processor is decomposed into four stages that closely match the 
structure of the software simulation: Alpha Calculation and Respread 1302, raised-cosine fil- 
tering 1304, de-spreading 1306, and MRC 1308. Each pass through the MUD processor is 
40 equivalent to one processing stage of the software implementation. The design is pipelined and 
"parallelized." In the illustrated embodiment, the clock speed can be 132 MHz resulting in a 
throughput of 2.33 ms/frame, however, the clock rate and throughput varies depending on the 
requirements. The illustrated embodiment allows for three-pass MUD processing with addi- 

58 



Express Mail Label: US EV 093 931 8 06 US 



tional overhead from external processing, resulting in a 4-times real-time processing through- 
put. 

The alpha calculation and respread operations 1302 are carried out by a set of thirty-two 
5 processing elements arranged in parallel. These can be processing elements within an ASIC, 
FPGA, PLD or other such device, for example. Each processing element processes two users 
of four fingers each. Values for b are stored in a double-buffered lookup table. Values of a(hat) 
and ja(hat) are pre-multiplied with beta by an external processor and stored in a quad-buffered 
lookup table. The alpha calculation state generated the following values for each finger, where 
1 0 subscripts indicate antenna identifier: 

a 0 =$ 0 (Ca Q -jC ja 0 ) 
ya 0 =p 0 O*C a 0 +C ja 0 ) 
15 al = $ r (Ca y -jCja l ) 

These values are accumulated during the serial processing cycle into four independent 
8-times oversampling buffers. There are eight memory elements in each buffer and the element 
20 used is determined by the sub-chip delay setting for each finger. 

Once eight fingers have been accumulated into the oversampling buffer, the data is 
passed into set of four independent adder-trees. These adder-trees each termination in a single 
output, completing the respread operation. 

30 

The four raised-cosine filters 1 304 convolve the alpha data with a set of weights deter- 
mined by the following equation: 



35 




The filters can be implemented with 97 taps with odd symmetry. The filters illustrated 
run at 8-times the chip rate, however, other rates are possible. The filters can be implemented 
in a variety of compute elements 220, or other devices such as ASICs, FPGAs for example. 
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The despread function 1306 can be performed by a set of thirty-two processing ele- 
ments arranged in parallel. Each processing element serially processes two users of four fin- 
gers each. 

For each finger, one chip value out of eight, selected based on the sub-chip delay, is 
accepted from the output of the raised-cosine filter. The despread state performs the following 
calculations for each finger (subscripts indicate antenna): 

SF-\ 

> ; o = Zcr 0 +ycyr 0 

o 

SF-l 
0 

SF-l 

yi = y L c ' r i+J c 'J r i 

0 

SF-\ 
0 

The MRC operations are carried out by a set of four processing elements arranged in 
parallel, such as the compute elements 220 for example. Each processor is capable of serially 
processing eight users of four fingers each. Values for y are stored in a double-buffered lookup 
table. Values for b are derived from the MSB of the y data. Note that the b data used in the 
MUD stage is independent of the b data used in the respread stage. Values of a and ja <are 
pre-multiplied with p by an external processor and stored in a quad-buffered lookup table. 
Also, +J a > for each channel is stored in a quad-buffered table. 

The output stage contains a set of sequential destination buffer pointers for each chan- 
nel. The data generated by each channel, on a slot basis, is transferred to the RACEway™ 
destination indicated by these buffers. The first word of each of these transfers will contain a 
counter in the lower sixteen bits indicating how many y values were generated. The upper six- 
teen bits will contain the constant value 0xAA55. This will allow the DSP to avoid interrupts 
by scanning the first word of each buffer. 

In addition, the DSP UPDATE register contains a pointer to single RACEway™ loca- 
tion. Each time a slot or channel data is transmitted, an internal counter is written to this loca- 
tion. The counter is limited to 10 bits and will wrap around with a terminal count value of 
1023. 
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The method of operation for the long-code multiple user detection algorithm (LCMUD) 
is as follows. Spread factor for four-channels requires significant amount of data transfer. In 
order to limit the gate count of the hardware implementation, processing an SF4 channel can 
result in reduced capability. 

5 

A SF4 user can be processed on certain hardware channels. When one of these special 
channels is operating on an SF4 user, the next three channels are disabled and are therefore 
unavailable for processing. This relationship is as shown in the following table: 



SF4Chan 


Disabled Channels 


SF4 Chan 


Disabled Channels 


0 


1,2, 3 


32 


33, 34, 35 


4 


5, 6, 7 


36 


37, 38, 39 


8 


9, 10, 11 


40 


41,42, 43 


12 


12, 14, 15 


44 


45, 46, 47 


16 


17, 18, 19 


48 


49, 50, 51 


20 


21,22, 23 


52 


53, 54, 55 


24 


25, 26, 27 


56 


57, 58, 59 


28 


29, 30, 31 


60 


61,62, 63 



The default y and b data buffers do not contain enough space for SF4 data. When a 
channel is operating on SF4 data, the y and b buffers extend into the space of the next channel 
in sequence. For example, if channel 0 is processing SF data, the channel 0 and channel 1 b 
buffers are merged into a single large buffer of 0x40 32-bit words. The y buffers are merged 
similarly. 

In typical operation, the first pass of the LCMUD algorithm will respread the control 
channels in order to remove control interference. For this pass, the b data for the control chan- 
nels should be loaded into BLUT while the y data for data channels should be loaded into 
YDEC. Each channel should be configured to operate at the spread factor of the data channel 
stored into the YDEC table. 

Control channels are always operated at SF 256, so it is likely that the control data will 
need to be replicated to match the data channel spread factor. For example, each bit (b entry) of 
control data would be replicated 64 times if that control channel were associated with an SF 4 
data channel. 

Each finger in a channel arrives at the receiver with a different delay. During the 
Respread operation, this skew among the fingers is recreated. During the MRC stage of MUD 
processing, it is necessary to remove this skew and realign the fingers of each channel. 
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This is accomplished in the MUD processor by determining the first bit available from 
the most delayed finger and discarding all previous bits from all other fingers. The number of 
bits to discard can be individually programmed for each finger with the Discard field of the 
MUDPARAM registers. 

5 

This operation will typically result in a 'short' first slot of data. This is unavoidable 
when the MUD processor is first initialized and should not create any significant problems. The 
entire first slot of data can be completely discarded if 'short' slots are undesirable. 

1 0 A similar situation will arise each time processing is begun on a frame of data. To avoid 

losing data, it is recommended that a partial slot of data from the previous frame be overlapped 
with the new frame. Trimming any redundant bits created this way can be accomplished with 
the Discard register setting or in the system DSP. In order to limit memory requirements, the 
LCMUD FPGA processes one slot of data at a time. Doubling buffering is used for b and y data 

15 so that processing can continue as data is streamed in. Filling these buffers is complicated by 
the skew that exists among fingers in a channel. 

Figure 13 illustrates the skew relationship among fingers in a channel and among the 
channels themselves. The illustrated embodiment allows for 20us (77.8 chips) of skew among 
20 fingers in a channel and certain skew among channels, however, in other embodiments these 
skew allowances vary. 

There are three related problems that are introduced by skew: Identifying frame & slot 
boundaries, populating b and y tables and changing channel constants. 

30 

Because every finger of every channel can arrive at a different time, there are no univer- 
sal frame and slot boundaries. The DSP must select an arbitrary reference point. The data 
stored in b & y tables is likely to come from two adjacent slots. 

35 Because skew exists among fingers in a channel, it is not enough to populate the b & y 

tables with 2,560 sequential chips of data. There must be some data overlap between buffers to 
allow lagging channels to access "old" data. The amount of overlap can be calculated dynami- 
cally or fixed at some number greater than 78 and divisible by four (e.g. 80 chips). The starting 
point for each register is determined by the Chip Advance field of the MUDPARAM register. 

40 

A related problem is created by the significant skew among channels. As can be seen in 
Figure 13, Channel 0 is receiving Slot 0 while Channel 1 is receiving Slot 2. The DSP must 



62 



Express Mail Label: US EV 093 931 8 06 US 



take this skew into account when generating the b and y tables and temporally align channel 
data. 

Selecting an arbitrary "slot" of data from a channel implies that channel constants tied 
5 to the physical slot boundaries may change while processing the arbitrary slot. The Constant 
Advance field of the MUDPARAM register is used to indicate when these constants should 
change. 

Registers affected this way are quad-buffered. Before data processing begins, at least 
10 two of these buffers should be initialized. During normal operation, one additional buffer is 
initialized for each slot processed. This system guarantees that valid constants data will always 
be available. 

The following two tables shown the long-code MUD FPG A memory map and control/ 
15 status register: 



Start Addr j 


End Addr 


Name 


Description 


0000_0000 


0000_0000 


CSR 


Control & Status Register 


0000 0008 


0000JJ00C 


DSPJJPDATE 


Route & Address for DSP updating 


0001J)000 


0001_FFFF 


MUDPARAM 


MUD Parameters 


0002_0000 


0002_FFFF 


CODE 


Spreading Codes 


0003_0000 


0004_FFFF 


BLUT 


Respread: b Lookup Table 


"0005_0000 


0005_FFFF 


BETA_A 


Respread: Beta * a_hat Lookup Table 


0006J)000 


0007_FFFF 


YDEC 


MUD & MRC: y Lookup Table 


0008.0000 


0008_FFFF 


ASQ 


MUD & MRC: Sum a__hat squared LUT 


OOOA_0000 


000A FFFF 


OUTPUT 


Output Routes & Addresses 



35 



Bit 


31 


30 I 29 


28 


27 1 26 | 25 


24 


23 1 22 


21 


20 1 19 


18 


17 |16 


Name 


Reserved 


RAV 


RO 


Reset 


X 


X X 


X 


X 


X 


X 


X 


X 


X 


x ! 


x |x 


X |x 


X 




Bit 


15 


14 | 13 |l2 


11 | 10 |9 


8 


7 |6 


5 


4 


3 


2 


1 


0 


Name 


Reserved 


YB 


CBUF 


A1 


AO 


R1 


R0 


Lst 


Rst 


RAV 


RO 


RO 


RO 


RO 


RO 


Rw 


Rw 


Rw 


Rw 


Reset 


X 


X |x 


X 


X 


X 


X 


0 


0 


0 


0 


0 


0 


0 


0 


0 



The register YB indicates which of two y and b buffers are in use. If the system is cur- 
rently not processing, YB indicates the buffer that will be used when processing is initiated. 
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CBUF indicates which of four round-robin buffers for MUD constants (a A beta) is cur- 
rently in use. Finger skew will result in some fingers using a buffer one in advance of this indi- 
cator. To guarantee that valid data is always available, two full buffers should be initialized 
before operation begins. 

5 

If the system is currently not processing, CBUF indicates the buffer that will be used 
when processing is restarted. It is technically possible to indicate precisely which buffer is in 
use for each finger in both the Respread and Despread processing stages. However, this would 
require thirty-two 32-bit registers. Implementing these registers would be costly, and the infor- 
1 0 mation is of little value. 

Al and AO indicate which y and b buffers are currently being processed. A 1 and AO will 
never indicate T at the same time. An indication of '0' for both A 1 and AO means that MUD 
processor is idle. 

15 

Rl and R0 are writable fields that indicate to the MUD processor that data is available. 
Rl corresponds to y and b buffer 1 and R0 corresponds to y and b buffer 0. Writing a 1 F into 
the correct register will initiate MUD processing. Note that these buffers follow strict round- 
robin ordering. The YB register indicates which buffer should be activated next. 

20 

These registers will be automatically reset to '0' by the MUD hardware once processing 
is completed. It is not possible for the external processor to force a '0' into these registers. 

A ' F in this bit indicates that this is the last slot of data in a frame. Once all available 
30 data for the slot has been processed, the output buffers will be flushed. 

A ' F in this bit will place the MUD processor into a reset state. The external processor 
must manually bring the MUD processor out of reset by writing a '0' into this bit. 

35 DSPJJPDATE is arranged as two 32-bit registers. A RACEway™ route to the MUD 

DSP is stored at address 0x0000 0008. A pointer to a status memory buffer is located at address 
0x0000 J)00C. 

Each time the MUD processor writes a slot of channel data to a completion buffer, an 
40 incrementing count value is written to this address. The counter is fixed at 1 0 bits and will wrap 
around after a terminal count of 1023. 
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A quad-buffered version of the MUD parameter control register exists for each finger to 
be processed. Execution begins with buffer 0 and continues in round-robin fashion. These buf- 
fers are used in synchronization with the MUD constants (Beta * a_hat, etc.) buffers. Each 
finger is provided with an independent register to allow independent switching of constant 
5 values at slot and frame boundaries. The following table shows offsets for each MUD chan- 
nel: 



10 



15 



20 



Offset 


User 




Offset 


User 




Offset 


User 




Offset 


User 


0x0000 


0 


0x0400 


16 


0x0800 


32 


OxOCOO 


48 


0x0040 


1 


0x0440 


17 


0x0840 


33 ; 


0x0C40 


49 


0x0080 


2 


0x0480 


4 O 

1o 




OxOooU 


34 ; 


UXUUoU 




OxOOCO 


3 


0x04C0 


19 




0x08C0 


35 


OxOCCO 


51 


0x0100 


4 


0x0500 I 


20 


0x0900 


36 


OxODOO 


52 


0x0140 


5 


0x0540 


21 


0x0940 


37 


0x0D40 


53 


0x0180 


6 


0x0580 


22 


0x0980 


38 


0x0D80 


54 


0x01 CO 


7 


0xO5CO 


23 


0x09C0 


39 


OxODCO 


55 


0x0200 


8 


0x0600 


24 


OxOAOO 


40 


OxOEOO 


56 


0x0240 


9 


0x0640 


25 


0x0A40 


41 


OxOE40 


57 


0x0280 


10 


0x0680 


26 


0x0A80 


42 


0x0E80 


58 


0x02C0 


11 


0x06C0 


27 


OxOACO 


43 


OxOECO 


59 


0x0300 


12 


0x0700 


28 


OxOBOO 


44 


OxOFOO 


60 


0x0340 


13 


0x0740 


29 


0x0B40 


45 


0x0F40 


61 


0x0380 


14 


0x0780 


30 


0x0B80 


46 


0x0F80 


62 


0x03C0 


15 


0x07C0 


31 


OxOBCO 


47 


OxOFCO 
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30 The following table shows buffer offsets within each channel: 



35 



40 



Offset 


Finger 


Buffer 


0x0000 


0 


0 


0x0004 




1 


0x0008 




2 


OxOOOC 




3 


0x0010 


1 


0 


0x0014 




1 \ 


0x0018 




2 


0x001 C 




3 


0x0020 


2 


0 


0x0024 




1 


0x0028 




2 


0x002C 




3 
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0x0030 


3 


0 


0x0034 




1 


0x0038 




2 


0x003C 




3 



5 

The following table shown details of the control register: 



10 



15 



Bit 


31 


30 


29 


28 


27 


26 


25 


24 


23 


22 


21 


20 


19 | IS 


17 


16 


Name 


Spread Factor 


Subchip 
Delay 


Discard 


RAV 


RW 


RW 


RW 


Reset 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


x 


X 


X 


X 


X 


Bit 


15 


14 


13 


12 


11 


10 


9 


8 


7 


6 


5 


4 


3 


2 


1 


0 


Name 


Chip Advance 


Constant Advance 


RAV 


RW 


RW 


Reset 


X 


x |x 


X 


X 


X 


X 


X 


X 


X 


X |x 


x |x 


x |x 



The spread factor field determines how many chip samples are used to generate a data 
bit. In the illustrated embodiment, all fingers in a channel have the same spread factor setting, 
however, it can be appreciated by one skilled in the art that such constant factor setting can be 
variable in other embodiments. The spread factor is encoded into a 3-bit value as shown in the 
following table: 



35 



SF Factor 


Spread Factor 


000 


256 


001 


128 


010 


64 


011 


32 


100 


16 


101 


8 


110 


4 


111 


RESERVED 



40 The field specifies the sub-chip delay for the finger. It is used to select one of eight 

accumulation buffers prior to summing all Alpha values and passing them into the raised- 
cosine filter. 
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Discard determines how many MUD-processed soft decisions (y values) to discard at 
the start of processing. This is done so that the first y value from each finger corresponds to the 
same bit. After the first slot of data is processed, the Discard field should be set to zero. 

5 The behavior of the discard field is different than that of other register fields. Once a 

non-zero discard setting is detected, any new discard settings from switching to a new table 
entry are ignored until the current discard count reaches zero. After the count reaches zero, a 
new discard setting may be loaded the next time a new table entry is accessed. 

10 All fingers within a channel will arrive at the receiver with different delays. Chip 

Advance is used to recreate this signal skew during the Respread operation. Y and b buffers are 
arranged with older data occupying lower memory addresses. Therefore, the finger with the 
earliest arrival time has the highest value of chip advance. Chip Advanced need not be a mul- 
tiple of Spread Factor. 

15 

Constant advance indicates on which chip this finger should switch to a new set of con- 
stants (e.g. a' ) and a new control register setting. Note that the new values take effect on the 
chip after the value stored here. For example, a value of 0x0 would cause the new constants to 
take effect on chip 1. A value of OxFF would cause the new constants to take effect on chip 0 
20 of the next slot. The b lookup tables are arranged as shown in the following table. B values 
each occupy two bits of memory, although only the LSB is utilized by LCMUD hardware. 



30 



35 



40 



67 



Express Mail Label: US EV 093 931 8 06 US 



5 



10 



Offset 


Buffer 


OfTset 


Buffer 


Offset 


Buffer 


Offset 


Buffer 


0x0000 


UO BO 


0x0400 


U16B0 


0x0800 


U32 BO 


OxOCOO 


U48 BO 


0x0020 


U1B0 


0x0420 


U17B0 


0x0820 


U33 BO 


OxOC20 


U49 BO 


0x0040 


U0B1 


0x0440 


U16B1 


0x0840 


U32B1 


0x0C40 


U48B1 


0x0060 


U1B1 


0x0460 ; 


U17B1 


0x0860 


U33B1 


OxOC60 


U49B1 


0x0080 


U2 BO 


0x0480 


U18B0 


0x0880 


U34 BO 


OxOC80 


U50 BO 


OxOOAO 


U3 BO 


0x04A0 


U19B0 


0x08A0 


U35 BO 


OxOCAO 


U51 BO 


OxOOCO 


U2B1 


OxO4C0 


U18B1 


OxOSCO 


U34B1 


OxOCCO 


U50B1 


OxOOEO 


U3B1 


0x04E0 


U19B1 


Ox08EO 


U35B1 


OxOCEO 


U51 Bl 


0x0100 


U4 BO 


0x0500 


U20 BO 


0x0900 


U36 BO 


OxODOO 


U52 BO 


0x0120 


U5 BO 


0x0520 


U21 BO 


0x0920 


U37 BO 


OxOD20 


U53 BO 


0x0140 


U4B1 


0x0540 


U20B1 


0x0940 


U36B1 


0x0D40 


U52B1 


0x0160 


U5B1 


0x0560 


U21 Bl 


0x0960 


U37B1 


OxOD60 


U53B1 


0x0180 


U6 BO 


0x0580 


U22 BO 


0x0980 


U38B0 


0x0D80 


U54 BO 


0x01 AO 


U7 BO 


0xO5AO 


U23 BO 


0x09A0 


U39 BO 


OxODAO 


U55 BO 


0x0 ICO 


U6B1 


Ox05CO 


U22B1 


0x09C0 


U38B1 


OxODCO 


U54B1 


0x01 EO 


U7B1 


Ox05EO 


U23B1 


Ox09EO 


U39B1 


OxODEO 


U55B1 


0x0200 


U8 BO 


0x0600 


U24 BO 


OxOAOO 


U40 BO 


OxOEOO 


U56B0 


0x0220 


U9 BO 


0x0620 


U25B0 


0x0A20 


U41 BO 


OxOE20 


U57 BO 


0x0240 


U8B1 


0x0640 


U24B1 


0x0A40 


U40B1 


0x0E40 


U56B1 


0x0260 


U9B1 


0x0660 


U25B1 


0x0A60 


U41 Bl 


OxOE60 


U57B1 


0x0280 


U10B0 


0x0680 


U26 BO 


0x0A80 


U42 BO 


OxOE80 


U58B0 


0x02A0 


UMBO 


0x06A0 


U27B0 


OxOAAO 


U43 BO 


OxOEAO 


U59B0 


0x02C0 


U10B1 


0x06C0 


U26B1 


OxOACO 


U42B1 


OxOECO 


U58B1 


Ox02EO 


Ull Bl 


0x06E0 


U27B1 


OxOAEO 


U43B1 


OxOEEO 


U59B1 


0x0300 


U12B0 


0x0700 


U28B0 


OxOBOO 


U44 BO 


OxOFOO 


U60 BO 


0x0320 


U13B0 


0x0720 


U29 BO 


0x0B20 


U45 BO 


OxOF20 


U61 BO 


0x0340 


U12B1 


0x0740 


U28B1 


0x0B40 


U44B1 


0x0F4O 


U60B1 


0x0360 


U13B1 


0x0760 


U29B1 


0x0B60 


U45B1 


OxOF60 


U61 Bl 


0x0380 


U14B0 


0x0780 


U30B0 


OxOBSO 


U46 BO 


OxOF80 


U62 BO 


0x03A0 


U15B0 


0x07A0 


U31 BO 


OxOBAO 


U47 BO 


OxOFAO 


U63 BO 


Ox03CO 


U14B1 


0x07C0 


U30B1 


OxOBCO 


U46B1 


OxOFCO 


U62BI 


Ox03EO 


U15B1 


Ox07EO 


U31 Bl 


OxOBEO 


U47B1 


OxOFEO 


U63B1 



The following table illustrates how the two-bit values are packed into 32-bit words. 
Spread Factor 4 channels require more storage space than is available in a single channel buffer. 
To allow for SF4 processing, the buffers for an even channel and the next highest odd channel 
are joined together. The even channel performs the processing while the odd channel is dis- 
abled. 



40 



Bit 


31 | 30 


29 | 28 


27 | 26 


25 | 24 


23 | 22 


21 | 20 


19 | 18 


17 | 16 


Name 


b(0) 


b(l) 


b(2) 


b(3) 


b(4) 


b(5) 


b(6) 


b(7) 


Bit 


15 14 


13 12 


11 10 


9 8 


7 6 


5 4 


3 2 


1 0 


Name 


b(8) 


b(9) 


b(10) 


b(H) 


b(12) 


b(13) 


b(14) 


b(15) 



68 



Express Mai! Label: US EV 093 931 8 06 US 



The beta*a-hat table contains the amplitude estimates for each finger pre-multiplied by 
the value of Beta. The following table shows the memory mappings for each channel. 



5 



Offset 


User 


Offset 


User 


Offset 


User 


Offset 


User 


0x0000 


0 


0x0800 


16 


0x1000 


32 


0x1800 


48 


0x0080 


1 


0x0880 


17 


0x1080 


33 


0x1880 


49 | 


0x0100 


2 


0x0900 


18 


0x1100 


34 


0x1900 


50 


0x0180 


3 


0x0980 


19 


0x1180 


35 


0x1980 


51 j 


0x0200 


4 


OxOAOO 


20 


0x1200 


36 


OxIAOO 


52 


0x0280 


5 


0x0A80 


21 


0x1280 


37 


0x1A80 


53 


0x0300 


6 


OxOBOO 


22 


0x1300 


38 


0x1 BOO 


54 


0x0380 


7 


0x0B80 


23 


0x1380 


39 


0x1 B80 


55 


0x0400 


8 


OxOCOO 


24 


0x1400 


40 | 


0x1 COO 


56 


0x0480 


9 


0x0C80 


25 


0x1480 


41 I 


0x1 C80 


57 


0x0500 


10 


OxODOO 


26 


0x1500 


42 


0x1 D00 


58 


0x0580 


11 


0x0D80 


27 


0x1580 


43 


0x1 D80 


59 ! 


0x0600 


12 


OxOEOO 


28 


0x1600 


44 


0x1 E00 


60 


0x0680 


13 


0x0E80 


29 


0x1680 


45 


0x1 E80 


61 


0x0700 


14 


OxOFOO 


30 


0x1700 


46 


0x1 F00 


62 


0x0780 


15 


0x0F80 


31 


0x1780 


47 


0x1 F80 
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The following table shows buffers that are distributed for each channel: 



30 



Offset 


User Buffer 


0x00 


0 


0x20 


1 


0x40 


2 


0x80 


3 







The following table shows a memory mapping for individual fingers of each antenna. 



40 



Offset 


Finger 


Antenna 


0x00 


0 


0 


0x04 


1 




0x08 


2 




OxOC 


3 




0x10 


0 


1 


0x14 


1 




0x18 


2 




0x1 C 


3 
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The y (soft decisions) table contains two buffers for each channel. Like the b lookup 
table, an even and odd channel are bonded together to process SF4. Each y data value is stored 
as a byte. The data is written into the buffers as packed 32-bit words. 



5 



Offset 


Buffer 


Offset 


Buffer 


Offset 


Buffer 


Offset 


Buffer 


0x0000 


UO BO 


0x4000 


U16B0 


0x8000 


U32 BO 


OxCOOO 


U48 BO 


0x0200 


Ul BO 


0x4200 


U17B0 


0x8200 


U33 BO 


OxC200 


U49 BO 


0x0400 


U2B1 


0x4400 


U18B1 


0x8400 


U34B1 


0xC400 


U50B1 


0x0600 


U3B1 


0x4600 


U19B1 


0x8600 


U35B1 


OxC600 


U51 Bl 


0x0800 


UO BO 


0x4800 


U16B0 


0x8800 


U32 BO 


OxC800 


U48 BO 


OxOAOO 


Ul BO 


0x4A00 


U17B0 


Ox8AO0 


U33 BO 


OxCAOO 


U49 BO 


OxOCOO 


U2B1 


0x4C00 


U18B1 


Ox8COO 


U34B1 


OxCCOO 


U50B1 


OxOEOO 


U3B1 


0x4E00 


U19B1 


Ox8EOO 


U35B1 


OxCEOO 


U51 Bl 


0x0000 


U4 BO 


0x5000 


U20 BO 


0x9000 


U36 BO 


OxDOOO 


U52 BO 


0x0200 


U5 BO 


0x5200 


U21 BO 


0x9200 


U37 BO 


OxD200 


U53 BO 


0x0400 


U6B1 


0x5400 


U22B1 


0x9400 


U38B1 


0xD400 


U54 Bl 


0x0600 


U7B1 


0x5600 


U23B1 


0x9600 


U39B1 


OxD600 


U55 Bl 


0x0800 


U4 BO 


0x5800 


U20 BO 


0x9800 


U36B0 


OxD800 


U52 BO 


OxOAOO 


U5 BO 


Ox5AOO 


U21 BO 


0x9A00 


U37B0 


OxDAOO 


U53 BO 


OxOCOO 


U6B1 


Ox5COO 


U22B1 


0x9C00 


U38B1 


OxDCOO 


U54 Bl 


OxOEOO 


U7B1 


Ox5E0O 


U23B1 


Ox9EOO 


U39B1 


OxDEOO 


U55B1 1 


0x0000 


U8 BO 


0x6000 


U24 BO 


OxAOOO 


U40 BO 


OxEOOO 


U56 BO 


0x0200 


U9 BO 


0x6200 


U25 BO 


0xA200 


U41 BO 


0xE200 


U57B0 


0x0400 


U10B1 


0x6400 


U26B1 


0xA400 


U42B1 


0xE400 


U58B1 


0x0600 


Ull BI 


0x6600 


U27B1 


OxA600 


U43 Bl 


0xE600 


U59B1 


0x0800 


U8 BO 


0x6800 


U24 BO 


OxA800 


U40 BO 


OxE800 


U56 BO 


OxOAOO 


U9 BO 


0x6A00 


U25 BO 


OxAAOO 


U41 BO 


OxEAOO 


U57 BO 


OxOCOO 


U10B1 


0x6C00 


U26B1 


OxACOO 


U42B1 


OxECOO 


U58B1 


OxOEOO 


Ull Bl 


0x6E00 


U27B1 


OxAEOO 


U43B1 


OxEEOO 


U59B1 


0x0000 


U12B0 


0x7000 


U28 BO 


OxBOOO 


U44 BO 


OxFOOO 


U60 BO 


0x0200 


U13B0 


0x7200 


U29 BO 


OxB200 


U45 BO 


0xF200 


U61 BO 


0x0400 


U14B1 


0x7400 


U30B1 


0xB400 


U46B1 


0xF400 


U62B1 


0x0600 


U15B1 


0x7600 


U31 Bl 


OxB600 


U47B1 


0xF600 


U63B1 


0x0800 


U12B0 


0x7800 


U28 BO 


OxB800 


U44 BO 


OxF800 


U60 BO 


OxOAOO 


U13B0 


0x7A00 


U29 BO 


OxBAOO 


U45 BO 


OxFAOO 


U61 BO 


OxOCOO 


U14B1 


Ox7COO 


U30B1 


OxBCOO 


U46B1 


OxFCOO 


U62B1 


OxOEOO 


U15B1 


Ox7E00 


U31 Bl 


OxBEOO 


U47B1 


OxFEOO 


U63B1 



The sum of the a-hat squares is stored as a 1 6-bit value. The following table contains a 
35 memory address mapping for each channel. 



40 



0x0000 


0 


0x0200 


16 


0x0400 


32 


0x0600 


48 


Offset 


User 


Offset 


User 


Offset 


User 


Offset 


User 


0x0020 


1 


0x0220 


17 


0x0420 


33 


0x0620 


49 | 


0x0040 


2 


0x0240 


18 


0x0440 


34 


0x0640 


50 


0x0060 


3 


0x0260 


19 


0x0460 


35 


0x0660 


51 


0x0080 


4 


0x0280 


20 


0x0480 


36 


0x0680 


52 


OxOOAO 


5 


0x02A0 


21 


0x04AO 


37 


0x06A0 


53 


OxOOCO 


6 


0x02C0 


22 


0x04C0 


38 


0x06C0 


54 
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5 



OxOOEO 


7 


0x02E0 


23 


0x04E0 


39 


0x06E0 


55 


0x0100 


8 


0x0300 


24 


0x0500 


40 


0x0700 


56 


0x0120 


9 


0x0320 


25 


0x0520 


41 


0x0720 


57 


0x0140 


10 


0x0340 


26 


0x0540 




0x0/40 


DO 


0x0160 


11 


0x0360 


27 


0x0560 


43 


0x0760 


59 


0x0180 


12 


0x0380 


28 


0x0580 


44 


0x0780 


60 


OxOIAO 


13 


0x03A0 


29 


0x05A0 


45 


0x07A0 


61 


0x01 CO 


14 


0x03C0 


30 


0x05C0 


46 


0x07C0 


62 


0x01 EO 


15 


0x03E0 


31 


OxOSEO 


47 


0x07E0 


63 



10 

Within each buffer, the value for antenna 0 is stored at address offset 0x0 with the value 
for antenna one stored at address offset 0x04. The following table demonstrates a mapping for 
each finger. 



15 



Offset 


User Buffer 


0x00 


0 


0x08 


1 


0x10 


2 


0x1 C 


3 



Each channel is provided a RACE way™ route on the bus, and a base address for buffer- 
ing output on a slot basis. Registers for controlling buffers are allocated as shown in the fol- 
lowing two tables. External devices are blocked from writing to register addresses marked as 
reserved. 
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Offset 


User 


Offset 


User 


Offset 


User 


Offset 


User 


0x0000 


0 


0x0200 


16 


0x0400 


32 


0x0600 


48 


0x0020 


1 


0x0220 


17 


0x0420 


33 


0x0620 


49 


0x0040 


2 


0x0240 


18 


0x0440 


34 


0x0640 


50 


0x0060 


3 


0x0260 


19 


0x0460 


35 


0x0660 


51 


0x0080 


4 


0x0280 


20 


0x0480 


36 


0x0680 


52 


OxOOAO 


5 


0x02A0 


21 


0x04A0 


37 


0x06A0 


53 


0x0000 


6 


0x02C0 


22 


0x04C0 


38 


0x06C0 


54 


OxOOEO 


7 


0x02E0 


23 


0x04E0 


39 


0x06E0 


55 


0x0100 


8 


0x0300 


24 


0x0500 


40 


0x0700 


56 


0x0120 


9 


0x0320 


25 


0x0520 


41 


0x0720 


57 


0x0140 


10 


0x0340 


26 


0x0540 


42 


0x0740 


58 


0x0160 


11 


0x0360 


27 


0x0560 


43 


0x0760 


59 


0x0180 


12 


0x0380 


28 


0x0580 


44 


0x0780 


60 


OxOIAO 


13 


0x03A0 


29 


0x05A0 


45 


0x07A0 


61 


0x01 CO 


14 


0x03C0 


30 


OxOSCO 


46 


0x07C0 


62 
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15 |0x03E0 1 31 |0x05E0 j 47 1 0x07E0 [63 | 



5 



Offset 


Entry 


0x0000 1 


Route to Channel Destination 


0x0004 


Base Address for Buffers 


0x0008 


Buffers 


OxOOOC 


RESERVED 


0x0010 


RESERVED 


0x0014 


RESERVED 


0x0018 


RESERVED 


0x001 C 


RESERVED 



Slot buffer size is automatically determined by the channel spread factor. Buffers are 
used in round-robin fashion and all buffers for a channel must be arranged contiguously. The 
j 5 buffers control register determines how many buffers are allocated for each channel. A setting 
of 0 indicates one available buffer, a setting of 1 indicates two available buffers, and so on. 
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30 



35 
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Methods for Estimating Symbols Embodied In Short-Code User Wave- 
forms 

As discussed above, systems according to the invention perform multi-user detection by 
5 determining correlations among the user channel-corrupted waveforms and storing these cor- 
relations as elements of the R-matrices. The correlations are updated in real time to track con- 
tinually changing channel characteristics. The changes can stem from changes in user code 
correlations, which depend on the relative lag among various user multi-path components, as 
well as from the much faster variations of the Rayleigh-fading multi-path amplitudes. The 
10 relative lags among multi-path components can change with a time constant, for example, of 
about 400 ms whereas the multi-path amplitudes can vary temporally with a time constant of, 
for example, 1.33 ms. The R-matrices are used to cancel the multiple access interference 
through the Multi-stage Decision-Feedback Interference Cancellation (MDFIC) technique. 

1 5 In the preceding discussion and those that follow, the term physical user refers to a 

CDMA signal source, e.g., a user cellular phone, modem or other CDMA signal source, the 
transmitted waveforms from which are processed by a base station and, more particularly, by 
MUD processing card 118. In the illustrated embodiment, each physical user is considered to 
be composed of a one or more virtual users and, more typically, a plurality of virtual users. 

20 

A virtual user is deemed to "transmit" a single bit per symbol period, where a symbol 
period can be, for example, a time duration of 256 chips (1/15 ms). Thus, the number of virtual 
users, for a given physical user, is equal to the number of bits transmitted in a symbol period. 
In the illustrated embodiment, each physical user is associated with at least two virtual users, 
30 one of which corresponds to a Dedicated Physical Control Channel (DPCCH) and the other of 
which corresponds to a Dedicated Physical Data Channel (DPDCH). Other embodiments may 
provide for a single virtual user per physical user, as well, of course, to three or more virtual 
users per physical user. 

35 In the illustrated embodiment, when a Spreading Factor (SF) associated with a physical 

user is less than 256, the J = 256/SF data bits and one control bit are transmitted per symbol 
period. Hence, for the r th physical user with data-channel spreading factor SF, there are a total 
of 1 + 256/SF virtual users. The total number of virtual users can then be denoted by: 



40 




r-\ 

The waveform transmitted by the rth physical user can be written as: 



(1) 
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x r [<]=lhI,s k [t-mT]b k [m} 

s k [t]^h[t-pN e ]p k [p] (2) 

5 

where / is the integer time sample index, T = NN c represents the data bit duration, N = 
256 represents short-code length, N c is the number of samples per chip, and where fi k =J$ e if the 
*th virtual user is a control channel and fi k = fi d if the k\h virtual user is a data channel. The 
multipliers fi c and fi d are utilized to select the relative amplitudes of the control and data chan- 
IQ nels. In the illustrated embodiment, at least one of the above constants equals 1 for any given 
symbol period, m. 

The waveform sk[t], which is herein referred to as the transmitted signature waveform 
for the kth virtual user, is generated by the illustrated system by passing the spread code 
j 5 sequence ck[n] through a root-raised-cosine pulse shaping filter h[t]. When the kth virtual user 
corresponds to a data user with a spreading factor that is less than 256, the code ck[n] retains a 
length of 256, but only Nk of these 256 elements are non-zero, where Nk is the spreading factor 
for the kth virtual user. The non-zero values are extracted from the code C ch 256 M [n]-s sh [n]. 

20 The baseband received signal can be written as: 

r[t] = ^J,s k [t-mT]b k [m]^t] 

*=1 m 
L 

30 

where wftj is receiver noise, s k [t] is the channel-corrupted signature waveform for 
virtual user k t L is the number of multipath components, and a kq . are the complex multipath 
amplitudes. The amplitude ratios fi k are incorporated into the amplitudes a k(f ,. If A and / are two 
virtual users that correspond to the same physical user.then, aside from scaling factors fi k and 
35 fir > a kq- and a i q - are e 9 ual - Th* s is due t0 * e ^ act ^ { ^ e signal waveforms of all virtual users 
corresponding to the same physical user pass through the same channel. Further, the waveform 
sj'ij represents the received signature waveform for the fc* h virtual user, and it differs from the 
transmitted signature waveform given in Equation (2) in that the root-raised-cosine pulse h[t] 
is replaced with the raised-cosine pulse gftj. 

40 

The received signal that has been match-filtered to the chip pulse is also match-filtered 
in the illustrated embodiment to the user code sequence in order to obtain detection statistic, 
herein referred to asy k , for the k xh virtual user. Because there are K v codes, there are K v such 
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detection statistics. For each virtual user, the detection statistics can be collected into a column 
vector y[m] whose entry corresponds to the m ih symbol period. More particularly, the 
matched filter output y/mj for the I th virtual user can be written as: 



*[m] S Rej£fl-~ 



(4) 



where a is an estimate of ai , f is an estimate of \ , N, is the (non-zero) length of 

lq 

10 code c f [nj, and r\/m] represents the match-filtered receiver noise. Substituting the expression 
for r[t] from Equation (3) in Equation (4) results in the following equation: 



15 



«• *=l [»=l „ J 



20 



30 



= X X r « t m i 6 * [ m - w i + 



m ' k=\ 



r, k [m ] h Re j £ a* ~ £ 5, [n # f + T„ + m T] • c,* [n] 

= £ X ^ J ^ c„ . • -i- X ** [«N e + « T + \ - V .] • c, [n] 

The terms for m * = 0 result from asynchronous users. 



(5) 



Calculation of the R-matrix 



35 Determination of the R-matrix elements defined by Equation (5) above can be divided 

into two or more separate calculations, each having an associated time constant or period of 
execution corresponding to a time constant or period during which a corresponding character- 
istic of the user waveforms are expected to change in real time. In the illustrated embodiment, 
three sets of calculations are employed as reflected in the following equations: 

40 
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L L 



10 



15 



20 



30 



35 



40 



»=■ 9=1 
L L 



«;«*,-~II«i(«-pW+*T+t % -v]cj/»].c;[ji]} 

^/V/ op J 



£"l It p 

= ^TrZ 8l mN < + ">r + X lq - \.]^c k [n - m]c,[n] 

= ^rI«i" iy . + " T+l i»-vi r iW 



r « M = S c * l" ~ m l * c # W (6) 



where the hats ( A ), indicating parameter estimates, have been omitted. 

With reference to Equation (6), the r-matrix, whose elements vary with the slowest 
time constant, represents the user code correlations for all values of offset m. For the case of 
100 voice users, the total memory requirement for storing the r-matrix elements is 21 Mbytes 
based on two bytes (e.g., the real and imaginary parts) per element. In the illustrated embodi- 
ment, the r-matrix matrix is updated only when new codes associated with new users are 
added to the system. Hence, the r-matrix is effectively a quasi-static matrix, and thus, its com- 
putational requirements are minimal. 

The selection of the most efficient method for calculating the r-matrix elements 
depends on the non-zero length of the codes. For example, the non-zero length of the codes in 
case of high data-rate users can be only 4 chips long. In such a case, a direct convolution, e.g., 
convolution in the time domain, can be the most efficient method of calculating the elements 
of the r-matrix. For low data-rate users, it may be more efficient to calculate the elements of 
the r-matrix by utilizing Fast Fourier Transforms (FFTs) to perform convolutions in the fre- 
quency domain. 

In one method according to the teachings of the invention, the C-matrix elements are 
calculated by utilizing the r-matrix elements. The C-matrix elements need to be calculated 
upon occurrence of a change in a user's delay lag (e.g., time-lag). For example, consider a case 
in which each multi-path component changes on average every 400 ms, and the length of the 
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gfj function is 48 samples. In such a case, assuming an over-sampling by four, then forty-eight 
operations per element need to be performed (for example, 12 multiple accumulations, real x 
complex, for each element). Further, if 100 low-rate users (i.e., 200 virtual users) are utilizing 
the system, and assuming a single multipath lag of four changes for one user, a total of 

5 (1 .5)(2)K JLN v elements need to be calculated. The factor of 1 .5 arises from the three C-matri- 
ces (e.g., m' = -1,0,1) which is reduced by a factor two as a result of a conjugate symmetry 
condition. Moreover, the factor two arises based on the fact that both rows and columns need 
to be updated. The factor N v represents the number of virtual users per physical user, which for 
the lowest rate users is N = 2 as stated above. In total, this amounts to approximately 230,400 

1 0 operations per multipath component per physical user. Accordingly, it gives rise to 230 MOPS 
based on 100 physical users with four multipath components per user, each changing once per 
400 msec. Of course, in other embodiments these values can differ. 

The C-matrices are then utilized to calculate the R-matrices. More particularly, the ele- 
15 ments of the R-matrix can be obtained as follows by utilizing Equation (6) above: 

*[«1 = li> Ja; V .C^|ml| = Rc{fl / i, .C l >l.a 4 } (7) 

<H *=i I J 

where a k are L x 1 vectors, and C lk [m 'J are L x L matrices. The rate at which the above 
20 calculations need to be performed depends on the velocity of the users. For example, in one 
embodiment, the update rate is selected to be 1 .33 msec. An update rate that is too slow such 
that the estimated values of the R-matrix deviate significantly from the actual R-matrix values 
results in a degradation of the MUD efficiency. For example, Figure 14 presents a graph that 
depicts the change in the MUD efficiency versus user velocity for an update rate of 1 .33 msec, 
30 which corresponds to two WCDMA time slots. This graph indicates that the MUD efficiency 
is high for users having velocities that are less than about 1 00 km/h. The graph further shows 
that the interference corresponding to fast users is not canceled as effectively as the interfer- 
ence corresponding to slow users. Thus, for a system that is utilized by a mix of fast and slow 
users, the total MUD efficiency is an average of the MUD efficiency for the range of user 
35 velocities. Utilizing the above Equation (7), the R-matrix elements can be calculated in terms 
of an X matrix that represents amplitude-amplitude multiplies as shown below: 



40 
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r lk [m ] = Re{/r [a," • Q [m ] • a k ]} = Re [tr [Q [m ] • a t • a," ]} h Re {tr [C a [m\ X lk ]} 
= /r[C>] -**]-/r[C,>l ^] 

5 X lk =a k a" =X? k + jX' lk 

C n [m'\ = CZ[m'\+jC! k [m'\ < 8 ) 



10 



15 



20 



30 



35 



40 



The use of the X-matrix as illustrated above advantageously allows reusing the X- 
matrix multiplies for all virtual users associated with a physical user and for all m ' (i.e., m = 0, 
1 ). The remaining calculations can be expressed as a single real dot product of length 2L2 = 32. 
The calculations can be performed, for example, in 16-bit fixed point math. Then, the total 
operations can amount to 1.5(4)(KL)2 = 3.84 MOPS resulting in a processing requirement of 
2.90 GOPS. The X-matrix multiplies, when amortized, amount to an additional 0.7 GOPS. 
Thus, the total processing requirement can be 3.60 GOPS. 

The matched-filter outputs can be obtained from the above Equation (5) as follows: 

K K K 

y,M = r t [0ft[m]+£rJ-^ (9) 

wherein the first term represents a signal of interest, and the remaining terms represent 
Multiple Access Interference (MAI) and noise. The illustrated embodiment uses a Multistate 
Decision Feedback interference Cancellation (MDFIC) algorithm can be utilized to solve for 
the symbol estimates in accord with the following relationship: 

iK K K 

y,['n]-X^J- 1 ]^[™ +1 ]-2hJ 0 ]-^ 0 l 5 *]**[™]-£^[ , ]^['«- 1 ][ < 10 ) 

with initial estimates given by hard decisions on the matched-filter detection statistics, 
h i [m] = sign{y l [m]}^ 

A further appreciation of these and alternate MDFIC techniques may be attained by 
reference to An MDFIC technique which is described in an article by T.R. Giallorenzi and S. G. 
Wilson, titled, "Decision feedback multi-user receivers for asynchronous CDMA systems", 
published in IEEE Global Telecommunications Conference, pages 1 677- 1 682 (June 1 993), and 
herein incorporated by reference. Related techniques, known as , is closely related to Succes- 
sive Interference Cancellation (SIC) and Parallel Interference Cancellation (PIC), can be used 
in addition or instead. 
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In the illustrated embodiment, the new estimates b f [m] are immediately introduced back 
into the interference cancellation as they are calculated. Hence at any given cancellation step, 
the best available symbol estimates are used. In one embodiment, the above iteration can be 
performed on a block of 20 symbols, which represents two WCDMA time slots. The R-matri- 
5 ces are assumed to be constant over this period. The sign detector in Equation (10) above can 
be replaced by a hyperbolic tangent detector to improve performance under high input BER. A 
hyperbolic tangent detector has a single slope parameter which varies from one iteration to 
another. 

10 The three R-matrices (R[-l], R[0] and R[l]) are each K v x K v in size. Hence, the total 

number of operation per iteration is 6KJ. The computational complexity of the MDFIC algo- 
rithm depends on the total number of virtual users, which in turn depends on the mix of users 
at various spreading factors. For K v = 200 users (e.g. 100 low-rate users), the computation 
requires 240,000 operations. In one embodiment, two iterations are employed which require a 

15 total of 480,000 operations. For real-time applications, these operations must be performed in 
1/1 5 ms or less. Thus, the total processing requirement is 7.2 GOPS. Computational complexity 
is markedly reduced if a threshold parameter is set such that IC is performed only for those 
\y/mj\ below the threshold. If \yj[m]\ is large, there is little doubt as to the sign of b/m], and 
IC need not be performed. The value of the threshold parameter can be variable from stage to 

20 stage. 

C-Matrix Calculation 

As discussed above, the C-matrix elements are utilized to calculate the R-matrices, 
30 which in turn are employed by an MDF Interference Cancellation routine. The C-matrix ele- 
ments can be calculated by utilizing different techniques, as described elsewhere herein. In 
one approach, the C-matrix elements are calculated directly whereas in another approach the 
C-matrix elements are computed from the T-matrix elements, as discussed in detail below and 
illustrated elsewhere herein. 

35 

More particularly, in one method for calculating the C matrix elements, each C-matrix 
element can be calculated as a dot product between the kth user's waveform and the 1th user's 
code stream, each offset by some multipath delay. For this method of calculation, each time a 
user's multipath profile changes, all the C-matrix elements associated with the changed profile 
40 need to be recalculated. A user's profile can change very rapidly, for example, every 100 msec 
or faster, thereby necessitating frequent updates of the C-matrix elements. Such frequent 
updates of the C-matrix elements can give rise to a large amount of overhead associated with 
computations that need to be performed before obtaining each dot product. In fact, obtaining 
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the C-matrix elements by the above approach may require dedicating an entire processor for 
performing the requisite calculations. 

Another approach according to the teachings of the invention for calculating the C- 
5 matrix elements pre-calculates the code correlations up-front when a user is added to the 
system. The calculations are performed over all possible code offsets and can be stored, for 
example, in a large array (e.g., approximately 21 Mbytes in size), herein referred to as the T- 
matrix. This allows updating C-matrix elements when a user's profile changes by extracting the 
appropriate elements from the Gamma matrix and performing minor calculations. Since the T- 
10 matrix elements are calculated for all code offsets, FFT can be effectively employed to speed 
up the calculations. Further, because all code offsets are pre-calculated, rapidly changing mul- 
tipath profiles can be readily accommodated. This approach has a further advantage in that it 
minimizes the use of resources that need to be allocated for extracting the C-matrix elements 
when the number of users accessing system is constant. 

15 

C-matrix Elements Expressed in Terms of Code Correlations 

As discussed above, the R-matrix elements can be given in terms of the C-matrix ele- 
ments as follows: 

20 

p tt [m -\A,A k = X X Re \a lq C lkgq .[m ] 

q=\ q'=\ [ J 



30 



(11) 



where C }kqq ,[m'] is a five-dimensional matrix of code correlations. Both / and k range 
from 1 to K v , where K y is the number of virtual users. The indices q and q ' range from 1 to L, 
representing the number of multipath components, which in this exemplary embodiment is 
35 assumed to be 4. The symbol period offset m ' ranges from -1 to 1 . The total number of matrix 
elements to be calculated is then N C =3(K V L) 2 =3(800) 2 =1.92M complex elements, requir- 
ing 3.84 MB of storage if each element is a byte. The following symmetry property of the C- 
matrix elements can be utilized to halve the storage requirement, for example, in this case to 
1.92 MB: 



40 



N 

C w ^h^] = -^C^[/nl (12) 
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It is evident from the above Equation (12) that each element of C lkqq {m 'J is formed as a 
complex dot product between a code vector c, and a waveform vector s^.. In this exemplary 
embodiment, the length of the code vector is 256. The waveform sjt], herein referred to as the 
signature waveform for the kth virtual user, is generated by applying a pulse-shaping filter g[t] 
5 to the spread code sequence c k [n] as follows: 

where N = 256 and gftj is the raised-cosine pulse shape. Since gftj is a raised-cosine 
1 0 pulse as opposed to a root-raised-cosine pulse, the signature waveform sjt] includes the effects 
of filtering by the matched chip filter. For spreading factors less than 256, some of the chips 
cjp] are zero. The length of the waveform vector sjt] is L + 255N c , where L g is the length of 
the raised-cosine pulse vector gftj and N e is the number of samples per chip. The values for 
these parameters in this exemplary embodiment are selected to be L = 48 and N c = 4. The 
1 5 length of the waveform vector is then 1068, but for performing the dot product, it is accessed 
at a stride of N = 4, which results in an effective length of 267. 

In this exemplary embodiment, the raised-cosine pulse vector g[t] is defined to be non- 
zero from t = -L/2 + 1 :L g /2, with g[0] = 1 . With this definition the waveform s k [t] is non-zero 
20 in a range from t = -L/2 + 1 : L/2+ 255N, 

By combining Equations (11) and (13), the calculation of the C-matrix elements can be 
expressed directly in terms of the user code correlations. These correlations can be calculated 
up front and stored, for example, in SDRAM. The C-matrix elements expressed in terms of the 
30 code correlations T }k [m] are: 
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10 



15 



35 



40 



= 777 X Z + x] ■ c t [n - «] • c,* [n] 

= £ g[rotf c + x] ~ X c; [n] c k [n- m] 



T tt M = 77l c; [n] c k [n- m] 



Since the pulse shape vector g/w/ is of length L g , at most 2L g /N = 24 real macs need to 
be performed to calculate each element C lkqq .[m 'J (the factor of 2 arises because the code cor- 
relations TJm] are complex). For a given x, the method of the invention efficiently calculates 

2Q the range of m values for which g[mN c + x ] is non-zero as described below. The minimum 
value of m is given by /n fflj/r/ i¥ f + x = - LJ2 + 1 , and t is given by x = m f NN c + x /? - x kq ,. If each 
x value is decomposed as x = n N c + p lq , then m mM ~ ceil[ (- x - LJ2 + l)W c ] = -/n'JV- + 
w^. - U(2N c ) + ceil[ (p kq -p lq + 1)/N e ], where ceil[ (p kq -p lq + ] will be either 0 or 1. It 
is convenient to set this value to 0. In order to avoid accessing values outside the allocation for 

30 gM, gM = 0-0 for « = - L g /2: - L g /2 - (TV - 1). All but one of the Nf possible values for ceil[ 
(^.-^+1)/^ ]are0. 

Accordingly, the following relation holds: 



wherein L g is divisible by 2jV , and LJ(2N c ) is a system constant. 

Since, the maximum value of /w is given by rn maxJ N c + x = L/2, the following holds: 

m majl = floor[ (- t + U2)/N ] = -m'N- n lq + V + LJ(2N c ) + floor[ 0y- ]. 
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Further, floor[ ip kq -p l )IN c ] can be either -1 or 0. In this exemplary embodiment, it is 
convenient to set this value to 0. In order to avoid accessing values outside the allocation for 
g[n], gfn] is set to 0.0 {g[n]= 0.0) for n = -L/2 + 1: L/2 + N c . It is noted that half of the N/ 
possible values for floor[ (p^- p ftJ )fN c ] are 0. Accordingly, the following relation holds: 

5 

. = ~ m >N " % + V + L g l(2N € ) (1 6) 

The values of m minl and m majJ are quickly calculable. 

10 The calculation of the C-matrix elements typically requires a small subset of the T 

matrix elements. The T matrix elements can be calculated for all values of m by utilizing Fast 
Fourier Transform (FFT) as described in detail below. 



15 



Using FFT to Calculate the r-matrix Elements 

It was shown above that the r-matrix elements can be represented as a convolution. 
Accordingly, the FFT convolution theorem can be exploited to calculate the r-matrix elements. 
From the above Equation (14), the r-matrix elements are defined as follows: 

1 N ~ ] 

20 r !k [m] = — 5»]cJ/7~m] (17) 

2N, n=0 

where N = 256. Three streams are related by this equation. In order to apply the convo- 
lution theorem, these three streams are defined over the same time interval. The code streams 
cjn] and c/nj are non-zero from n = 0:255. These intervals are based on the maximum spread- 
30 ing factor. For higher data-rate users, the intervals over which the streams are non-zero are 
reduced further. The intervals derived from the highest spreading factor are of particular inter- 
est in defining a common interval for all streams because they represent the largest intervals. 
The common interval allows the FFTs to be reused for all user interactions. 

35 With reference to Figure 15, the range of values of m for which TJm] is non-zero can 

be derived from the above intervals. The maximum value of m is limited by n-m>0, which 
gives 



40 



255 -m^ =0 => =255 (18) 
and the minimum value of m is limited by n~ m ^ 255 ? which gives 
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O-*"™. = 255 =*/»„,„ =-255 

(19) 

To achieve a common interval for all three streams, an interval defined by m = -M/2: 
M/2 -1,M = 512 is selected. The streams are zero-padded to fill up the interval, if needed. 

Accordingly, the DFT and IDFT of the streams are given by the following relations: 



10 



30 



35 



2 



C,[r]= X c,M-e _y2 



2 

it-l 
2 



cM—tc.lrW™"" (20) 
13 M M 



which gives 



20 if- 

I „= 



M_ M_ ] A/_ 

2 2 2 " 



— ^ J y C 1 [r]V 2 " ( "-" ,)r ' w Y C]\r\e ilmr ^ 

2N > M \Zk r Zk 



M_ 
2 



1 



— -I — -1 — -1 



2N.M 2 * l " Af 

— -i 

i *~ 



X Q[r] C;[r]e-^ /A ' (21) 



2 



40 Hence, r^/m/ can be calculated for all values of m by utilizing FFT. Based on the 

analysis presented above, many of these values will be zero for high data rate users. In this 
exemplary embodiment, only the non-zero values are stored in order to conserve storage space. 
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The values of m for which rjmj is non-zero can be determined analytically, as described in 
more detail below and illustrated elsewhere herein. 



Storage and Retrieval ofT-matrix Elements 

5 

As discussed above, the values of the T-matrix elements which are non-zero need to be 
determined for efficient storage of the T-matrix. For high data rate users, certain elements c/nj 
are zero, even within the interval n = 0:N-1, N = 256. These zero values reduce the interval 
over which TJmJ is non-zero. In order to determine the interval for non-zero values consider 
1 0 the following relations: 



(22) 



The index /, for the /th virtual user is defined such that c/nj is non-zero only over the 
15 interval n = j,N, : j,N, +N, -1 . Correspondingly, the vector cjn] is non-zero only over the 
interval n = j k N k :j k N k +N k -l. Given these definitions, TJm] can be rewritten as 



20 



30 



T„ = — £ c]{n + j,N l )c k [n + j,N, - m) 

The minimum value of m for which TJm] is non-zero is 
i-vrt^-JM+jM-Nt+l 

and the maximum value of m for which TJm] is non-zero is 
The total number of non-zero elements is then 



(23) 



(24) 



(25) 



35 



St — + 1 

max 2 min 2 



= N,+N k -\ 



(26) 



The table below provides the number of bytes per l,k virtual-user pair based on 2 bytes 
40 per element - one byte for the real part and one byte for the imaginary part. 





N k = 256 


128 


64 


32 


16 


8 


4 


N, = 256 


1022 


766 


638 


574 


542 


526 


518 
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128 


766 


510 


382 


318 


286 


270 


262 


64 


638 


382 


254 


190 


158 


142 


134 


32 


574 


318 


190 


126 


94 


78 


70 


16 


542 


286 


158 


94 


62 


46 


38 


8 


526 


270 


142 


78 


46 


30 


22 


4 


518 


262 


134 


70 


38 


22 


14 



10 



15 



The memory requirements for storing the T matrix for a given number of users at each 
spreading factor can be determined as described below. For example, for K virtual users at 
spreading factor N g s 2*~ q 9 q = 0:6, where K is the gth element of the vector K (some ele- 
ments of K may be zero), the storage requirement can be computed as follows. Let Table 1 
above be stored in matrix Mwith elements M . For example, M M = 1022, and M QJ = 766. The 
total memory required by the Y matrix in bytes is then given by the following relation 

tyta Z—i I o 99 ^ 9 9 9 



q=0 



q'=q+\ 



•1 



z <7=0 [ ? =0 



(27) 



For example, for 200 virtual users at spreading factor N Q = 256, K q = 2005 0 , which in 
2Q turn results in M byta - V*K 0 (K 0 + l)M m = 1 00(20 1)( 1022) - 20.5 MB. For 10 384 Kbps users , 
K = K 0 b q0 + K 6 b q6 with = 10 and ^ 6 = 640, which results in a storage requirement that is 
given by the following relations: 



30 



M^ = ^(tf^ 5(11)(1022)+ 10(640)(518) + 

320(641)(14) = 6.2 MB. 



35 



The T-matrix data can be addressed, stored, and accessed as described below. In par- 
ticular, for each pair (l t k)> k>= I , there are 1 complex TJm] values for each value of m, where 
m ranges from m min2 to m max2 , and the total number of non-zero elements is m total - rn max2 - m mjn2 
+ 1. Hence, for each pair (l,k), k>= I , there exists 2w ww/ time-contiguous bytes. 



In one embodiment, an array structure is created to access the data, as shown below: 



40 



struct { 

int m_min2; 
int m_max2; 
int m total; 
char * Glk; 
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rithm: 



} G_info[N_VU_MAX][ N_VU_MAX]; 

The C-matrix data can then be retrieved by utilizing the following exemplary algo- 



m min2 = G_info[l][k].m_min2 
m max2 = G_info[l][k].m_max2 
N =L/N 

8 8 c 



Nl =m'*N-L/(2N c ) 
10 form' = 0:1 

forq = 0:L-l 



forq' = 0:L-l 

m mM = N] -\ + v 

15 m =m . ,+N 

max I mini g 

m . = maxf m . , , m . , 1 

mtn L mini 7 mtn2 J 

m = min[ m , , m . ] 

max 1 max/ ' max 2 J 

if m >=m . 

max mtn 

m = m - m . + I 

span max mm 

20 suml=0.0; 

ptrl=&G_info[l][k].Glk[m miri ] 
ptr2 = &gfm mln *N+xJ 
while m > 0 

span 

suml +=( *ptrl++ )*( *ptr2++ ) 
30 m — 

span 

end 



C[m'][l][k][q][q'] = suml 



end 



end 



35 end 
end 



Another method for calculating the T-matrix elements, herein referred to as the direct 
method, performs a direct convolution, for example, by employing the SALzconvx function, to 
40 compute these elements. This direct method is preferable when the vector lengths are small. 
As an illustration of the time required for performing calculations, The table below provides 
exemplary timing data based on a 400 MHz PPC7400 with 16 MHz, 2 MB L2 cache, wherein 
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the data is assumed to be resident in LI cache, 
is not severe. 
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The performance loss for L2 cache resident data 



5 





N, 


Timing {\is) 


GFLOPS 


1024 


4 


19.33 


1.70 


1024 


8 


29.73 


2.20 


1024 


16 


50.55 


2.59 


1024 


32 


92.32 


2.84 


1024 


64 


176.53 


2.97 


1024 


128 


346.80 


3.47 



As discussed above, FFT can also be utilized for calculating the T-matrix elements. 
The time required to perform a 512 complex FFT, with in-place calculation, on a 400 MHz 
PPC7400 with 16 MHz, 2 MB L2 cache is 10.94 (is for LI resident data. Prior to performing 
the final FFT, a complex vector multiplication of length 512 needs to be performed. Exemplary 
timings for this computation are provided in the following table: 

Length | Lo cati on | Timing (us) | GFLOPS 

T023 □ Ztt TT3B 

m& C2 £T27 0^53 

ORAM TO 



Further, exemplary timing data for moving data between memory and the processor is 
provided in the following table: 
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Length 


Location 


Timing (us) 


1024 


L1 


1.20 


1024 


L2 


15.34 


1024 


DRAM 


30.05 



Figure 16 illustrates the T-matrix elements that need to be calculated when a new 
physical user is added to the system. Addition of a new physical user to the system results in 
adding 1 + J virtual users to the systems: that is, 1 control channel + 7 = 256/5Fdata channels. 
The number K v represents the number of initial virtual users. Hence there are (K v + 1 ) elements 
added to the T-matrix as a result of increase in the number of the control channels, and J{K y + 
1) + J(J + l)/2 elements added as a result of increase in the number of the data channels. The 
total number of elements added is then (J + l)[A: v + 1 + J/2]. If FFT is utilized to perform the 
calculations, the total number of FFTs to be performed is (J + 1 ) + (J + 1 )[K v + 1 + J/2]. The 
first term represents the FFTs to transform cjn], and the second term represents the (J+ \)[K v 
+ 1 + J/2] inverse FFTs of FFT{c k [n]}*FFT{c'[n]}. The time to perform the complex 512 
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FFTs can be, for example, 10.94 us, whereas the time to perform the complex vector multiply 
and the complex 512 FFT can be, for example, 24.27/2 + 10.94 = 23.08 ^s. 

In order to provide illustrative examples of processing times, two cases of interest are 
5 considered below. In the first case scenario, a voice user is added to the system while K = 100 
users (K v = 200 virtual users) are accessing the system. Not all of these users are active. The 
control channels are always active, but the data channels have activity factor AF = 0.4. The 
mean number of active virtual users is then K + AF*K = 140. The standard deviation is 
a = y}K AF (\-AF) = 4.90 . Accordingly, there are K v < 140 + 3a < 1 55 active user with a 
10 high probability. 

The second case, which represents a more demanding scenario, arises when a single 
384 Kbps data user is added while a number of users are accessing the system. A single 384 
Kbps data user adds interference equal to (.25 + 0.125*100)/(.25 + 0.400*1) ~= 20 voice users. 
1 5 Hence, the number of voice users accessing the system must be reduced to approximately K = 
100-20-80 (K v = 160). The 3a number of active virtual users is then 80 + (0.125)80 +3(3.0) 
= 99 active virtual users. The reason this scenario is more demanding is that when a single 384 
Kbps data user is added to the system, J + 1 = 64 +1 = 65 virtual users are added to the 
system. 

20 

In the first case scenario in which there are K v = 200 virtual users accessing the system 
and a voice user is added to the system (J = 1), the total time to add the voice user can be (1 + 
1)(10.94 us) + (1 + 1)[200 + 1 + l/2](23.08 jas) = 9.3 ms. 

30 For the second scenario in which there are K v = 160 virtual users accessing the system 

and a 384 Kbps data user is added (J= 64), the total time to add the 384 Kbps user can be (64 
+ 1)( 10.94 |is) + (64 + 1)[1 60 + 1 + 64/2](23.08 us) = 290 ms, which is significantly larger than 
9.3 ms. Hence, at least for high data-rate user, the T-matrix elements are calculated via convo- 
lutions. 



35 
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In the direct method of calculating the T-matrix elements, the SAL zconvx function is 
utilized to perform the following convolution: 

1 N '~ ] 

r„ [m] = — £ c\ [n + j l N l ] • c k [n + j { N t - m] 

= ^T7 X c i I" + Jk N k + w] • cj/i + j h N h ] (28) 



89 



Express Mail Label: u £V 093 931 806 US 



For each value of m, N min = min{Af r N k ) complex macs (cmacs) need to be performed. 
Each cmac requires 8 flops, and there are m toui = N t + N k -\ /n-values to calculate. Hence, the 
total number of flops is &N mjn (N t + N k - 1). In the following, it is assumed that the convolution 
calculation is performed at 1.50 GOPs = 1500 ops/^s. The time required to perform the convo- 
lutions is presented in the table below 





N k = 256 


128 


64 


32 


16 


8 


4 


N=256 


697.69 


261.46 


108.89 


48.98 


23.13 


11.22 


5.53 


128 


261.46 


174.08 


65.19 


27.14 


12.20 


5.76 


2.79 


64 


108.89 


65.19 


43.35 


16.21 


6.74 


3.03 


1.43 


32 


48.98 


27.14 


16.21 


10.75 


4.01 


1.66 


0.75 


16 


23.13 


12.20 


6.74 


4.01 


2.65 


0.98 


0.41 


8 


11.22 


5.76 


3.03 


1.66 


0.98 


0.64 


0.23 


4 


5.53 


2.79 


1.43 


0.75 


0.41 


0.23 


0.15 



10 



15 



The total time to calculate the T-matrix is then given by the following relation: 



20 



« \K(K+\) « 

q=0 I z <? =<7+l 



9=0 [ </-=0 J 



= -[K diag{T) + K T T K] 



(29) 



30 



where T are the elements in the above Table 5. Now suppose K' = K + A, where A a = 



75 + 7 8 , and where x and y are not equal. Then 

x qx y qy> * 



AT r =T r (K)-T r (K) 



40 



(30) 



In the first scenario, there are K v = 200 virtual users accessing the system and a voice 
user is added to the system (7=1). Hence, K = Kb q0 (SF = 256), K y = 200, 7=7=2 and J = 
0. The total time is then 

l AJ{J+ \)T 00 + JKT 0Q = (0.5)(2)(3)(0.70 ms) + (2)(200)(0.70 ms) = 283 ms 
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This number is large enough to require that for voice users, at least, the r-matrix ele- 
mentst be calculated via FFTs. 

For the second scenario, there are K v = 160 virtual users accessing the system and a 384 
5 Kbps data user is added to the system (J = 64). Hence, K = Kb q0 (SF = 256), K v = 160,7 = 1 
(control) and J y = J=64 (data). The total time is then 

(K v + \)T 00 + J(K v + l)T 06 + (J + 1X7/2)^ 

10 = (161)(697.7 ms) + (64)(161)(5.53 |as) + (65)(32)(0.15 [is) = 

1 12.33 ms + 56.98 ms + 0.31 ms - 169.62 ms 

Accordingly, these calculations should also be performed by utilizing FFT, which can 
15 require, for example, 23.08 \is per convolution. In addition, 1 FFT is required to compute 
¥FT{c k '[nJ}) for the single control channel. This can require an additional 10.94 |is. The total 
time, then, to add the 384 Kbps user is 

10.94 (is + (161)(23.08) us + (64)(161)(5.53) us + (65)(32)(0. 1 5) \is = 

20 

= 61.02 mst 

r-matrix elements to SDRAM 

30 With reference to above Equation (27), the size of the r-matrix in bytes is given by the 

following relation: 



35 



,=0 

6 



11 i-i 1 9 91 



2 ,^,+1 
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= - Y | KM aa + T K K M . 

^ Aj\ q qq La q q qq 

= ^[K diag(M)+K T M K] ( 31 ) 
Now suppose K f = K + A, where A q = Jb qx + Jb^ and where x and y are not equal 



Then 



91 



10 



Express Mai) Label: 

u£V 093 931 806 US 



9=0 



Consider a first exemplary scenario in which K - 2008^ (SF = 256) and a single voice 
user is added to the system: J x = 2 (data plus control), and J = 0. The total number of bytes to 
be written to SDRAM is then 0.5(2)(3)(1022) + 200(2)(1022) = 0.412 MB. Assuming a 
SDRAM write speed of 1 33MHz*8 bytes * 0.5 = 532 MB/s, the time required to write T-matrix 
to SDRAM is then 0.774 ms. 



For additional illustration of the time required for storing the T-matrix, consider a 
second scenario in which K q = 1605^ (SF = 256), and a single 384 Kbps (SF = 4) user is added 
to the system \J = \ (control) and J = 64 (data). The total number of bytes i s then 0. 5( 1 )(2)( 1 022) 
+ 0.5(64)(65)(14) + 160(1(1022) + 64(518)} = 5.498 MB. The SDRAM write speed is 
133MHz*8 bytes * 0.5 = 532 MB/s. The time to write to SDRAM is then 10.33 ms. 



Packing the Gamma-Matrix Elements in SDRAM 

In this exemplary embodiment, the maximum total size of the T-matrix is 20.5 MB. If 
it is assumed that in order to pack the matrix, every element must be moved (this is the most 
demanding scenario), then for a SDRAM speed of 133MHz*8 bytes * 0.5 = 532 MB/s, the 
move time is then 2(20.5 MB)/(532 MB/s) = 77.1 ms. If the T-matrix is divided over three 
processors, this time is reduced by a factor of 3. The packing can be done incrementally, so 
there is no strict time limit. 



Extracting Gamma-Matrix Elements from SDRAM 

As described above, in this exemplary embodiment, the C-matrix data is retrieved by 
utilizing the following algorithm: 



40 
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m roin 2 = G_info[l][k].m_min2 
= G_info[l][k].m_max2 



N = L /N 

g g c 



Nl =m'*N-L/(2N) 
5 form' = 0:1 

forq = 0:L-l 

forq' = 0:L-l 

x = mT+T /? - x v 
m . = Nl -n, + n. . 

mm/ tq kq 

10 m =m . , + N 

meal mini g 

m — min[ m , , vn , 1 

max L mar/ 7 max/ J 

if m >= m . 

max nun 

m = m - m . + 1 

span max mm 

15 sum 1=0.0; 

ptrl=&G_info[l][k].Glk[m m J 
ptr2 = &g[m min *N+x] 
while m > 0 

span 

sum] += ( *ptrl++ )*( *ptr2++ ) 
20 int- 
end 

C[m')[l][k][q][q'] = suml 



end 



end 



30 end 
end 

The time requirements for calculating the T-matrix elements in this exemplary embodi- 
ment, when a new user is added to the system was discussed above. The time requirements for 
35 extracting the corresponding C-matrix elements are discussed below. 

The TJm] elements are accessed from SDRAM. It is highly likely that these values 
will not be contained in either LI or L2 cache. For a given (l,k) pair, however, the spread in x 
is likely to be, for most cases, less than 8 ps (i.e. for a 4 ps delay spread), which equates to (8 
40 |is)(4 chips/ps)(2 bytes/chip) = 64 bytes, or 2 cache lines. In an embodiment in which data is 
read in for two values of m \ a total of 4 cache lines must be read. This will require 16 clocks, 
or about 16/133 = 0.12 |is. However, in some embodiments, accesses to SDRAM may be per- 
formed at about 50% efficiency so that the required time is about 0.24 us. 
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If a user / = x is added to the system, the elements Cfm f ][x][k][q][q 'J for all m\k y q 
and q ' need to be fetched. As indicated above, all the m ', q and q ' values are typically contained 
in 4 cache lines. Hence, if there are K v virtual users, 4K v cache lines need to be read, thereby 
requiring 32K v clocks, where the number of clocks has been doubled to account for the 50% 
5 efficiency in accessing the SDRAM. In general, addition of 7+ 1 virtual users to the system at 
a time, requires 32K v (J + 1) clocks. 

In one example where there are 1 55 active virtual users and a new voice user is added 
to the system, the time required to read in the C-matrix elements can be 32(155)(1 + 1) clocks/ 
10 (133 clocks/ps) = 74.6 ps. The present industry standard hold time t h for a voice call is 140 s. 
The average rate X of users added to the system can be determined from Xt h = K, where K is the 
average number of users utilizing the system. For K = 100 users, X = 100/140 s = 1 user are 
added per 1 .4 s. 

1 5 In another example where there are 99 active virtual users and a 384 Kbps user is added 

to the system, the time required to read in the C-matrix elements can be 32(99)(64 + 1) clocks/ 
(133 clocks/ps) = 1.55 ms. However data users presumably will be added to the system more 
infrequently than voice users. 

20 Time to Extract Elements When r y Changes 

Now suppose, for example, that user / = x lag q -y changes. This necessitates fetching 
the elements C[m'][x][k][y][q'J for all m\ k and q\ All the q' values will be contained typi- 
cally in 1 cache line. Hence, 2{K){\) = 2K v cache lines need to be read in, thereby requiring 
30 \6K clocks, where the number of clocks has been doubled to account for the 50% efficiency 
in accessing the SDRAM. In general, when a time lag changes, there are J+ 1 virtual users for 
which the C-matrix elements need to be updated. Such updating of the C-matrix elements can 
require 16/^(7+ 1) clocks. 

35 In one example in which 1 55 active virtual users are present and a voice user's profile 

(one lag) changes, the time required to read in the C-matrix elements can be 16(155)(1 + 1) 
clocks/(133 clocks/|is) = 37.3 ps. As discussed above, for high mobility users, such changes 
should occur at a rate of about 1 per 100 ms per physical user. This equates to about once per 
1 .33 ms processing interval, if there are 100 physical users. Hence, approximately 37.3 ps will 

40 be required every 1.33 ms. 

In another example where there are 99 virtual users and a 384 Kbps data user's profile 
(one lag) changes, the time required to read in the C-matrix elements can be 16(99)(64 + 1) 
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c!ocks/(133 clocks/us) = 0.774 ms. However data users will have lower mobility and hence 
such changes should occur infrequently. 

Writing C-Matrix Elements to L2 Cache 

5 

Consider again the case where user / = jc is added to the system. In such a case, the ele- 
ments Cfm 'Jfxjfkjfqjfq '] for all m\k, q and q ' need to be written to cache. If there are K y 
active virtual users, 4KL 2 bytes need to be written, where the number of bytes have been 
doubled because the elements are complex. In general, addition of J + 1 virtual users to the 
10 system at a time will require 4KL 2 (J + 1) bytes to be written to L2 cache. 

In one example, there are 155 active virtual users and a new voice user is added to the 
system. In this case, the time required to write the C-matrix elements can be 4(155)(16)(1 + 1) 
bytes/(2128 bytes/^is) = 9.3 us. 

15 

In another example, there are 99 active virtual users and a 384 Kbps user is added to the 
system. In such a case, the time required to write the C-matrix elements can be 4(99)(16)(64 + 
1) bytes/(2128 bytes/^s) = 193.5 |as. Data users are typically added to the system more infre- 
quently than voice users. 

20 

Time to Extract Elements When Changes 

Consider a situation in which for user / = x lag q = y changes. In such a case, the ele- 
ments C[m 'Jfxjfkjfqjfq J for all m \ k and q ' need to be written. If there are K v active virtual 
30 users, 4KL bytes need to be written, where the number of the bytes has been doubled since the 
elements are complex. In general, addition of J + 1 virtual users the system at a time will 
require 4KL(J + 1) bytes to be written to L2 cache. 

In one example, there are 1 55 active virtual users and a voice user's profile (one lag) 
35 changes. In such a case, the time required to write the C-matrix elements will be 4(155)(4)(1 
+ 1) bytes/(2128 bytes/jis) = 2.33 us. 

In a second case, there are 99 active virtual users and a 384 Kbps data user's profile 
(one lag) changes. Then, the time required to write the C-matrix elements will be 4(99)(4)(64 
40 + l)bytes/(2128 bytes/fis) = 48.4 \is. However data users will have lower mobility and hence 
such changes typically occur infrequently. 
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Packing C-matrix Elements In L2 Cache 

In this exemplary embodiment, the C-matrix elements are packed in memory every 
time a new user is added to or deleted from the system, and every time a new user becomes 

5 active or inactive. In this embodiment, the size of the C-matrix is 2(3/2)(KL) 2 = 3(KL) 2 bytes. 
If three processors are utilized, the size per processor is (KL) 2 bytes. Hence, the total time 
required for moving the entire matrix within L2 cache is 2{KL) 2 bytes/(2128 bytes/^is), where 
the factor of 2 accounts for read and write. By way of example, if there are 1 55 active virtual 
users, the time required to move the C-matrix elements is 2(155*4) 2 bytes/(2128 bytes/fis) = 

10 0.361 ms, whereas if there are 99 active virtual users the time required to move the C-matrix 
elements is 2(99*4) 2 bytes/(2128 bytes/^s) = 0.147 ms. 

Hardware Calculation Of T-matrix Elements 

15 As discussed above, the C-matrix elements can be represented in terms of the underly- 

ing code correlations in accord with the following relation: 

20 1 

3U m LiSi n 

= XgK + T]T () [m] 



r «W s ^rI c ?W c k [n-m] 

tmT+vv (33) 



35 



The T-matrix represents the correlation between the complex user codes. The complex 
code for user / is assumed to be infinite in length, but with only N f non-zero values. The non- 
zero values are constrained to be ±l±y\ The T-matrix can be represented in terms of the real 
and imaginary parts of the complex user codes as follows: 
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5 



10 



15 



20 



30 



= ^rZK[«]-^;[n]}{cf[/»-m]+yc j ;[«-m]} 

= ^rZK[«]^[/i-m]+c;[n]c 4 '[/i-m] 
+ }c* [n] ■ c[ [n-m]- jcj [n] ■ cf [n - m)} 

= if [m] + r* w + y {r* [»] - 1* w} (34) 

where 



r 2M^lc;w. c ;[«-m] 

rrw-^rl^w-*—] (35) 



Consider any one of the above real correlations, denoted 

Tf [m] = — Y cf [n] • c[ [n-m] 
lk 2N l „ (36) 

whereof and Y can be either R or /. Since the elements of the codes are now constrained 

to be i 1 or 0, the following relation can be defined: 

35 cf[n] = (l-2yf[n])mf[ n ] (37) 

where Y? ["] and mf[n] are both either zero or one. The sequence mf[ri\ is a mask 
used to account for values of cf[n] that are zero. With these definitions, the above Equation 
(4) becomes 

40 
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10 



15 



r^[m}^^J j (l-2yf[n})mf[n](l-2Y k [n-m})ml[n-m] 



2N, 



X(l - 2tf [«]) (1 -2tf [« -m])- m?[n\m r k [n - m) 



= ^l^-^[rt}®i f l [n-m})\mf[n]ml[n-m] 



= 2^^' flWl " 2 ^ [m] ' 



^WaX(TffW®TfI[«-«])-VW-«r[»-«] 



(38) 



where © indicates modulo-2 addition (or logical XOR). 

20 

In addition to configurations discussed elsewhere herein, Figures 17, 18 and 19 illus- 
trate exemplary hardware configurations for computing the functions M™[rri\ and N™[m\ 
for calculating the T-matrix elements. Once the functions M™[rri\ and N™[m\ are obtained, 
the remaining calculations for obtaining the T-matrix elements can be performed in software, 
30 or hardware. In this exemplary embodiment, these remaining calculations are performed in 
software. More particularly, Figure 17 shows a register having an initial configuration subse- 
quent to loading a code and a mask sequences. Further, Figure 18 depicts a logic circuit for 
performing the requisite boolean functions. Figure 19 depicts the configuration of the register 
after implementing a number of shifts. 

35 

The four functions ^[m] corresponding to X, Y = R, J which are components of 
T/Jw] can be calculated in parallel. For K v = 200 virtual users, and assuming that 10% of all 
(/, k) pairs need to be calculated in 2 ms, then for real-time operation, 0.10(200) 2 = 4000 
r tt [/w] elements (all shifts) need to be computed in 2 ms, or about 2M elements (all shifts) per 
40 second. For K v = 1 28 virtual users, the requirement drops to 0.8 1 92M elements (all shifts) per 
second. 
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In this embodiment, the T lk [m] elements are calculated for all 512 shifts. However, not 
all of these shifts are needed- Thus, it is possible to reduce the number of calculations per 
T lk [m] elements by calculating only those elements that are needed. 

As described in more detail elsewhere herein, in one hardware implementation of the 
invention, a single processor is utilized for performing the C-matrix calculations whereas a 
plurality of processors, for example, three processors, are employed for the R-matrix calcula- 
tions, which are considerably more complex. In what follows, a load balancing method is 
described that calculates optimum R-matrix partitioning points in normalized virtual user 
space to provide an equal, and hence balanced, computational load per processor. More par- 
ticularly, it is shown that a closed form recursive solution exists that can be solved for an arbi- 
trary number of processors. 

Balancing Computational Load Among Processors For Parallel Calcu- 
lation Of R-matrix 

As a result of the following symmetry condition, only half of the R-matrix elements 
need to be explicitly calculated: 

R lk (m) = ^R kl (-m). (39) 

In essence, only two matrices need to be calculated. One of these matrices is combina- 
tion of R(l) and R(-l), and the other is the R(0) matrix. In this case, the essential R(0) matrix 
elements have a triangular structure. The number of computations performed to generate the 
raw data for the R(l)/R(-1) and R(0) matrices are combined and optimized as a single number. 
This approach is adopted due to the reuse of the X matrix outer product values (see the above 
Equation (8)) across the two R-matrices. Combining the X matrix and correlation values 
dominate the processor utilization since they represent the bulk of the computations. In this 
embodiment, these computations are employed as a cost metric for determining the optimum 
loading of each processor. 

The optimization problem can be formulated as an equal area problem, where the solu- 
tion results in equal partition areas. Since the major dimensions of the R-matrices are given in 
terms of the number of active virtual users, the solution space for the optimization problems 
can be defined in terms of the number of virtual users per processor. It is clear to those skilled 
in the art that the solution can be applicable to an arbitrary number of virtual users by normal- 
izing the solution space by the number of virtual users. 
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With reference back to Figure 10, the computations of the R(l)/R(-1) matrix can be 
represented by a square HJKM while the computations of the R(0) matrix can be represented 
by a triangle ABC. From elementary geometry, the area of a rectangle of length b and height 
h is given by: 

5 

A r =bh. (40) 
and the area of a triangle with a base width b and a height h is given by 
10 4=\bh. (41) 

Accordingly, a combined area of a rectangle A ri and a triangle A A having a common 
height a. is given by the following relation: 



20 



15 4 = 4, +4, 

1 2 

= a.+-af 

' 2 ' (42) 

wherein A. provides the area of a region below a given partition line. For example, A 2 
provides the area within the rectangle HQRM plus the region within triangle AFG. The differ- 
ence in the area of successive partition regions is employed to form a cost function. More 
particularly, 

B, = A i -A i _ i 

1 2 1 2 

= -Z a i +o,-- a,-.,-*,--, 

2 2 (43) 

35 For an optimum solution, B.'s corresponding to i = 1 , 2, N, where N is the number 

of processors performing the calculations, are equal. Because the total normalized load is 
equal to A N , the load per processor is equal to A^N. That is 

A^ 3 

Bi= T = T = 2N' (44) 



30 



40 



fori = 1,2, ...,N. 
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By combining the above equations for B., the solution for a. can be found by finding the 
roots of the following equation: 

+ I q] J_ = 0 . (45) 

2 4 ' 2 ' 1 11 2W 

Hence, the solution of a. is given as follows: 



a f .=-l±^? 



2 . . 3 (46) 



,+2a, .+— , 



10 The negative roots of the above solution for a. are discarded because the solution space 

falls in a range [0,1]. Although it appears that a solution of a, requires first obtaining values 
of a.- 1, expanding the recursion relations of the a. and utilizing the fact that a Q equals 0 results 
in obtaining the following solution for a, that does not require obtaining a r \: 



15 fl ' = " ,+ J^i (47) 

The table below illustrates the normalized partition values of two, three, and four pro - 
cessors. To calculate the actual partition values, the number of active virtual users is multi - 
plied by the corresponding table entries. Since a fraction of a user can not be allocated, a 
20 ceiling operation can be performed that biases the number of virtual users per processor 
towards the processors whose loading function is less sensitive to perturbations in the number 
of users 



30 



35 



I ,ocation 


Two processors 


Three processors 


Kourprocessors 




-1 = ^(0.5811) 


-1 + V2 (0.4142) 


_1+ /Z (0.3229) 






-1 + V3 (0.7321) 


-1 + J| (0.5811) 


a s 






-1 + J— (0.8028) 



The above methods for calculating the R-matrix elements can be implemented in hard - 
ware and/or software as illustrated elsewhere herein. With reference to Figure 20, in one 
40 embodiment, the above calculations are performed by utilizing a single card that is populated 
with four Power PC 7410 processors. These processors employ the AltiVec SIMD vector 
arithmetic logic unit which includes 32 128-bit vector registers. These registers can hold 
either four 32-bit float, or four 32-bit integers, or eight 1 6-bit shorts, or sixteen 8-bit characters. 
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Two vector SIMD operations can be performed by clock. The clock rate utilized in this 
embodiment is 400 Mz, although other clock rates can also employed. Each processor has 32 
KB of LI cache and 2 MB of 266 MHz L2 cache. Hence, the maximum theoretical perfor- 
mance level of these processors is 3.2 GFLOPS, 6.4 GOPS (16-bit), or 12.8 GOPS (8-bit). In 
5 this exemplary embodiment, a combination of floating-point, 1 6-bit fixed-point and 8-bit fixed- 
point calculations are utilized. 

With continued reference to Figure 20, the calculation of the C-matrix elements are 
performed by a single processor 220. In contrast, the calculation of the R-matrix elements are 
10 divided among three processors 222, 224, and 226. Further, a RACE++ 266 MB/sec 8-port 
switched fabric 228 interconnects the processors. The high bandwidth of the fabric allows 
transfer of large amounts of data with minimal latency so as to provide efficient parallelism of 
the four processors. 

15 



20 



30 



35 



40 
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Vector Processor-Based R-Matrix Generation 

Vector processing is beneficially employed, in one embodiment of the invention, to 
speed calculations performed by the processor card of Figures 2 and 3. Specifically, the 

5 AltiVec™ vector processing resources (and, more particularly, instruction set) of the Motorola 
PowerPC 7400 processor used in node processors 228 are employed to speed calculation of the 
R-matrix. These processors include a single-instruction multiple-data (SIMD) vector arithme- 
tic logic unit which includes 32 128-bit input vector units. These units can hold either four 
32-bit integers, or eight 16-bit integers, or even sixteen 8-bit integers. The clock rate utilized 

10 in this embodiment is 400 Mz, although other clock rates can be also employed. 

Of course, those skilled in the art will appreciate that other vector processing resources 
can be used in addition or instead. These can include SIMD coprocessors or node processors 
based on other chip sets, to name a few. Moreover, those skilled in the art will appreciate that, 
15 while the discussion below focuses on use of vector processing to speed calculation of the R- 
matrix, the techniques described below can be applied to calculating other matrices of the type 
described previously as well, more generally, to other calculations used for purposes of CDMA 
and other communications signal processing. 

20 In the illustrated embodiment, a mapping vector is utilized to create a mapping between 

each physical user and its associated (or "decomposed") virtual users. This vector is populated 
during the decomposition process which, itself, can be accomplished in a conventional manner 
known in the art. The vector is used, for example, during generation of the R-matrix as 
described below. 

30 

As further evident in the discussion below, the X-matrix (see Equation (8)), is arranged 
such that a "strip-mining" method of the boundary elements can be performed to further 
increase speed and throughput. The elements of that matrix are arranged such that successive 
ones of them can be stripped to generate successive elements in the R-matrix. This permits 
35 indices to be incremented rather than calculated. The elements are, moreover, arranged in a 
buffer such that adjacent elements can be multiplied with adjacent element of the C-matrix, 
thereby, limiting the number of required indices to two within the iterative summation loops. 

In the discussion that follows, a node processor 228 operating as a vector processor is 
40 referred to as vector processor 4 1 0. Figure 21 is a block diagram depicting the architecture and 
operation of one such node processor 228, and its corresponding vector processor 4 1 0, used in 
an embodiment of the invention to calculate the R-matrix 428 using integer representations of 
the C-matrix 424 and waveform amplitudes 426. To facilitate a complete understanding of the 
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illustrated embodiments, only a sampling of operands are illustrated, e.g., a few elements each 
of the C-matrix 424, complex amplitudes 426 and R-matrix 428. In actual operation of a 
system according to the invention, the vector processor 4 1 0 can used to process matrices con- 
taining hundreds or thousands of elements. 

5 

As shown in the drawing, the illustrated node processor 228 is configured via software 
instructions to execute a floating point to integer transformation process 406 and an integer to 
floating point transformation process 412, well as to serve as a vector processor 41 0. The rela- 
tionship and signalling between these modes is depicted in the drawing. 

10 

By way of overview, and as discussed above, one or more code-division multiple access 
(CDMA) waveforms or signals transmitted, e.g., from a user cellular phone, modem or other 
CDMA signal source are decomposed into one or more virtual user waveforms. The virtual 
user is deemed to "transmit" a single bit per symbol period of that received CDMA waveform. 
15 In turn, each of the virtual user waveforms is processed according to the methods and systems 
described above. 

In some embodiments, waveform processing is performed using floating-point math, 
e.g., for generating the gamma-matrix, C-matrix, R-matrix, and so on, all in the manner 

20 described above. However, in an embodiment of the invention, e.g., reflected in Figure 21, 
integer math is performed on the vector processor 4 1 0, taking advantage of block-floating point 
representation of the operands. This speeds waveform processing, albeit at the cost of accu- 
racy. However, in the illustrated embodiment, a balance is achieved by through use of 16-bit 
block-floating point representation, e.g., in lieu of conventional 32-bit floating-point represen- 

30 tations. Those skilled in the art will appreciate that block-floating representations of other bit 
widths could be used instead, depending on implementation requirements. 

Referring to Figure 21, the C-matrix 424 is generated by the node processor 228 as 
described above, and is stored in memory accordingly in a floating-point representation, e.g., 

35 C 0 401, Cj 402, and so on. Further, the amplitudes 426 are stored in memory as floating-point 
representations. Both sets of representations are transformed into floating-block format via a 
transformation process 406 which generates a common exponent 414 and a 16-bit integer for 
each operand. The transformation process 416 stores two integers in each word, e.g., C Q 408a, 
Cj 408b, and a Q 409a, a, 409b, and the corresponding block exponent 414. The transformation 

40 process 414 can be performed via special purpose function or through use of extensions to the 
C programming language, as can be seen in a programming listing that is further described. 
The integers stored in memory, e.g., 408, 409, are moved by the transformation process 406 to 
the vector processor 4 1 0 for processing. 
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The vector processor 4 1 0 includes two input vector units 4 1 6, 4 1 8, an output vector unit 
420, and an arithmetic processor 422. Each vector unit is 128-bits in length, hence, each can 
store eight of the 16-bit integer operands. The arithmetic processor 422 has a plurality of oper- 
ating elements, 422a through 422c. Each of the operating element 422a through 422c applies 
5 functionality to a set of operands stored in the input vector units 416, 418, and stores that pro- 
cessed data in the output vector 420. For example, the operating element 422a performs func- 
tionality on operands C 0 4 1 6a and a 0 4 1 8a and generates R Q 420a. The arithmetic processor 422 
can be programed via C programming instructions, or by a field programmable gate array or 
other logic. 

10 

Although vector processor 410 includes two input vector units 416, 418, in other 
embodiments it can have numerous vector units, that can be loaded with additional C-matrix 
and complex amplitude representations at the same time. Further, the operands can be stored in 
a non-sequential order to accommodate increased throughput via storing operands according to 
15 a first-used order. 



20 



As noted above, one way to program the arithmetic processor 422 is through exten- 
sions to a high level programming language. One such program, written in C, suitable for 
instruction the vector processor 422 to generate the R-matrix is as follows: 



#include "mudlib.h" 

#define DO_CALC_STATS 0 

#define DOJTRUNCATE 1 

#define DO_S ATURAT E 1 

#define DO_SQUELCH 0 

fldefine SQUELCH_THRESH 1 . 0 
30 #define TRUNCATE_BIAS 0.0 

#if DO_TRUNCATE 

tfdefilie SATURATEJTHRESH (128.0 + TRUNCATE_BIAS) 
#else 

#define SATURATE_THRESH 127.5 
# end if 

#define SATURATE { f ) \ 
35 { \ 

if ( (f) >= SATURATE_THRESH ) f = (SATURATEJTHRESH - 1.0); \ 
else if ( (f) < -SATURATE_THRESH ) f = - SATURATEJTHRESH ; \ 

} 

#if DO_TRUNCATE 
#if 0 

tfdefine BF8_FIX ( f ) ( (BF8) (FABS(f) <= TRUNCATE_BIAS) ? 0 : \ 

(((f) > 0.0) ? ((f) - TRUNCATE_B I AS ) : \ 
40 ((f) + TRUNCATE_BIAS) ) ) 

tfdefine BF8_FIX( f ) ((BF8)(f)) 
#else 

#defirte BF8_FIX( f ) ( (BF8 )((((( f ) < 0.0)) && ((f) == (float )( (int) (f ))) ) ?\ 

((f) + 1.0) : (f))) 

#endif 
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#else 

#define BF8_FIX ( f ) ((BF8)(((f) >= 0.0) ? ((f)+0.5) : ((£)-0.5))) 
#endif 

tfdefine UPDATE_MAX ( f, max ) \ 

if ( FABS { f ) > max ) max = FABS ( f ) ; 

tfdefine uchar unsigned char 
tfdefine ushort unsigned short 
#define ulong unsigned long 

#if DO_CALC_STATS 

static float max_R_value; 

#endif 

void gen_X_row ( 

C0MPLEX_BF16 *mpathl_bf , 
C0MPLEX_BF16 *mpath2 _bf , 
C0MPLEX_BF16 *X_bf ( 
int phys_index, 
int tot_phys_users 

) ; 

void gen_R_sums ( 

COMPLEX_BF16 *X_bf , 
COMPLEX_BFB *corr_bf , 
uchar *ptov_map, 
BF32 *R_sums, 
int num__phys_users 

) ; 

void gen_R_sums2 ( 

COMPLEX_BF16 *X_bf , 

COMPLEX_BF8 *corra_bf , 

COMPLEX_BF8 * cor rb_bf , 

uchar *ptov_map, 

BF32 *R_sumsa, 

BF32 *R_sumsb, 

int num_j)hys_users 

) ; 

void gen_R_mat rices ( 
BF32 *R_sums, 
float *bf_scalep, 
float *inv_scalep, 
float *scalep, 
BF8 *no_scale_row_bf , 
BF8 *scale_row_bf , 
int num_yirt_users 

) ; 



void mudlib_gen_R ( 

C0MPLEX_BF16 *mpathl_bf , 



BC0MPLEX_BF16 *mpath2_bf , 
C0MPLEX_BF8 *corr_0_bf , 



C0MPLEX_BF8 *corr_l_bf , 



/* ANTENNA DATA 1: TWO AMPLITUDE 
DATA VALUES a hat FOR EACH USER 
*/ 

/* ANTENNA DATA 2 */ 
/* adjusted for starting physical 
user */ 

/* C MATRIX, I.E., C(0), SYMBOL YOU 
ARE ON VERSUS OTHER SYMBOLS */ 

/* adjusted for starting physical 
user */ 

/* C MATRIX. THIS IS A VIRTUAL 
USER BY VIRTUAL USER MATRIX. 
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EACH USER HAS 16 VALUES THAT 
CORRELATE THAT USER TO OTHER 
USERS */ 



uchar *ptov__map, 



/* no more than 256 virts. per phys 



MAPPING OF PHYSICAL TO 
VIRTUAL USERS MAP . IN FURTHER 
EMBODIMENTS , THIS COULD 
DYNAMICALLY CHANGE AS USERS 
ENTER INTO AND LEAVE SYSTEM */ 



float *bf_scalep, 



/* scalar: always a power of 2 */ 



float *inv_scalep / 
float *scalep, 
char *Ll_cachep, 



BF8 *R0_upper_bf , 

BF8 *R0_lower_bf , 

BF8 *Rl_transJbf , 

BF8 *Rlm_bf, 

int tot_phys_users , 

int tot_virt_users , 

int start_j?hys_user , 



int start virt_user, 



/* VECTOR WITH SCALAR FOR EACH 

VIRTUAL USER -- NOTWITHSTANDING 
*/ 

/* start at O'th physical user */ 
/* start at O'th physical user */ 
/* temp: 32K bytes, 32-byte aligned 
*/ 

/* OUTPUTS (BEGINNING AT NEXT LINE) 
*/ 

/+ UPPER PART OF R(l) MATRIX -- A 
TRIANGULAR PACKED MATRIX */ 

/* LOWER PART OF R(0) MATRIX */. 

/* TRANSPOSED FORM OF R{0) */ 

/* R(-l) --> w m" STANDS FOR -1 */ 

/* TOTAL PHYSICAL USERS */ 

/* SUM OF VIRTUAL USERS */ 

/* zero-based starting row 
(inclusive) */ 

/* STARTING PHYSICAL USER TO WHICH 
THIS PROCESSOR IS ASSIGNED */ 

/* relative to start_phys_user */ 



int end_phys_user, 



int end virt user 



/* STARTING VIRTUAL USER TO WHICH 
THIS PROCESSOR IS ASSIGNED +/ 

/* NOTE: THIS IS AN ADVANTAGE 

IN ALLOWING US TO PARTITION A 
GIVEN PHYSICAL USER TO MULTIPLE 
PROCESSORS */ 

/* zero-based ending row inclusive) 
*/ 

/* SAME AS ABOVE, BUT END VALUES */ 
/* relative to end_phys_user */ 



COMPLEX_BF16 *X_bf ; 

BF32 *R_sumsO, *R_sumsl; /* BEGINNING OF PARTITIONING AND 

PARAMETER SET-UP LOGIC +/ 

uchar *R0_ptov_map 

int bump, byte_offset, i, iv, last_virt_user ; 

int R0_align, RO_skipped_virt_users, R0_tcols, R0_virt_users, Rl_tcols; 
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#if DO_CALC_STATS 

max_R_value = 0.0; 
tfendif 

X_bf = (COMPLEX_BF16 *)Ll_cachep; 

byte_offset = tot^Dhys_users * NUM_FINGERS_SQUARED * sizeof (COMPLEX^ 

R sumsO = (BF32 *) (((ulong)X bf + byte offset + R MATRIX ALIGN MASK) & 
- R_MATR I X_AL I GN_MASK ) ; 

byte_offset - tot_virt_users * sizeof (BF32) ; 

R_sumsl = (BF32 * ) ( ( (ulong) R_sums0 + byte_offset + R_MATRIX_ALIGN_MASK) & 
- R_MATR I X_AL I GN_MAS K ) ; 

R0 ptov map = (uchar *) (((ulong)R sumsl + byte offset + R MATRIX ALIGN_ 
~ MASK) & ~ R_MATR I X_AL I GN_MAS K ) ; " 

Rl_tCOls = (tOt_virt_users + R_MATR I X_AL I GN_MAS K ) & ~ R_MATR I X_AL I GN_MAS K ; 

R0 virt users = 0; 

for ( i = start_phys_user; i < tot phys_users ; i++ ) { 

R0_virt_users += ( mt ) ptov_map [ij; 
^ R0jptov_map [i] = ptov_map[i]; 

R0 ptov map [start phys user] -= start_virt user; 

R0 skipped virt_users = tot virt_users - RU_virt_users + start_virt_user ; 
R0_virt_users -= (start_virr_user + 1) ; 

- -inv_scalep; /* predecrement to allow for common 

indexing */ 

for ( i = start_phys_user; i <= endj?hys_user ; i + + ) { /* LOOP OVER ALL 

PHYSICAL USERS 
(ASSIGNED 
TO THIS 
PROCESSOR) */ 

gen_X_row ( /* FIND C CODE THAT PERTAINS TO THIS */ 

mpathl_bf , 
mpath2_bf , 
X_bf , 
i, 

totjphys_users 

>; 

--R0_ptov_map [i] ; /* excludes R0 diagonal */ 

last_virt_user = (i < end_phys_user) ? ( (int) ptov_map [i] - 1) : 

end_virt_user ; 

for ( iv = start_virt_user; (iv + 1) <= last_virt_user ; iv 2 ) { 

gen_R_sums2 ( 

X_bf + {i * NUM_F I NGERS_SQUARED ) , 
corr_0_bf , 

corr_0_bf + < (R0_virt_users - 1) * NUM_FINGERS_SQUARED) , 
R0_ptov_map + i, 

R_sums0 + {R0_skipped_virt_users + 1) , 
R_sumsl + (R0_skipped_virt_users + 1) , 
tot_phys_users - i 

) ; 

R0_tcols = Rl_tcols - (R0_skipped_virt_users & ~R_MATRIX_ALIGN_MASK) ; 
R0_align = (R0_skipped_virt_users & R_MATRIX_ALIGN_MASK) + 1; 

gen_R_matrices ( 

R_sums0 + (R0_skipped_virt_users + 1) , 
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bf_scalep, 

inv_scalep + (R0 skipped_virt_users + 1) , 
scalep + (RO_skipped_virt_users + 1), 
R0_lower_bf + R0_align, 
R0_upperjbf + R0_align, 
R0_virt_users 

); 

R0_upperjbf [ R0_align - 1 ] = 0; /* zero diagonal element */ 

R0_lower_bf += R0_tcols; 
R0_upper_bf += R0_tcols; 

R0_tcols = Rl_tcols - { (RO_skipped_virt_users + 1) & 

~ R_MATR I X_AL I GN_MAS K ) ; 
R0_align = ( (RO_skipped_virt_users + 1) & R_MATRIX_ALIGN_MASK) + 1; 

gen_R_ma trices ( 

R_sumsl + (RO_skipped_virt_users + 2) , 
bf_scalep, 

inv_scalep + (RO_skipped_virt_users + 2) , 
scalep + (RO_skipped_virt_users + 2) , 
R0_lower_bf + R0_align, 
R0_upper_bf + R0_align, 
R0_virt_users - 1 

) ; 

RO_upper_bf [ R0_align - 1 ] = 0; /* zero diagonal element */ 

R0_lower_bf += R0_tcols; 
R0_upper_bf += R0_tcols; 



* create ptov_map[i] number of 32-element dot products involving 

* X bf [i] and corr_l_bf [i] [j] where 0 < j < ptov_map[i] 

*/ " 
gen_R_sums2 ( 
X_bf , 

corr_l_bf , 

corr_l_bf + (tot_virt_users * NUM_FINGERS_SQUARED) , 

ptovjnnap, 

R_sums0 , 

R_sumsl , 

tot_phys_users 

); 



* scale the results and create two output rows (1 per matrix) 
*/ 

gen_R_mat rices ( 
R_sums0 , 
bf_scalep, 

inv_scalep + (RO_skipped_virt_users + 1) , 
scalep, 
Rl_trans_bf , 
Rlm_bf , 

tot_virt_users 

) ; 

Rl_trans_bf += Rl_tcols; 
Rlm_bf += Rl_tcols; 

gen_R_mat rices ( 
R_sumsl , 
bf_scalep, 

inv_scalep + (RO_skipped_virt_users + 2), 

scalep, 

Rl trans_bf, 
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Rlm_bf , 

tot_virt_users 

) ; 

Rl_trans_bf += Rl_tcols; 
Rlm__bf += Rl_tcols; 

corr_0_bf += ( ( (2 * R0_virt 
corr_l_bf += ( (2 * tot_virt 
RO_ptov_map [i] - = 2; 
R0_virt_users - = 2; 
RO_skipped_virt_users += 2 ,- 



f ( iv <= last_virt_user ) { 

bump = R0_ptov_map [ i ] ? 0 : 1; 
gen_R_sums { 

X_bf + ( (i + bump) * NUM_FINGERS_SQUARED) , 

corr_0_bf , 

R0_ptov_map + i + bump, 

R_sums0 + (RO_skipped_virt_users + 1) , 

tot_phys_users - i - bump 

) ; 

R0_tcols = Rl_tcols - <RO_skipped_virt_users & ~ R_MATR I X_AL I GN_MAS K ) ; 
R0_align = (RO_skipped_virt_users & R_MATRIX_ALIGN_MASK) + 1; 

gen_R_mat rices ( 

R_sums0 + (RO_skipped_virt_users + 1) , 
bf_scalep, 

inv_scalep + (RO_skipped_virt_users + 1) , 
scalep + (RO_skipped_virt_users + 1) , 
R0_lower_bf + R0_align, 
R0_upperjbf + R0_align, 
RO_virt_users 

) ; 

R0_upper_bf [ R0_align - 1 ] = 0; /* zero diagonal element */ 

R0_lower_bf += R0_tcols; 
R0_upper_bf += R0_tcols; 

/* 

* create ptov_map[i] number of 32-element dot products involving 

* X_bf[i] and corr_l_bf [i] [ j ] where 0 < j < ptov_map[i] 
*/ 

gen_R_sums ( 
X_bf , 

corr_l_bf , 
ptov_map, 
R_sums0 , 
tot_phys_users 

)/ 
/* 

* scale the results and create two output rows (1 per matrix) 
*/ 

gen_R_mat rices ( 
R_sums0 , 
bf_scalep, 

inv_scalep + (RO_skipped_yirt_users + 1) , 
scalep, 
Rl_trans_bf , 
Rlm_bf , 

tot_virt_users 

) ; 



users) - 1) * NUM_FINGERS_SQUARED) ; 
users) * NUM_FINGERS_SQUARED) ; 
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Rl_trans_bf += Rl_tcols; 
Rlm_bf += Rl_tcols; 

corr_0_bf += (R0_virt_users * NUM_FINGERS — SQUARED) ; 
corr_l_bf += (tot_virt_users * NUM_FINGERS_SQUARED) ; 
R0_ ptov_map [ i 3 - = 1 ; 
RO_virt_users -= 1; 
RO_skipped_virt_users += 1; 

start virt user =0; /* for all subsequent passes */ 

} " ~ 

#if DO_CALC_STATS 

printf( "max_R_value = %f\n", max_R__value ); 

if ( max_R_value > 127.0 ) 

printf ( ****** OVERFLOW *****\n" ); 
#endif 
} 

#if COMPILE_C 

/* OUTPUT PRODUCT OF TWO 
ANTENNAS */ 



void gen_X_row ( 



C0MPLEX_BF16 *mpathl_bf , 
COMPLEX_BF16 *mpath2_bf , 

COMPLEX_BF16 *X_bf , 
int phys_index, 
int tot_phys_users 



/* EACH ANTENNA HAS TWO 

VALUES PER PHYSICAL USER 
*/ 

/* 2ND ANTENNA IS DIVERSITY ANTENNA */ 
/* RESULTING OUTPUT PRODUCT IS REP' D BY 
X sub l,k */ 



COMPLEX_BFl 6 * injnpathlp , 
COMPLEX_BFl 6 *out_mpa thlp , 
int i, j, q, ql ; 
BF32 sir, sli, s2r, s2i; 
BF32 air, ali, a2r, a2i; 
BF32 cr, ci; 



^in_mpath2p; 
*out_mpath2p; 



out_mpathlp = mpathl_bf + (phys_index * NUM_FINGERS) ; 
out_mpath2p = mpath2_bf + (phys_index * NUM_F I NGERS ) ; 

for ( i = 0; i < tot_phys_users; i++ ) { 

in_mpathlp « mpathl_bf + (i * NUM_FINGERS) ; /* 4 complex values */ 
in_rnpath2p = mpath2_bf + (i * NUM_FINGERS) ; /* 4 complex values */ 

j = 0; 

for ( ql = 0; ql < NUM_FINGERS; ql++ ) { 

sir = (BF32) out_mpathlp [ql] .real; 

sli = (BF32)out_mpathlp[qlj .imag; 

s2r = (BF32)out_mpath2p[ql] .real; 

s2i = (BF32) out_mpath2p [ql] .imag; 

for ( q = 0; q < NUM_FINGERS; q++ ) { 

air = (BF32) in_mpathlp[q) .real; 

ali = (BF32) in_mpathlp [q] . imag; 

a2r = (BF3 2) in_mpath2p[q] .real; 

a2i = (BF32) in_mpath2p(q] .imag; 
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BLOCK 



cr = (air * sir) + (ali * sli) ; /* COMBO OF TWO ANTENNAS 

COULD BE MORE, OF COURSE 
*/ 

ci = (air * sli) - (ali * sir); /* cr IS REAL PART OF 

ELEMENT OF X-MATRIX */ 

cr += (a2r * s2r) + (a2i * s2i) ; 
ci += (a2r * s2i) - (a2i * s2r) ; 

X_bf(i * NUM_FINGERS_SQUARED + jj.real = (BF16) (cr >> 16); /* 

X_bf [i * NUM_FINGERS_SQUARED + j ] . imag = (BF16) (ci » 16); 
++ j ; 



void gen_R_sums ( 

COMPLEX_BF16 *X_bf , 
COMPLEX_BF8 *corr_bf , 
uchar *ptov_map, 
BF32 *R_sums, 
int num_phys_users 

) 



{ 

int i, j, k; 
BF32 sum; 

for ( i = 0; i < num_j>hys_users; i++ ) { 
for ( j = 0; j < (int)ptov_map[i] ; j++ ) { 
sum = 0; 

for ( k = 0; k < 16; k+ + ) { 

sum += (BF32)X_bf [k] .real * (BF32) corr_bf - >real ; 
sum += (BF32)X_bf [k] .imag * (BF3 2) corrjbf - >imag ; 
++corr bf; 

} 

*R sums++ = sum; 
} " 

X bf += NUM FINGERS_SQUARED; 



void gen_R__sums2 ( 

COMPLEX_BF16 *X_bf , 

COMPLEX_BF8 *corra_bf , 

COMPLEX_BF8 *corrb_bf , 

uchar *ptov_map, 

BF3 2 *R_sumsa, 

BF32 *R_sumsb / 

int numjphys_users 

) 



{ 

int i, j, k; 
BF32 suma, sumb; 

for ( i = 0; i < num_j>hys_users ; i + + ) { 
for ( j = 0; j < (int)ptovjnap [i] ; j++ ) { 
suma = 0 ; 
sumb = 0; 

for ( k = 0; k < 16; k++ ) { 

suma += <BF32)X_bf [k] .real * (BF32) corra_bf - >real ; 
suma += (BF32)X_bf tk] .imag * (BF32) corrajbf ->imag; 
sumb += (BF32)X_bf [k] .real * (BF32) corrb_bf ->real; 
sumb += (BF3 2)X_bf [k] .imag * (BF32) corrb_bf - >imag; 
++corra_bf ; 
++corrb_bf ; 
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} 



} 

*R_sumsa++ = suma; 
*R sumsb++ = sumb; 

} " 

X bf += NUM FINGERS SQUARED; 
} " 



void gen_R_mat rices ( 
BF32 *R_sums, 
float *bf_scalep, 
floa t * inv_s c a 1 ep , 
float *scalep, 
BF8 *no_scale_row_bf , 
BF8 *scale_row_bf , 
int num_virt_users 

io ) 
{ 

int i ; 

float bf_scale, fsum, fsum_scale, inv_scale, scale; 

bf_scale = *bf_scalep; 
inv_scale = *inv_scalep; 

J5 for ( i = 0; i < num_virt_users; i++ ) { 
scale = scalep[i]; 
fsum = (float) (R_sums [i] ) ; 
fsum *= bf_scale; 

fsum_scale = fsum * inv_scale; 
fsum_scale *■ scale; 

#if DO_CALC_STATS 
20 UPDATE_MAX{ fsum_scale, max_R_value ) 

UPDATE_MAX( fsum, max_R_value ) 
# end if 

#if DO_SQUELCH 

if ( FABS ( fsum_scale ) <= SQUELCH JTHRESH ) fsum_scale = 0.0; 
if ( FABS ( fsum ) <= SQUELCH_THRESH ) fsum - 0.0; 
#endif 

30 #if DO_SATURATE 

SATURATE ( fsum_scale ) 
SATURATE ( fsum ) 
#endif 

no_scale_row_bf [i] = BF8_FIX ( fsum ); 
scale row bf[i] = BF8 FIX ( fsum scale ); 



35 



#endif /* COMPILE_C */ 



A transformation process 412 transforms the block-floating representations stored in 
40 the output vector unit 420 into floating point representations and stores those to a memory 428. 
Here, the R-matrix elements stored in the output vector unit 420 are transformed into floating- 
point representations, which can then be used in the manner described above for estimating 
symbols in the physical user waveforms. 
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In summary, sufficient throughput can be achieved with necessary accuracy using a 
vector processor 410 applying integer math on 16-bit block-floating integers. Of course, in 
other embodiments, different block-floating sizes can be used depending on such criteria as the 
number of users, speed of the processors, and necessary accuracy of the symbol estimates, to 
5 name a few. Further, like methods and logic described can be used to generate other matrices 
(e.g., the gamma-matrix and the C-matrix) and to perform other calculations within the illus- 
trated embodiment. 

A further understanding of the operation of the illustrated and other embodiments of the 
1 0 invention may be attained by reference to (i) US Provisional Application Serial No. 60/275,846 
filed March 14, 2001, entitled "Improved Wireless Communications Systems and Methods"; 
(ii) US Provisional Application Serial No. 60/289,600 filed May 7, 2001, entitled "Improved 
Wireless Communications Systems and Methods Using Long-Code Multi-User -Detection"' 
and (iii) US Provisional Application Serial Number. 60/295,060 filed June 1, 2001 entitled 
15 "Improved Wireless Communications Systems and Methods for a Communications Com- 
puter," the teachings all of which are incorporated herein by reference, and a copy of the latter 
of which may be filed herewith. 

The above embodiments are presented for illustrative purposes only. Those skilled in 
20 the art will appreciate that various modifications can be made to these embodiments without 
departing from the scope of the present invention. For example, the processors could be of 
makes and manufactures and/or the boards can be of other physical designs, layouts or archi- 
tectures. Moreover, the FPGAs and other logic devices can be software or vice versa. More- 
over, it will be appreciated that while the illustrated embodiments decomposes physical user 
30 waveforms to virtual user waveforms, the mechanisms described herein can be applied, as 
well, without such decomposition, and that, accordingly, the terms "waveform" or "user wave- 
form" should be treated as referring to either physical or virtual waveforms unless otherwise 
evident from context. 

35 Therefore, in view of the foregoing, what we claim is: 



40 
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Wireless Communications Systems And Methods For 
Contiguously Addressable Memory Enabled Multiple 
Processor Based Multiple User Detection 

A communications device for detecting user transmitted symbols encoded in spread 
spectrum waveforms (hereinafter "user waveforms") comprising 

a first memory, 

a set of one or more first processing elements, coupled to the first memory, the set of 
first processing elements generating a matrix (hereinafter "R-matrix") representative of 
cross correlations among user waveforms and storing that R-matrix to contiguous loca- 
tions within the first memory. 

The device of claim 1, comprising 

a second processing element coupled with the first memory, 

the second processing element accessing the R-matrix from contiguous locations within 
the first memory and generating symbol estimates as a composition of the R-matrix. 

The device of claim 1 , comprising 

a second memory coupled with the set of first processing elements and a third process- 
ing element, 

the third processing element generating a matrix (hereinafter "gamma-matrix") repre- 
sentative of a correlation between a code associated with one user and those associated 
with one or more other users, 

the third processing element places the gamma-matrix in the second memory. 
The device of claim 3, comprising 

the third processing element generating the gamma-matrix and placing that matrix in 
contiguous location within the second memory, 
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the set of first processing elements accessing the gamma-matrix from contiguous loca- 
tions within the second memory and generating the R-matrix. 

5. The device of claim 3, comprising 

a multi-port switch coupled to the third processing element and to the second memory, 

the third processing element places the gamma-matrix in the second memory via the 
data switch. 

6. The device of claim 3, wherein the gamma-matrix is a composition of a complex con- 
jugate of the code associated with one user and a complex conjugate of the codes asso- 
ciated with one or more other users. 

7. The device of claim 3, wherein the third processing element updates the gamma-matrix 
as users are added or removed from the spread spectrum system. 

8. The device of claim 3, wherein the set of first processing elements generate the R- 
matrix as a composition of the gamma-matrix. 

9. The device of claim 1 , comprising 

a host controller coupled to each of the set of first processing elements, the host control- 
ler generates a partitioning of the R-matrix, that partitioning divides the R-matrix into 
one or more portions based on a number of users and a number of available processing 
elements, 

the host controller assigns to each first processing element a portion of the R-matrix to 
generate according to the partitioning, 

each of the first processing elements generating the assigned portion of the R-matrix 
according to the partitioning, 

the host re-calculates the partitioning of the R-matrix when a user is added or removed 
from the spread spectrum system, and assigns a new portion of the R-matrix to each 
first processing element according to that new partitioning. 
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10. The device of claim 9, wherein each first processing element places its respective por- 
tion of the R-matrix in the first memory according to its respective partition such that 
each portion of the R-matrix is contiguous with respect to the other portions. 

11. A communications device for detecting user transmitted symbols encoded in spread 
spectrum waveforms (hereinafter "user waveforms") comprising 

a first memory, 

a set of one or more first processing elements, coupled to a direct memory access engine 
(hereinafter "DMA engine"), the set of first processing elements generating a matrix 
(hereinafter "R-matrix") representative of cross correlations among user waveforms, 

the DMA engine coupled with the first memory, the DMA engine storing that R-matrix 
to contiguous locations within the first memory. 

1 2. The device of claim 1 1 , comprising 

a second processing element coupled with the first memory, 

the second processing element accessing the R-matrix from contiguous locations within 
the first memory and generating symbol estimates as a composition of the R-matrix. 

13. The device of claim 11 comprising 

a second memory coupled with the set of first processing elements and a third process- 
ing element, 

the third processing element generating a matrix (hereinafter "gamma-matrix") repre- 
sentative of a correlation between a code associated with one user and those associated 
with one or more other users, 

the third processing element places the gamma-matrix in the second memory. 

1 4. The device of claim 1 1 , comprising 

the third processing element generating the gamma-matrix and placing that matrix in 
contiguous location within the second memory, 
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the set of first processing elements accessing the gamma-matrix from contiguous loca- 
tions within the second memory and generating the R-matrix. 

15. The device of claim 1 3, comprising 

a multi-port switch coupled to the third processing element and to the second memory, 

the third processing element places the gamma-matrix in the second memory via the 
data switch. 

1 6. The device of claim 1 3, wherein the gamma-matrix is a composition of a complex con- 
jugate of the code associated with one user and a complex conjugate of the codes asso- 
ciated with one or more other users. 

17. The device of claim 13, wherein the third processing element updates the gamma- 
matrix as users are added or removed from the spread spectrum system. 

18. The device of claim 13, wherein the set of first processing elements generate the R- 
matrix as a composition of the gamma-matrix. 

1 9. The device of claim 1 1 , comprising 

a host controller coupled to each of the set of first processing elements, the host control- 
ler generating a partitioning of the R-matrix, that partitioning divides the R-matrix into 
one or more portions based on a number of users and a number of available processing 
elements, 

the host controller assigns to each first processing element a portion of the R-matrix to 
generate according to the partitioning, 

each first processing element generating the assigned portion of the R-matrix according 
to the partitioning, 

the host re-calculating the partitioning of the R-matrix when a user is added or removed 
from the spread spectrum system, and assigning a new portion of the R-matrix to each 
first processing element according to that new partitioning. 
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20. The device of claim 1 9, wherein DMA engine places each portion of the R-matrix in the 
first memory according to the partitioning such that the each portion of the R-matrix is 
contiguous with respect to the adjacent portions. 

21. A communications device for detecting user transmitted symbols encoded in spread 
spectrum waveforms (hereinafter "user waveforms") comprising 

a first memory, 

a set of one or more first processing elements, coupled to the first memory, the set of 
first processing elements generating a matrix (hereinafter "R-matrix") representative of 
cross correlations among user waveforms and storing that R-matrix to contiguous loca- 
tions within the first memory 

a second processing element coupled to the first memory, the second processing ele- 
ment accessing R-matrix from contiguous location within the first memory and generat- 
ing symbol estimates as a composition of the R-matrix. 

22. The device of claim 2 1 comprising 

a second memory coupled with the set of first processing elements and a third process- 
ing element, 

the third processing element generating a matrix (hereinafter "gamma-matrix") repre- 
sentative of a correlation between a code associated with one user and those associated 
with one or more other users, 

the third processing element places the gamma-matrix in the second memory. 

23. The device of claim 22, comprising 

the third processing element generating the gamma-matrix and placing that matrix in 
contiguous location within the second memory, 

the set of first processing elements accessing the gamma-matrix from contiguous loca- 
tions within the second memory and generating the R-matrix. 



24. 



The device of claim 22, comprising 
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a multi-port switch coupled to the third processing element and to the second memory, 

the third processing element places the gamma-matrix in the second memory via the 
data switch. 

25. The device of claim 22, wherein the gamma-matrix is a composition of a complex con- 
jugate of the code associated with one user and a complex conjugate of the codes asso- 
ciated with one or more other users. 

26. The device of claim 22, wherein the third processing element updates the gamma- 
matrix as users are added or removed from the spread spectrum system. 

27. The device of claim 22, wherein the set of first processing elements generate the R- 
matrix as a composition of the gamma-matrix. 

28. The device of claim 2 1 , comprising 

a direct memory access engine (hereinafter "DMA engine") coupled with the set of 
second processing elements and the first memory, 

the DMA engine placing the R-matrix in contiguous locations within the first memory. 

29. The device of claim 2 1 , comprising 

a direct memory access engine (hereinafter "DMA engine") coupled with the third pro- 
cessing element and the second memory, 

the DMA engine placing the gamma-matrix in contiguous locations within the second- 
memory. 
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Abstract 

The invention provides methods and apparatus for multiple user detection (MUD) pro- 
cessing that have application, for example, in improving the capacity CDMA and other wire- 
less base stations. One aspect of the invention provides a multiprocessor, multiuser detection 
system for detecting user transmitted symbols in CDMA short-code spectrum waveforms. A 
first processing element generates a matrix (hereinafter, "gamma matrix") that represents a cor- 
relation between a short-code associated with one user and those associated with one or more 
other users. A set of second processing elements generates, e.g., from the gamma matrix, a 
matrix (hereinafter, "R-matrix") that represents cross-correlations among user waveforms 
based on their amplitudes and time lags. A third procesing element produces estimates of the 
user transmitted symbols as a function of the R-matrix. 
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